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Polycrystalline  CuInSej  (CIS)  and  related  materials,  are  promising  candidates  of 
low-cost  high  efficiency  solar  cells.  This  thesis  describes  the  applications  of  phase 
diagrams  to  suggest  alternative  absorber  precursor  structures,  explores  the  use  of  Rapid 
Thermal  Processing  (RTP),  and  elucidates  the  structure  of  CuInSej  and  related  materials. 

A critical  assessment  of  the  thermochcmical  and  phase  diagram  data  for  the  Sc 

an  evaluated  CuiSe-InjSej  pseudobinary  phase  diagram  and  suggest  liquidus  lines  in  the 
Cu-ln-Se  ternary  system.  The  resulting  information  was  used  to  propose  a novel  bi-layer 
precursor  structure  consisting  of  a co-deposited  In-Se  (50  to  55  Se  at  %)  top  layer  and  a 
Cu-Se  (50  to  55  Se,  al.%)  bottom  layer  This  precursor  structure  was  designed  to  form 
CIS  by  reacting  a liquid  CuSe  with  a solid  InSe  phase  in  the  process  temperature  window 


524"C  to  61  rc.  Two  different  precursor  structures  were  fabricated  and  tested-  The  first 
one  consists  of  the  three  layers,  ltuSes/CulnSej/CujSe,  and  the  second  one  was  a simple 
two-layer  structure,  InSe/CuSe.  After  RTF,  the  first  precursor  structure  underwent  a 
liquid-solid  reaction  and  formed  large  grained  CulnSej  films.  It  is  believed  the  melting  of 
top  ItuSe.  layer  occurs  through  a monotetic  reaction.  The  second  precursor  structure  also 
lead  to  the  formation  of  large  grained  CIS(-5|im)  films,  as  suggested  in  the  reaction 
pathway  analysis. 

X-ray  Absorption  Fine  Structure  (XAFS)  measurements  were  performed  on 
crystals  of  CuInSei  with  different  compositions  close  to  the  1:1 :2  stoichiometry,  and 
CultliSe.  It  was  found  that  the  Cu-Sc  bond  length  is  independent  of  the  Cu  occupation 
number  in  this  composition  range  within  a length  of  0.01  A,  Least-square  fitting  of  the  Se- 
K edge  EXAFS  spectra  indicated  CulnSe?  consists  of  Se-centercd  tetrahedron  with  an 
average  of  2 Cu  and  2 In  nearest  neighbors.  On  the  other  hand,  CuInjSes  consists  of  Se- 
centered  tetrahedron  with  an  average  of  0.8  Cu  and  2.4  In  nearest  neighbors.  The  Cu-Se 
and  In-Se  bond  lengths  were  determined  to  be  2.424  and  2.598  A,  respectively.  This  result 
directly  proves  that  CuIrljScj  is  indeed  a defect  tetrahedral  structure  and  suggests  it 
consists  of  three  types  of  local  tetrahedral  cationic  clusters  around  each  Se:  20%  2Cu+2ln 
(k=8).  and  40%  0+Cu+2In  (k= 7),  and  40%  0+3  In  (k= 9),  where  □ denotes  vacancy  and  k 


denotes  the  sum  of  the  cation  valence  electrons 


CHAPTER  I 
INTRODUCTION 

The  search  for  a cheap,  clean  and,  renewable  energy  source  has  been  a 
fundamental  issue  for  mankind.  The  major  energy  source  today  is  that  generated  by 
burning  fossil  fuels.  These  valuable  resources  of  the  earth  are  limited  and  will  be  depleted 
in  the  relatively  near  future.  In  addition  to  their  limited  supply,  fossil  fuels  produce 
combustion  gases  that  trap  heat  and  are  the  largest  contributors  to  the  global  warming 
effect  Various  renewable  energy  sources  do  exist  today,  including  photovoltaic  solar 
cells,  bio-fuels,  hydroelectric  power,  and  wind  power.  The  conversion  of  sunlight  directly 
into  electricity  using  the  photovoltaic  properties  of  suitable  materials  is  the  most  elegant 

The  photovoltaic  effect  was  discovered  by  Becquerel  in  1839.  The  first  practical 
photovoltaic  device  was  the  silicon  solar  cell  developed  at  Bell  Laboratories  in  1954  for 
space  application  The  energy  crisis  in  the  early  1970s  raised  the  serious  consideration  for 
photovoltaic  solar  cells  as  an  alternative  energy  source  for  terrestrial  applications. 
Increased  research  efforts  resulted  in  further  improvements  and  the  development  of  new 
photovoltaic  materials  and  devices  These  research  efforts  focused  on  improving 
performance,  lowering  cost,  and  increasing  reliability.  The  last  few  decades  of  research 
and  development  have  brought  photovoltaic  technology  to  a high  level  of  performance 
These  efforts  have  been  reflected  in  a steady  growth  at  1 5 to  20%  per  annum  for  more 
than  a decade  to  reach  SI  billion  sales  in  1996.  The  industiy  expanded  worldwide  by  a 


38%  growth  in  1997.  Despite  the  expanding  terrestrial  market  and  the  progress  in  cell 
technology,  serious  challenges  remain.  For  photovoltaics  to  be  widely  used  the  cost  must 
be  competitive  with  conventional  energy  sources,  where  the  key  issue  is  development  of 
low-cost  manufacturing  processes. 

The  low  energy  density  of  the  incoming  sunlight  requires  large  areas  of 
semiconductors  for  the  production  of  considerable  amounts  of  electricity  Thin  film 
photovoltaics  were  developed  as  a means  of  substantially  reducing  the  cost  of 
photovoltaic  systems  due  to  their  reduced  material  costs,  energy  costs,  handling  costs, 
and  capital  costs.  Many  compound  semiconductor  materials  have  been  examined  for 
application  in  thin  polycrystalline  film  solar  cells  These  materials  include  CujS,  CulnSej 
and  related  alloys,  CdTe,  WScj,  ZnjPi,  and  ZnSiAs.  Among  these  materials  CdTe  and 
CulnSe;  and  related  alloys  have  made  significant  technical  progress.  Thin  film  CdTe  and 
CulnSej  based  solar  cells  have  reported  efficiencies  of  over  1 5%  and  1 8.8%.  respectively, 
for  small  area  laboratory  devices,  and  over  12%  Ibr  prototype  module  Many  companies 
world  wide  are  now  pursuing  commercialization  efforts  based  on  these  two  technologies. 

I I Photovoltaic  Devices 

In  a photovoltaic  device,  electron-hole  (e-h)  pairs  are  generated  by  absorption  of 
photons.  The  e-h  pairs  are  then  separated  by  the  electric  field  present  in  the  cell.  This 
electric  field  can  be  produced  by  using  a p-n  junction,  Schottky  barrier,  MIS  (mctal- 
insulator-semiconductor)  structure,  or  semiconductor-electrolyte  junction.  Most  of  the 
photovoltaic  devices  are  p-n  junction  solar  cells,  including  CuInSct-based  photovoltaic 


The  device  I-V  characteristics  can  be  described  by  the  Poisson  and  current 
continuity  equations: 

V»D  = q(/>-n  + Wj  + A^)  (|.j) 

V.J„  = q(«-G)  (12) 

V.J ,,  =q(G-«)  (l3) 

where  D is  the  displacement  vector,  n and  p are  the  tree  electron  and  hole  densities,  N'a 
and  N [>  denote  the  ionized  donor  and  acceptor  densities,  J„  and  Jp  are  electron  and  hole 
current  densities,  G is  the  free  carrier  generation  rate,  and  RP  and  Kn  are  the  hole  and 
electron  recombination  rates.  Giving  appropriate  boundary  conditions,  these  equations 
can  be  solved  numerically  [Gra  96)  or  analytically  [Fon  80]  for  small  device  structures. 
The  fundamental  material  parameters  associated  with  these  carrier  transport 
equations,  and  thus  the  device  performance,  includes  the  band  gap.  electron  affinity, 
optical  absorption  and  reflection  coefficients,  refraction  index,  dielectric  constants, 
carrier  mobility,  carrier  lifetime,  and  doping  density  (intrinsic  and  extrinsic).  An 
understanding  of  the  correlation  between  the  structural,  electrical,  and  optical  properties 
in  the  constituent  materials  and  the  processing  chemistries  is  essential  in  reproducibly 
manufacturing  high-efliciency  and  low  cost  photovoltaic  modules. 

1.2  A Review  ofCulnSe;-based  Photovoltaic  Technology 

The  first  CuInSei-based  photovoltaic  device  was  reported  by  Shay, 
Wagner  and  Kasper  [Shn  75]  in  1975.  It  was  made  by  evaporating  CdS  on  the  selenium 
(112)  plane  of  a single  crystal  of  p-CulnSe2.  Subsequently  then  thin  film  polycrystalline 
devices  have  been  developed  and  the  first  efficient  thin  film  CdS/CulnSej  solar  cell  was 
fabricated  by  Mickelsen  and  Chen  [Mic  81],  Remarkable  progress  has  been  made  in  the 


past  few  years  by  incorporating  Ga  into  the  absorber  for  optimized  absorption  of  sunlight. 


The  current  state-  of-the-art  device  was  fabricated  in  NERL  with  a total  area  efficiency  of 
18.8%  [Con  99].  Figure  1 . 1 shows  a cross  section  of  a common  thin  film  CulnSei-based 
photovoltaic  device  structure.  The  low  cost  soda-lime  glass  serves  as  a substrate.  A 
molydenum  back  contact,  usually  around  lpm  thick,  is  deposited  by  sputtering  or  c-beam 
evaporation.  The  Cu(ln,Ga)Se:  absorber  layer  is  normally  grown  by  co-evaporation  of 
the  constituent  elements.  Following  the  absorber  layer,  a thin  buffer  layer  (CdS  or  other 
Cd-free  alternatives)  is  prepared  by  chemical  bath  deposition.  The  n-type  window  layer 
normally  consists  of  a thin  layer  of  resistive  ZnO,  followed  by  a thicker  and  conductive 
layer  of  ZnO  and  finished  with  an  Anti-Reflection  Coating  (ARC)  (e.g.  MgFj).  For  small 
cells,  a metal  grid  is  formed  on  top  of  tile  ARC.  In  a typical  grid,  a thin  Ni  layer  is  first 
evaporated  to  form  an  Ohmic  contact,  followed  by  a thick  Al  layer  to  form  the  major 
body  of  the  contact- 


Figure  1 I Schematic  of  a typical 


CulnSerbased 


film  photovoltaic 


For  large-area  devices  (modules),  Ihe  individual  cells  are  series-connecied  by  monolithic 
integration  rather  than  use  of  metallic  top-contact  grids.  The  series  interconnection 
between  adjacent  cells  relies  on  the  front  contact  TCO  and  the  back  contact  Mo  layers. 
Figure  1.2  displays  a schematic  cross-section  of  a typical  module  structure  with  the 
equivalent  circuit.  A review  of  module  design  analysis  is  given  in  Stanbery  el  al  [Sta  85]. 
The  focus  of  this  thesis  is  on  the  CIS  absorber  layer. 


Figure  1 .2  Schematic  < 


i of  the  module  structure  with  the  equivalent  circuit. 


13  Current  Underslandinn  and  Future  R&D  Challenges 


1.3.1  SjiJislBteandjiolybdenunt 

Several  materials  have  been  tested  be  used  as  a substrate  for  CIS-based 
photovoltaic  device,  including  glass,  alumina,  selected  polymers,  and  stainless  steel.  The 
high-efficiency  cells  have  been  fabricated  using  low  cost  soda-lime  glass  (SLG), 
probably  a result  of  the  positive  impact  of  Na  out-diffiision  and  its  influence  on  grain 
growth  during  formation  of  the  absorber.  Experiments  [Ull  97J  have  shown  the  CIS- 
based  PV  performance  strongly  depends  on  the  properties  of  the  Mo  contact  The  Mo 
back  contact  is  normally  fabricated  by  dc  magnetron  sputtering.  Adhesion  and  resistivity 
are  the  two  properties  affecting  the  Mo  contact  quality  [Sco  95]  Another  important  role 
of  the  Mo  layer  is  the  interplay  between  the  Mo  and  underlying  soda-lime  glass  substrate. 
The  microstructure  of  the  Mo  film  controls  the  extent  of  Na  diffusion  into  the  absorber 
layer.  Current  research  topics  in  this  area  include  obtaining  a better  understanding  of  the 
relation  between  the  microstructure  of  Mo  film  and  Na  out-difiusion  from  SLG  substrate, 
optimizing  the  Mo  deposition  process  to  tailor  a desirable  microstructure  [Sco  95],  and 
the  search  for  a good  Na  diffusion  barrier  [Pro  96]. 

1.3.2  Absorber 

A wide  variety  of  techniques  have  been  used  to  fabricate  CulnSei-based  thin 
films.  There  are  two  main  types  of  deposition  processes:  co-evaporation  and  sequential. 


Co-evaporation  is  the  most  successful  technique  to  fabricate  high  efficient 
CulnSc2-based  solar  cells  since  the  demonstration  of  an  1 1%  efficiency  cell  by  the 


Boeing  group  [Mic  82],  The  basic  co-evaporaiion  process  is  a single  layer  process.  This 
layer  is  deposited  by  first  heating  the  elemental  sources  and  the  substrates  to  the  desired 
temperature,  then  opening  the  shutter  to  initiate  the  deposition.  A high  substrate 
temperature  (~S50”C)  is  required  to  form  device  quality  films.  This  process  produces  a 
homogeneous  but  small  grained  polycrystalline  film  across  the  thickness.  The  co- 
evaporation  technique  allows  for  changing  the  constituent  elemental  flux  ratio  during  the 
deposition.  The  bilayer  process,  first  developed  by  Mickelsen  and  Chen  [Mic  80],  utilizes 
this  feature  to  produce  CIS  films  with  large  grains  Deposition  of  CIS  begins  with  a 
slight  copper  excess  flux  to  form  a CulnSej  and  Cur-xSe  mixed-phase  film  at  low 
temperature,  and  followed  by  an  indium  excess  flux  to  forthcr  re-crystallize  the  first  layer 
at  a higher  temperature.  The  success  of  this  recipe  is  attributed  to  a liquid  assisted  re- 
growth process,  which  is  a result  of  the  melting  of  Cu!.,Se  in  the  presence  of  excess 
selenium  A further  development  lead  to  a three  stage  process  [Gab  94],  This  process 
begins  with  the  formation  of  an  InjSei  layer.  This  In-Se,  layer  is  then  subjected  to  a Cu 
and  In  flux  to  form  CuInSe-(s)  and  Cu.Sefl)  mixture  in  the  second  stage  of  this  process. 
Finally  In  and  Se  fluxes  are  used  to  consume  the  excess  Cu,Sc  phase  to  form  a large 
grain  size  device  quality  film.  A similar  process  was  also  developed  for  CIGS  thin  film 
formation. 

The  characteristic  of  the  sequential  process  can  be  described  as  low  temperature 
precursor  film  deposition  followed  by  high  temperature  annealing  There  are  many  kinds 
of  sequential  processes  Selenization  and  Stacked  Elemental  Layer  (SEL)  processes  are 
the  two  most  studied.  Selenization  is  a two-stage  process.  First,  sequential  layers  of 


copper  and  indium  are  deposited  on  a molydenum  coaled  substrate.  This  metal  precursor 
is  then  exposed  to  either  HiSe  or  Se„  gas  at  high  temperature  to  form  CulnSe?  films.  The 
SEL  process  differs  from  the  selenization  one  by  incorporating  a Se  layer  into  the 
precursor  structure  [Oum  89]  CIS-based  PV  devices  with  efficiencies  > 10%  had  been 
achieved  by  many  groups  using  these  two  processes.  The  metal  precursor  layers  for  these 
two  processes  can  be  deposited  at  high  rate  by  a variety  of  deposition  techniques,  such  as 
sputtering,  evaporation,  electro-deposition,  e-beam  evaporation.  This  is  an  advantage  for 
large-scale  production,  however,  these  processes  are  difficult  to  control.  The  search  for  a 
belter  precursor  structure  is  an  active  research  area  [Cha  97], 

1-3.3  The  Role  of  CulniSe; 

CIS-based  solar  cells  have  been  thought  to  be  p-n  structures  (p-CIS.  n- 
CdS).  However,  such  a model  does  not  predict  the  failure  of  high-Ga-content  devices  A 
novel  model  has  been  proposed  by  Schmid  el  al  [Sch  93]  that  the  growth  kinetics  of 
PVD-grown  CulnSei  result  in  the  in-silu  formation  of  a surface  layer  of  CulniSe*  and 
that  this  CuInSe2/CuInjSej  heterojunction  is  the  true  heterojunction  partner  in  CIS 
devices  This  surface  material  Is  reported  to  show  compositions  near  CuInjScs  and  has 
been  referred  to  as  an  "Ordered  Vacancy  Compound”  (OVC).  This  significant  finding 
resembles  earlier  observation  of  a buried  junction  elucidated  by  electron-beam  induced 
current  (EBIC)  [Mat  86]  measurements.  However,  the  crystal  structure  of  CuInjSes  and 
its  phase  relation  with  CulnSe;  is  still  not  clear  at  this  point. 


1 ,3.4  Junction  Formation 


Junction  formation  is  a critical  step  in  forming  a good  photovoltaic  device.  The 
best  results  today  for  CIS-based  PV  have  incorporated  deposition  of  a thin  CdS  layer 
(-500  A)  by  chemical  bath  deposition  (CBD)  [Kan  98],  Considerable  effort  has  been 
directed  towards  understanding  the  mechanism  responsible  for  this  dramatic 
improvement  of  cell  performance.  Elucidation  of  the  role  of  CBD  CdS  will  lead  to  better 
device  engineering  and  suggesting  alternative  processes.  Several  explanations  have  been 
proposed  to  account  for  the  improved  photovoltaic  performance  over  other  processes: 

1.  A pinhole-free  conformal  layer  of  evaporated  CdS  requires  thicker  films  to 

completely  cover  the  rough  CIGS  absorber  surface 

2.  The  liquid  electrolyte  chemical  modifies  the  absorber  surface. 

- Cd  diffuses  into  CIS  lattice  forming  Cdcu  (Cd  on  Cu  sites)  donor  and  induces 
a type  conversion  (p  to  n)  [Kan  98]. 

- Cd  reacts  with  CIS  surface  forming  CdSe,  or  CdIn,Se,.  which  produced  a 
graded  interface  structure  [Kes92). 

The  first  effect  has  been  commonly  accepted  as  beneficial  impacts.  The  chemical 
modification  explanation  appears  to  play  a very  crucial  role  in  addition  to  the  first  three 
mechanisms.  Intensive  studies  have  been  carried  out  in  this  area.  Several  proposed 
explanations  based  on  the  interfacial  reactions  between  the  CIS  surface  and  the  chemical 
bath.  For  example,  it  is  suggested  that  exchange  reactions  at  the  near-surface  region 
between  the  CIS  and  the  Cd  or  S ions  in  the  chemical  bath.  Another  proposed  mechanism 
is  washing  or  etching  effect  that  removes  soluble  components,  (eg,  sodium  complex,  a 
superficial  oxide)  The  oxide  removal  postulate  upon  dipping  in  the  chemical  bath  has 


been  carefully  studied  [Kcs  92].  In  a more  recent  set  of  studies  on  the  impact  of  the  buffer 
layer  process,  the  beneficial  effects  of  a Cd  partial  electrolyte  treatments  on  the  CIS 
films.  There  were  two  interactions  suggested  for  the  Cd  partial  electrolyte  treatment.  One 
model  [Kes  92]  suggested  that  Cd  ions  react  with  the  CIS  surface  and  form  a CdSe  or 
Cdln.Se,  compound,  thus  reducing  the  density  of  interface  defects  The  other  suggested 
interaction  is  that  Cd  ions  difliisc  into  the  CulnSe:  lattice  and  form  Cdcu  donors,  thus 
inducing  a type  conversion  [Kan  98],  The  Cd  doping  model  had  been  proposed  by  [Mat 
87]  for  an  evaporated-CdS/single  crystal  CulnSei  device.  The  true  mechanism  for  the 
observed  type  conversion  was  not  clearly  resolved  due  to  the  complexity  of  the  defect 
chemistry  Furthermore,  the  interplay  of  the  In-rich  compounds  (eg.  CuIniSes)  with  the 
chemical  bath  could  also  be  a significant  factor  for  the  reported  high  efficient  cell  with 
the  buried  junction  between  CuhtsSes  and  CulnSe;. 

Impurities  (e  g.  O,  N)  in  CBD  CdS  films  have  been  reported  [Kyi  96],  Several 
compounds  such  as,  Cd(OH);,  CdO,  CdNCN,  CdSOr,  CdCOj.  Cl,  and  H;0  have  been 
suggested  The  disagreement  for  the  impurity  compounds  depend  on  the  recipes  used  for 
the  bath  and  die  difficulty  in  interpretation  ol  the  analvsis  results  The  effects  of 
impurities  on  the  cell  performance,  however,  are  not  clear  so  far. 

1.3.5  Alternative  buffer  layers 

Considerable  effort  has  been  devoted  to  find  a replacement  for  the  CBD  CdS  process 
in  CIGS-based  solar  cells.  This  search  has  been  directed  towards  developing  a Cd-free 
process  for  environmental  reasons,  and  a dry  deposition  process  to  make  it  more 
compatible  with  the  remaining  of  manufacturing  process.  A variety  of  buffer  layer 
materials  have  been  tested  as  replacements  for  CdS,  including  ZnO,  ZnS,  Zn(0,S,0H)„ 


SnOj,  Sn(S,0)2,  In:Si,  In(OH)j,  Iri,(OH.S),.  ZnSe , and  ZrO;  all  deposited  by  CBD.  Other 
approaches  include  ZnO  by  atomic  layer  epitaxy  and  ZnSe  by  MOCVD  CIGS  solar  cells 
with  continuously  evaporated  Cd-ffee  buffer  layers  ln,Se,.  Gasses,  and  ZnIn,Sey  have 
also  been  demonstrated.  Although,  Significant  progress  has  been  achieved  in  this  area, 
these  alternative  buffer  layer  materials  and  processes  still  produce  inferior  CIS-based  PV 
cells,  and  some  suffer  from  stability  problems  (e.g.,  Zn(0,S,0H), , Znln.Se, ) The  search 
for  a better  alternative  buffer-layer  process  is  still  needed. 

1.3.6  TCP  layer 

The  role  of  the  TCO  layer  is  to  effectively  transport  photo-carriers  generated  by 
the  absorber  layer  to  the  external  circuit  without  the  sacrifice  of  incident  sunlight.  ZnO  is 
the  common  choice  for  the  top  contact.  Criteria  for  the  optimization  of  the  ZnO  film  are 
low  sheet  resistance  and  high  optical  transmission.  The  figure  of  merit  can  be  defined  as 
the  ratio  of  the  electrical  conductivity  to  the  optical  absorption  coefficient  of  the  film.  RF 
sputtering  or  MOCVD  are  the  two  widely  used  deposition  processes.  The  ongoing 
research  for  an  improved  TCO  layer  includes  elucidation  of  the  relation  between 
microstructure  (point  and  planar  defects)  and  processing  parameters,  novel  TCO 
materials  and  processes,  and  a thinner  TCO  layer. 

I 4 Statement  of  Thesis  Work 
Tile  objectives  of  this  thesis  are: 

I.  Develop  a detailed  thermodynamic  description  of  the  Cu-ln-Se  system  for  better 
understand  the  processing  chemistry  of  CulnSei  absorber  layer  formation. 


Apply  the  obtained  knowledge  of  thermochemistry  and  phase  diagram  to  design 
novel  precursor  structure  for  synthesis  of  CulnSe2  absorber  layer  by  rapid  thermal 

Investigate  the  local  structure  of  important  copper  indium  selenidc  semiconductor 
alloys  (CuInSej,  CulnjSe»,  Culm  Sc*)  by  EXAFS. 


CHAPTER  2 

THEORY  AND  PROCEDURES  FOR  ASSESSING 
THERMOCHEMISTRY  AND  PHASE  DIAGRAM  DATA 


2.1  Imroduction 

Van  Laar  [Van  08]  established  the  basis  of  calculating  phase  diagrams  in  1908  by 
using  a solution  model  to  describe  the  Gibbs  energy  of  solution  phases  and  incorporate 
I them  into  the  calculation,  -ith  the  introduction  of  large  scale  computational  capabilities, 
these  approaches  can  now  be  applied  to  large  scale  systems.  The  CALPHAD  method  is 
characterized  by  the  combination  of  Computer  CALculations  of  PHAse  Diagram  with 
experimental  thermochemical  and  phase  diagram  data.  The  main  output  of  this  procedure 
is  to  generate  useful  phase  equilibria  and  thermochemical  information  to  solve  problems 
in  materials  processing.  In  1973,  a group  (CALPHAD)  was  organized  to  share 
experiences  and  further  developed  this  approach  [Sau  98].  Initially,  systematic  trial  and 
error  methodologies  were  used  to  adjust  model  parameters.  As  the  computational  power 
increased,  it  became  possible  to  treat  large  data  set  using  formal  parameter  estimation 
algorithm  The  first  of  these  was  created  by  Dr.  H.L.  Lukas  [BINGSS  and  BINFKT]  and 
applicable  to  binary  systems  [Luk  77]  These  programs  allowed  simultaneous  adjustment 
of  the  model  parameters  to  all  types  of  measured  data  such  as  phase  diagrams  and 
thermochemical  data,  and  to  account  for  errors  of  measurement.  Several  sophisticated 
programs  have  been  developed  since  then,  including  Lukas  [Luk  92],  Thermo  Calc  [JOn 


90],  MTDATA  [Dav  90]  and  FACT  [Pel  90], 


2.2  Principle  of  Phase  Piagram  Calculalic 


Most  computer-based  calculations  of  phase  diagrams  are  based  on  the  principle  of 
Gibbs  energy  minimization  and  Maxwell’s  construction  (or  common  tangent 
construction).  The  Maxwell  construction  is  a geometric  interpretation  of  the  solution  for 
equilibrium  conditions.  Which  are: 

T.=T. 

(2.0 

P°=P>  (2.2) 

(23) 

/>!=?;  (24) 

for  a binary  (components  A and  B)  two-phase  (phases  a and  (i)  system.  An  illustration  of 
this  construction  is  given  in  Figure  2. 1 . 


Figure  2. 1 An  illustration  of  Maxwell's  construction. 


2 3 CALPHAD  Annrnaches 


The  CALPHAD  approaches  is  based  on  three  elements:  data,  models,  and 
software,  as  shown  in  Figure  2.2. 


Figure  2.2  Elements  and  functions  of  the  Calphad  approach. 


a (e.g,  AHf,  AH„ix,  AH,™,,  Cp) 


2.  Enthalpy  data 


2.3.2  The  models 

A thermodynamic  description  of  a system  requires  the  assignment  of 
thermodynamic  functions  for  each  phase.  The  CALPHAD  method  employs  various 
models  to  describe  the  Gibbs  energy  as  a function  of  temperature,  pressure,  and 
composition. 

2.3.2. 1 Reference  stale 

Energy  quantities  (eg,  enthalpy.  Gibbs  energy)  can  not  be  defined  absolutely. 
Only  their  differences  have  a physical  meaning.  By  convention  the  enthalpy  for  end 
element  in  its  stable  state  at  298.15  K and  I bar  H»>  is  assigned  the  value  zero.  The  use 
of  this  reference  for  data  is  called  Stable  Element  Reference  and  denoted  SER  and  the 
enthalpy  is  sometimes  denoted  Hser. 

2.3.2  2 Models  .for  phases  wjlh  fixed  composition 

A pure  element  and  a line  compound  cannot  vary  in  composition  and  thus  its 
Gibbs  energy  is  as  a function  of  temperature  and  pressure  only.  For  such  phases  the 
Gibbs  energy  is  expressed  as: 

G-Hsbi  = a*hT*  cTlnT + «/7*  + e/T+ fl‘  (2.5) 

The  Gibbs  energy  for  a solution  phase  can  be  divided  into  three  contributions  as 
the  following  equation: 

G = G°  + G'J  + Gu  (2.6) 

Where  G°  corresponds  to  the  Gibbs  energy  of  a mechanical  mixture  of  the  constituents  of 
the  phase;  0““.  is  the  entropy  of  mixing  for  an  ideal  solution,  and  the  term;  and  G"  is 
the  excess  Gibbs  energy,  which  accounts  for  the  departure  from  ideal  mixture  behavior. 


Many  models  have  been  developed  to  describe  this  non  ideality,  including  the 
regular  solution  model,  associated  model,  ionic  model,  -agner  Schottky  model, 
interstitial  model,  and  sublattice  models. 

One  commonly  used  model  is  the  Redlich  Kister  formalism.  The  Gibbs  energy, 

G,  fora  binary  solution  is  expressed  as: 

G - «“*  =CGa -Hf*)xA  +(‘GB-HfR)xB  + «'/'(xj  In  xA  + x„  In  .<■„)+ xAxB  1K,(xa  - xB)" 
(2.7) 

where  xA  and  xn  are  mole  fractions  of  elements  A and  B,  and  K,  is  a function  of 
temperature,  according  to: 

K,  = a,  + b»T  + c,T  InT  (2.8) 

Detailed  discussions  of  these  various  models  can  be  found  in  (Sun  90).  Recently 
efforts  have  been  directed  to  introduce  ab  inlio  calculations  of  phase  diagrams  into  the 
CALPHAD  community  (And  96J,  like  the  Cluster  Variation  Method  (CVM),  or  Monte 
Carlo  (MC)  calculation  [Fon  96].  These  methods  are  based  on  knowledge  of  the  atomic 
configuration  and  accompanying  electronic  structure  calculations,  promise  to  bring  more 
insight  into  appropriate  solution  models  for  systems  that  are  difficult  to  probe 
experimentally. 

2.3,3  Software 

Two  major  algorithms  are  involved  in  the  calculation,  one  is  the  optimizer  based 
on  the  Gaussian  least  squares  method,  and  the  other  is  Newton  Raphson  method.  The 
software  calculates  the  phase  diagrams  as  well  as  thermodynamic  functions  by  using  the 
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optimized  Gibbs  energy  function  from  the  optimizer  The  input  to  the  assessment  is 
thermodynamic  and  phase  equilibrium  data,  the  output  is  optimized  expressions  for  the 
Gibbs  energy  of  each  phase.  The  calculation  flowchart  is  given  in  Figure  2.3. 


Figure  2.3  Calphad  thermodynamic  assessment  flowchart. 


CHAPTER  3 

THERMODYNAMIC  ASSESSMENT  OF  SELENIUM  UNARY  SYSTEM 

3.1  Introduction 

For  the  thermodynamic  assessment  of  phase  diagrams,  the  phase  stability  of  the 
pure  components  are  required.  Furthermore  knowledge  of  the  vapor  phase  composition  is 
important  in  process  planning  Therefore  reliable  G-Hsim  expressions  for  the  condensed 
and  gaseous  selenium  species  are  necessary.  Critically  assessed  G-Hskr  expressions  for 
the  selenium  species,  however,  are  lacking  in  the  literature.  For  this  reason,  the 
thermodynamic  properties  of  Sc  (c.  hexagonal).  Se  (I)  and  Se„  (g)  (n=  I -S).  including 
AHuu9s.i5k,  S°298  i5K.  and  Cp(T)  were  assessed.  Using  these  properties  the  G-Hser 
expressions  for  each  phase  of  selenium  were  produced.  The  temperature  dependence  of 
vapor  pressures  and  compositions  for  Se„(g)  in  equilibrium  with  condensed  selenium 
were  also  calculated  from  equilibrium  constants  for  Se„(g)  and  compared  with 
experimental  data. 


3.2  Literature  review 

3.2.1  Solid  state 

Three  modifications  of  solid  selenium  have  been  reported;  hexagonal,  monoclinic, 
and  amorphous  selenium.  Hexagonal  selenium  is  taken  as  the  stable  standard  state.  The 
heat  capacity  of  hexagonal  selenium  had  been  reported  by  many  authors,  which  were 


reviewed  by  Gaur  el  ai  [Gau  81].  The  earlier  assessments  by  Hulgren  [Hul  73],  Mills 
[Mil  74]  and  SGTE  [Din  91]  did  not  include  the  heat  capacity  data  from  [Gro  73]  and 
[Cha  74],  No  experimental  data  are  available  for  the  solid  selenium  superheated  above  its 
melting  point  therefore  only  an  estimation  can  be  made  for  the  superheated  heat  capacity 

The  value  of  S°2<m  was  evaluated  by  both  [Gro  84]  and  [Gau  81],  with  the  value 
by  Gaur  era/.  [Gau  81]  lower.  This  discrepancy  was  produced  by  the  fact  that  the  heat 
capacity  data  measured  by  Chang  and  Bestul  [Cha  74]  in  the  range  of  4 to  30  K was  used 
by  Gronvold  e/  ai  [Gro  84]  but  not  by  Gaur  et  al  [Gaur  81],  who  used  Breston  and 
Lasjaunias's  [Bre  73]  data  instead  The  assessed  value  by  Gronvold  cl  ai  [Gro  84]  is 
adopted  in  this  work. 

3.2.2  Liquid  state 

Several  investigators  have  measured  the  heat  capacity  of  liquid  selenium.  The 
data  from  Mondain-Monval  [Mon  26],  Brizgys  [Bri  47],  and  Gattow  and  Heinrich  [Gat 
74]  are  in  poor  agreement.  The  later  measurements  by  [Gro  73],  [Cha  74],  [Meh  78],  and 
[Shu  80]  are  in  better  agreement.  They  investigated  the  heat  capacity  from  the 
supercooled  region  (300  to  494  K)  to  1000  K The  assessed  liquid  selenium  heal 
capacity,  CP(T),  along  with  the  experimental  data  are  shown  in  Figure  3.2.  The  heat 
capacity  decreases  in  the  temperature  range  300  to  760  K,  with  a subsequent  increase  to 
1000  K.  This  complex  behavior  is  due  to  vibrational  energy  contributions,  free  volume 
effects  and  heats  of  reactions  from  the  change  in  ring-chain  equilibria  and 
depolymerization  equilibria,  which  have  been  studied  by  Shu  cl  ai  [Shu  80].  No 
experimental  data  were  available  for  temperature  above  1000  K.  Constant  heat  capacity 


values  were  used  above  this  temperature  in  both  SGTE  [Din  91]  and  IVTAN  [Gur  78] 
databases. 


3.2.3  Gaseous  stale 

Similar  to  sulphur,  selenium  exhibits  a complex  composition  in  the  vapor  phase. 
The  mass-spectrometric  measurements  have  shown  the  presence  of  all  the  molecules 
from  monomer  to  Seio  species,  the  more  abundant  ones  being  Se?  through  Seg.  Studies  of 

These  include  vapor  density  [111  57],  and  [Rau  74],  torsion  or  Knudsen  eflusion 
measurements  [Yam  57],  [Rat  62],  [Kel  71],  optical  absorption  [Bre  65],  and  mass- 
spectrometry  [Ber  66],  [Ber  68],  [Yam  57],  [Hor  75],  [Fuj  66],  [Rat  62],  [Kno  68],  [Kel 
71],  [Gri  82],  [Hua  84],  pro  84],  [Bre  65],  [Sau  71],  and  [Gni  74],  MUls  [Mil  74] 
reviewed  the  data  available  before  1974.  Then  Gronvold  el  aL  [Gro  84]  assessed  the 
AHcot.  molecular  constants  and  thermodynamic  functions  for  selenium  monomers  and 
dimers.  Steudel  [Ste  8 1]  calculated  the  thermodynamic  functions  for  Ses  molecules  from 
spectroscopic  and  structural  data.  The  entropy  and  heat  capacity  data  for  cyclic  Sea 
(n-5,6,7)  molecules  were  obtained  from  linear  relationships  between  ring  size  n,  entropy, 
and  heat  capacity  values  of  Sez  and  Ses.  Rus  [Rus  91]  calculated  the  thermodynamic 
functions  for  Sea  (q33  to  7)  from  molecular  constants  derived  from  the  data  of  Ser,  Seg 


IS„(n=2to8). 


The  CP(T)  expression  for  solid  hexagonal  selenium  was  modeled  as  Cp  = a + 
bT  + cT!.  and  the  parameters  were  estimated  using  the  experimental  data  of  [Gro  73),  and 
[Cha  74).  The  coefficients  are  given  in  Table  3.10.  For  the  superheated  solid  selenium  the 
heal  capacity  was  simply  extrapolated  to  a constant  value  based  on  the  polynomial 
expression  at  760  K where  the  Cr(T)  curve  reaches  its  maximum  value  at  28.55  J/K-mol. 
The  coefficients  for  G-Hser  derived  from  CP(T)  expressions  using  the  entropy  value 
Scmi>  = 42.27  J/K-mol  by  Gronvold  el  al.  (Gro  84]  are  given  in  Table  3.12.  The  plot  of 
Cp  vs.  T along  with  experimental  data  is  shown  in  Figure  3.1.  The  reported  values  for  the 
heat  of  fusion  are  widely  divergent  [Mil  74).  The  experimental  value  AHllvm  - 6I59±4 
J/mol  reported  by  Gronvold  (Gro  73]  was  used  here. 

3.3.2  Liquid  stale 

the  liquid  state  heat  capacity.  CP(T),  was  expressed  as  one  equation  ( Cp  = a + bT 
+ cT! ) using  the  experimental  data  measured  by  (Gro  73],  (Cha  74],  and  [Shu  80]  from 
the  super-cooled  state  298  to  494  K to  1000  K.  No  experimental  data  are  available  for 
temperature  higher  than  1000  K.  The  extrapolation  was  carried  out  by  smoothing  the 
CP(T)  based  the  polynomial  expression  curve  to  a constant  value,  36  J/K-mol,  at  1200  K. 
The  fitted  heal  capacity  expression  is  shown  in  Table  3.10.  Cp  and  G-Hser  vs.  T plots 
with  experimental  data  and  SGTE  (UNARY)  database  are  given  in  Figure  3.2  and  3.3 
respectively. 


v [Gro  731 
|And  37| 


Figure  3. 1 Fitted  heat  capacity  data  for  hexagonal  selenium  along  with  experimental  data 
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3.3.3  Gaseous  siaie 

Gronvold  el  al  [Gro  84]  have  critically  assessed  the  standard  enthalpy  of 
formation  and  thermodynamic  functions  for  the  monomer  and  dimer.  Their  data  are 
adopted  here  For  Sc„  (n=3  to  8),  Drowart  and  Smoes  [Dro  84]  and  Keller  [Kel  70] 
calculated  the  value  of  S°2w.  The  molecular  constants  calculated  by  Keller  [Kel  70]  arc 
based  on  rough  approximations  No  information  was  given  in  [Dro  84]  concerning  the 
molecular  structures  and  vibratonal  frquencies  used  in  their  calculation.  The  heat  capacity 
data  listed  in  Mills  [Mil  74]  and  Rau  [Rau  74]  are  based  on  the  data  in  [Dro  84]  and  [Kel 
70]  respectively.  Rusin  [Rus  91]  estimated  the  vibrational  frequencies,  product  of  inertia 
of  the  Se„  molecules  (n=3  to  7)  from  the  molecular  constants  of  S„  and  Se2,  Ses  [Gttr  78], 
[Gro  84],  and  [Stc  81]  based  on  the  Molecular-Orbital  as  a Linear  Combination  of 
Atomic  Orbitals  (MO-LCAO)  method.  Nagata  et  al  [Nag  81]  assigned  the  complete 
vibrational  spectra  for  Se,:  from  Raman  and  IR  measurements.  The  estimated  frequencies 
by  Rusin  [Rus  91]  and  the  experimental  values  determined  by  Nagata  el  al  [Nag  81]  for 
Set  are  in  good  agreement.  The  molecular  constants  listed  by  Rusin  [Rus  91]  together 
with  the  data  by  Nagata  el  a I [Nag  81)  are  adopted  to  calculate  the  thermodynamic 
functions  for  Se,,  (n=3  to  8),  assuming  harmonic  oscillator,  rigid  rotator  and  ideal  gas 
behavior. 

The  values  of  standard  enthalpy  of  formation  for  Sen  (n=3  to  8)  were  investigated 
by  several  authors  employing  various  techniques.  The  partial  pressures  resulting  from  the 
electrochemical  Knudsen  cell  by  Keller  el  al  [Kel  71]  are  more  accurate  than  those  from 
conventional  mass-spectrometric  investigations,  especially  for  small  molecules  (So.  So 
and  So  ) which  can  be  formed  by  fragmentation  of  large  molecules  even  at  the  lowest 
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energies.  The  electrochemical  Knudsen  cells  allow  parent  ions  and  fragments  to  be  better 
distinguished.  When  the  logarithm  of  the  ion  intensity  in  the  mass-spectrometry  is  plotted 
against  the  potential,  the  different  species  should  have  different  slopes.  Only  they  when 
are  in  agreement  with  the  theoretical  value,  are  then  thought  to  be  the  parent  ion  intensity 
The  vapor  density  measurements  by  [III  57]  only  included  the  even  numbered  species  in 
the  vapor,  thus  the  values  are  disregarded.  Rau  [Rau  74]  measured  the  vapor  density  and 
total  vapor  pressure  over  liquid  selenium  and  optimized  the  thermodynamic  properties, 
the  enthalpies  and  entropies  of  Se„  (n=2  to  8)  by  fitting  the  vapor  pressure  and  density 
data  Therefore  the  set  of  thermodynamic  properties  for  Se„  listed  by  [Rau  74]  do  not  lead 
to  independent  values. 

To  evaluate  the  standard  enthalpies  of  formation  for  Sc„  (n=3  to  8),  the  calculated 
thermodynamic  functions  are  used  to  perform  second  law  and  third  law  analysis  of  the 
experimental  data  The  average  total  pressure  by  Neumann  and  Lichtenberg  [Neu  39]  and 
Niwa  and  Shibata  [Niw  40]  over  solid  selenium  is  used  to  calculate  the  partial  pressure  in 
third  law  analysis  The  assessment  for  each  species  is  stated  below. 
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Figure  3.2  Assessed  heat  capacity  for  liquid  selenium  along  with  experimental  data 
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3.3.3. 1 Se 


The  evaluated  values  of  standard  enthalpy  of  formation  for  Sei  (g)  from  [Stu  56], 
[Mil  74].  and  [Gro  84]  are  listed  in  Table  3.1.  Gronvold  ct  al  [Gro  84]  had  critically 
evaluated  the  standard  enthalpy  of  formation  and  dissociation  energy  data  of  Se2  (g)  gas 
They  had  used  these  data  to  calculate  the  standard  enthalpy  of  formation  for  Sei  (g)  gas. 
Their  assessed  value,  AHi;.,*  = 237.6I±0.67,  is  adopted. 


Table  3.1.  Standard  enthalpy  of  formation  data  for  Se,,) 


Authors 

Stull  and  Sinke  fStu  S6| 

Mills  [Mil  741 
Gronvold  el  al  [Gro  84] 


AtrYwann  (kJ/mol) 

206.69 


237.61  tO.67 


The  entropy  data  listed  in  [Stu  56],  [Hul  73],  [Mil  74],  and  [Gro  84]  are 
consistent.  The  value,  S°29g=  176.62  J/K-mol.  was  adopted  in  this  work.  The  heat  capacity 
data  from  different  authors  are  shown  in  Figure  3.5.  The  heat  capacity  data  from  [Gro  84] 
have  been  fitted  with  CP(T)  expressions  with  one  break  point  according  to  the  shape  of 
the  CP(T)  curve.  Results  of  CP(T)  and  G-HssaCO  are  presented  in  Tables  3.10,  3.13  and 
Figure  3.6  respectively. 

3.3.32  Ssa 

The  value  of  the  heat  of  formation  of  Se2(g)  has  been  assessed  by  many  authors 
[Stu  56],  [Hul  73],  [Mil  74],  [Gro  84],  [Ven  82],  The  critical  assessed  value,  AHr.29*  = 
144.14  II  K*mol,  by  Gronvold  el  al  [Gro  84]  is  adopted  here.  The  entropy  and  heat 
capacity  data  calculated  from  the  molecular  constants  were  performed  by  [Stu  56),  (Hul 


73],  [Mil  74),  [Gro  84],  [Rau  74],  and  [Ven  82],  The  reported  heat  capacity 


shown  in  Figure  3.5.  The  latest  assessed  data  from  [Gro  84)  are  adopted  here  Two 
equations  are  used  to  fit  the  complex  shape  of  the  heat  capacity  curve.  The  equations  and 
plots  of  Cp(T)  and  G-Hsm(T)  are  given  in  Table  3.10, 3. 13  and  Figure  3.6. 


T(KJ 


|Gro  84 1 

|Hul  73| 

| Mil  74| 


Figure  3 .4  reported  heat  capacity  for  Se(g)  in  the  range  T=298  to  2000  K. 


30 


T [KJ 

[84Gro] 

[73Hul] 

[82Ven] 

— [74MM] 

[73Rau] 


Figure  3.5  He 


city  for  Sej(g)  in  the  range  T=298  to  2000  K. 
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Figure  3.6  Fitted  heat  capacity  for  Se(g)  and  Sei(g)  in  the  range  T=298  to  2000  K 


3,3.33  S® 


The  values  of  ihc  standard  enthalpy  of  formation  obtained  from  various 
investigations  for  Sej  are  listed  in  Table  3.2.  The  data  are  calculated  from  assessed 
thermodynamic  functions  by  second  and  third  law  methods. 


Table  3 .2.  Standard  enthalpy  of  formation  data  for  Se)l8,. 


Berkowit2  and  Chupka  [Ber  66] 


Keller  et  ai  [Kel7l] 


Huang  el  ai  [Hua  84] 


Mass  spectrometry 


II:  Second  Law  Value,  III:  Third  Law  Value 


The  values  derived  from  conventional  mass-spectrometric  studies  [Ber  66],  [Hua 
84]  are  lower  than  those  by  Keller  et  al.  [Kel  71],  The  reason  for  such  discrepancy  can  be 
ascribed  to  fragmentation  of  larger  molecules.  The  electrochemical  Knudsen  cell  used  by 
Keller  el  a I [Kel  71]  is  more  reliable  for  the  study  of  smaller  Se  gas  molecules  than 
conventional  mass  spectrometer  The  third  law  values  by  [Kel  71]  are  self-consistent  with 
a standard  deviation  of  0.25  kJ/mol  and  also  agree  with  the  second  law  value.  Thus,  the 


third  law  value,  AHnra  = 181 .44  ± 0.26  kJ/mol,  by  [Kel  71]  was  selected.  The  entropy 
value.  S°2M  = 313.90  J/K-mol,  calculation  by  [Rus9l]  is  adopted  in  this  work.  Figure  3,7 
shows  the  assessed  C,,(T)  along  with  data  from  other  authors  The  results  of  assessed 
CP(T)  and  G-Hser(T)  are  given  in  Tables  3. 1 0 and  3. 13,  respectively. 
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7 Assessed  heal  capacity  for  Sej(g)  and  Se,(g)  in  the  range  T-298  lo  2000  K. 


Figure  3.' 


3.3.3  4 Se* 


The  values  of  the  standard  enthalpy  of  formation  obtained  from  various 
investigations  are  listed  in  Table  3.3.  The  entropy  value,  S°2m  = 343.32  J/K-mol,  reported 
by  [Rus  91]  is  used  in  this  work. 


Table  3.3  Standard  enthalpy  of  formation  data  for  Se*  (g). 


Keller  era/.  [Kol  7 1 J 


Mass  spectrometry  and 


Vapor  density  measurements 


Grimley  el  al.  [Gri  82] 


Huang  el  al  [Hua  84] 


II:  Second  law  value.  III:  Third  law  value 


As  with  Sc)  the  standard  enthalpy  of  formation  values  reported  by  [Gri  82]  and 
[Hua  84]  are  lower  than  the  values  of  Keller  er  al.  [Kel  71]  due  to  the  fragmentation  of 
larger  molecules  Although  the  experimental  methods  used  by  [Kel  71]  arc  more  reliable 
in  general,  the  third  law  analysis  shows  a standard  deviation  3.59  kj/mol  And  the  second 
law  value  seems  to  be  too  large.  The  standard  enthalpy  of  formation  value  for  Sej  needs 
frirthcr  refinement.  Se*.  however,  is  not  a major  component  in  the  vapor  phase  in 
equilibrium  with  condensed  selenium.  It  is  not  sensitive  to  either  the  total  vapor  pressure 
or  density.  An  optimization  through  total  vapor  pressure  or  density  data  is  not  reliable. 
Thus,  the  third  law  value,  AHu»  = I67.27±3.77  kJ/mol,  reported  by  [Kel  71]  is  selected 
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here.  Figure  3.7  shows  the  assessed  CP(T)  along  with  data  from  other  authors  The  results 
of  the  assessed  Cr(T)  and  G-Hsek(T)  are  given  in  Tables  3. 10  and  3.13,  respectively. 
3.3.35 

The  values  of  the  standard  enthalpy  of  formation  obtained  from  various 
investigations  are  listed  in  Table  3.4.  The  entropy  value,  S°i9s  = 4 1 5.50  J/K-mol,  reported 
by  [Rus  91]  is  adopted  in  this  work. 

The  third  law  values  derived  from  [Ber  68],  [Fuj  66],  [Gri  82],  [Hua  84],  (Dro  84] 
are  in  good  agreement.  The  third  law  value  derived  from  [Kel  71]  with  a standard 
deviation  of  2.71  kj/mol  seems  to  be  too  high.  From  the  mass-spectrometric  studies  it  is 
known  that  Se>  is  a major  species  in  equilibrium  with  liquid  selenium  in  the  temperature 
range  500-1000  K.  The  standard  enthalpy  of  formation  value  for  Sej  can  be  assessed  by 
the  total  vapor  pressure  data  in  equilibrium  with  liquid  selenium  in  the  temperature  range 
500  to  1000  K.  The  vapor  pressure  data  by  [III  57],  [Rau  74],  and  [Bak  68]  are  in  good 
agreement  within  this  temperature  range.  Their  data  are  used  to  perform  the  optimization. 
The  calculation  is  performed  by  using  the  assessed  G-Hsek  equations  for  liquid  selenium 
and  gaseous  Se„(n=l-8.except  5)  to  calculate  the  partial  pressures  for  each  species  And 
the  partial  pressure  for  Sej  is  obtained  by  subtracting  the  partial  pressures  of  Se*  (n=  I to 
8,  except  5)  from  the  total  vapor  pressure.  Since  the  enthalpy  of  formation  value  for  Sej  is 
uncertain,  the  value  was  changed  from  163  to  171  kJ/mol  in  calculating  the  partial 
pressure  for  Sej.  The  Pcng-Robinson  equation  of  state  is  used  for  real  gas  corrections  and 
the  ideal  mixture  behavior  is  assumed.  The  calculated  results  show  the  AHt»s  for  Sej 
equal  153.63  to  153.70  kJ/mol  when  AHrjw  for  Sej  change  from  171  to  163  kj/mol.  The 
mean  value  of  AHf.2»,=153.66±0.80  kJ/mol  is  selected 


Tabic  3.4.  Standard  enthalpy  of  formation  data  for  Sej  (g). 


Berkowitz  and  Chupka  [Ber  66] 


Fujisaki  el  al.  [Fuj  66] 


Yamadagni  and  Porter  [Yam  57] 


Keller  elal.  [Kel  71] 


Hoareau  el  al.  [Hor  75] 


Grimley  el  al.  [Gri  82] 


Huang  el  al.  [Hua  84] 


Drowait  elal.  [Dro  84] 


Mass  spectrometry 


II:  Second  law  value,  III:  Third  law  value 


Figure  3.8  showed  the  assessed  CP(T)  along  with  data  from  other  authors.  The  results  of 
evaluated  Cp(T)  and  G-HSeb(T)  are  given  in  Tables  3.10  and  3. 13,  respectively 


3.3.36  S® 

The  values  of  the  standard  enthalpy  of  formation  obtained  from  various 
investigations  are  listed  in  Table  3.5.  The  entropy  value,  S°»i  = 438.71  J/K-mol,  was 
calculated  from  the  molecular  constants  listed  by  [Rus91]  and  [Nag  81]. 

The  third  law  values  of  AHqm  derived  from  the  data  of  [Ber  66],  [Fuj  66],  [Rat 
62],  [Gri  82],  [Hua  84],  and  [Dro  84]  arc  in  relatively  good  agreement.  The  third  law  and 
second  values  by  [Kel  71]  are  self-consistent  with  a standard  deviation  value  of  0.25 
kJ/mol,  however,  smaller  than  others'  are.  Se„  (g)  is  the  most  prominent  species  in 
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equilibrium  with  solid  selenium.  It  was  found  the  calculation  Sc,.  <g)  vapor  pressure  from 
the  value  by  [Kel  71]  is  higher  than  the  total  pressure  in  equilibrium  with  solid  selenium 
as  measured  by  (Yam  57],  [Neu  39],  [Niw  40],  and  [Dod  20]  using  a torsion-effusion 
measurement  On  the  other  hand  the  Ses  (g)  vapor  pressure  calculated  from  the  mean 
AHr.;y*  value  from  [Ber  66],  [Hor  75],  [Fuj  66],  [Gri  82],  [Hua  84],  and  [Dro  84]  is 
consistent  with  those  experimental  measurements  Thus,  the  mean  value,  = 14 1 .26 
± 0.47  kJ/mol,  is  selected  here.  Figure  3.8  showed  the  evaluated  Cp(T)  along  with  data 
from  other  authors  The  results  of  C,(T)  and  G-Hsrjt  (T)  are  given  in  Tables  3.10  and 
3.13,  respectively. 


Table  3.5  Standard  enthalpy  of  formation  data  for  Se<,  (g). 


Authors 

Method 

(kJ/mol) 

Berkowilz  and  Chupka  [Ber  66] 

Mass  spectrometry 

145.55  (II) 
141.67  (111) 

Fujisaki  et  al  [Fuj  66] 

Mass  spectrometry 

144  66  (11) 

Yamadagni  and  Porter  [Yam  57] 

Mass  spectrometry 

149  47  (II) 

Ratchford  and  Ricken  [Rat  62] 

Mass  spectrometry  and 
electrochemical  Knudsen  cell 

137.32  ( 11) 

Keller  era/ [Kel  71] 

Mass  spectrometry  and 

electrochemical  Knudsen  cell 

137.60  (II) 
138.80  (III) 

Rau  [Rau  741 

vapor  density  measurements 

Moreau  et  al.  [Hor  75] 

Mass  spectrometry 

Grimley  el  al.  (Gri  82] 

Angular  distribution  mass 

spectrometry 

135.42  (II) 

Huang  cl  al.  (Hua  84] 

Mass  spectrometry 

143.13  (11) 
141.20  (III) 

Drowatt  el  al.  [Dro  84] 

Mass  spectrometry 

137.38  (11) 
140.94  ( HI ) 

11:  Second  law  value,  III:  Third  law  value 
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Figure  3.8  Heal  Capacity  of  Sej  (g)  and  Ses  (g)  in  the  range  T=298  to  2000  K. 


3.3  3.7  Ser 


The  values  of  the  standard  enthalpy  of  formation  obtained  from  various 


mvesugations  are  listed  in  Table  3.6.  The  entropy  value.  S°!M  = 491.35  J/K-mol,  from 
(Rus  91]  is  selected  here.  The  third  law  values  derived  from  [Kel  71]  are  self-consistent 


with  a standard  deviation  value  0.38  kj/mol  and  agree  well  with  the  second  law  value. 
The  value  AH|;>x  =149.84  ± 0.40  kJ/mol  is  solccted.  Figure  3.9  shows  the  evaluated 
Cp(T)  data  along  with  data  from  other  authors  The  results  of  Cp(T)  and  G-Hser  <T)  are 
given  in  Tables  3. 10  and  3.13,  respectively. 


Table  3.6  The  standard  enthalpy  of  formation  data  for  Se^g. 


Authors 

Method 

AH"r.Wi,K 

(kJ/mol) 

Berkowitz  and  Chupka  [Ber  66] 

Mass  spectrometry 

159  48  (II) 

152.48  (III) 

Fujisaki  el  al.  [Fuj  66] 

Mass  spectrometry 

164.15  (II) 

152.11  ( III  ) 

> amadagni  and  Porter  [\am  57] 

Mass  spectrometry 

Keller  el  al.  [Kel  71] 

Mass  spectrometry  and 

148.45  (11  ) 

Hoareau  el  at.  [Hor  75] 

Mass  spectrometry 

152  80  ( n ) 

Grimley  el  al.  [Gri  82] 

Angufardi5tribution  mass 

spectrometry 

149.23  (II) 

155.18  (UI) 

Huang  el  al.  [Hua  84] 

Mass  spectrometry 

159  46  ( U ) 

152.19  ( 111 ) 

Drowatt  el  al.  [Dro  84] 

Mass  spectrometry 

150.81  (II) 

152.57  (111) 

II:  Second  law  value.  III:  Third  law  value 


3.3.38  £e> 

The  values  of  the  standard  enthalpy  of  formation  obtained  from  various 
investigations  are  listed  in  Table  3.7.  The  entropy  value.  S°»*  = 532.24  J/K-mol,  from 
[Rus  91]  is  used  here. 

The  third  law  values  derived  from  [Kel  71]  are  self-consistent  with  a standard 
deviation  of  0.20  kj/mol,  and  agree  well  with  the  second  law  value.  Thus,  the  value, 
AHqm  =158.02  ± 0.21  kJ/mol,  is  selected  in  this  work.  Figure  3.9  shows  the  evaluated 


40 

Cp(T)  data  along  with  data  from  other  authors.  The  results  of  CP(T)  and  G-Hser(T)  arc 
given  in  Tables  3.10  and  3. 13,  respectively. 


I:  Second  law  value,  HI.  Third  law  value 
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Figure  3 9 Heat  capacity  of  Se?  (g)  and  Ses  (g)  in  the  temperature  range  298  to  2000  K 

3 4 Vapor  Pressure  and  Composition  in  Equilibrium  with  the  Condensed  Phases 
Several  authors  have  investigated  the  vapor  pressure  in  equilibrium  with 
condensed  selenium  using  a variety  of  techniques.  Gronvold  el  ai  [Gro  84]  has  reviewed 
the  data  published  through  1984.  The  equilibrium  partial  pressures  in  equilibrium  with 


Ilie  condensed  phases  are  calculated  using  the  assessed  G-Hsna  expressions  for  the 
condensed  and  gaseous  states  listed  in  Tables  3.11  and  3.13,  respectively.  The  Peng- 
Robinson  equation  of  state,  which  is  commonly  used  for  inorganic  gases,  is  included  to 
calculate  the  fugacity  coefficient  and  compressibility  for  real  gas  corrections.  The 
generalized  form  of  the  Peng-Robinson  equation  of  state  is: 


/» am 

I v-h  V{V*b)*(V-b ) 

With 

ii-  7'- 

0(7)  = 0.45724 — J-o(7) 

P 

* = 0.07780-^- 

K = 0.37464  + 1.5422®  - 0.26992®2 
® = -l.0-log10(pra/’(7>  =0.7  )/pc] 


(3.1) 


(32) 

(33) 
(3.4) 
(35) 
(3.6) 


Here  P'ap  (T,  = 0.7)  is  the  equilibrium  vapor  pressure  of  the  fluid  at  reduced  temperature 
Tt=T/T«=0.7.  The  critical  properties  of  selenium  were  determined  by  Hoshino  el  al  [Hos 
76],  The  ideal  mixture  behavior  is  assumed.  Hence  the  ftigacity  coefficient,  <J>|  (1=1-8), 
and  compressibility,  Zi  0=1-8),  are  equal  for  all  different  species  and  depend  on  only  the 
total  pressure  and  temperature  The  results  for  compressibility  and  fugacity  coefficient  at 
different  temperature  are  listed  at  Table  3.8. 


Table  3.8  Compressibility  and  fugacity  coefficient  for 
selenium  vapor  in  the  temperature  range  =500  to  1500  K. 


Temperature  Compressibility 


The  calculated  partial  pressures  over  solid  and  liquid  selenium  are  given  in 
Figures  3. 10  and  3 .11.  respectively  The  evaluated  total  pressure  in  equilibrium  with  solid 
and  liquid  selenium  with  experimental  data  is  shown  in  Figures  3.12  and  3.13  The 
calculated  vapor  compositions  show  that  Se,.  is  the  major  species  in  equilibrium  with 
solid  selenium  and  Sec.  Scs.  along  with  Se;  are  the  major  species  over  liquid  selenium  up 
to  900  K At  higher  temperature  Se;  is  the  dominant  species.  This  result  is  consistent  with 
the  experimental  reports  by  [Ber  66],  [Ber  68],  [Hor  75],  [Fuj  66],  [Gri  82],  [Hua  84], 
[Dro  84],  [Gol  59],  [Sau  71],  [Gru  74],  and  [Hos  92],  The  calculated  equilibrium  vapor 
compositions  over  solid  and  liquid  selenium  are  given  in  Figures  3.14  and  3 15, 
respectively. 
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Figure  3.10  Equilibrium  vapor  pressure  over  solid  selenium(hexagonal)  in  the 
temperature  range  298-494  K 
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Figure  3.11  Equilibrium  vapor  pressure  over  liquid  selenium  in  ihc  temperature  range 
T=494  to  1500K. 
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Figure  3.12  Experimental  and  assessed  total  vapor  pressure  over  hexagonal  solid 
selenium  at  T=29S  to  494  K. 
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Figure  3.13  Experimental  and  assessed  total  vapor  pressure  over  liquid  selenium  in  the 
temperature  range  494  to  1 500  K. 
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Figure  3.15  Equilibrium  selenium  vapor  composition 
temperature  range  T=494  to  1500  K. 


liquid  selenium  in  the 


Figure  3 16  Calculated  free  evaporation  Se„  (g)  flux  from  liquid  selenium  source. 
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3.5  Flux  Calculation  for  a Freely  F.vaporatinu  Source 
The  thermodynamic  functions  presented  in  this  chapter  were  used  to  estimate  the 
Langmuir  evaporation  flux  These  calculation  are  usefiil  in  characterizing  the 
performance  of  an  evaporative  source.  The  Hertz-Langmuir  equation 

J*.  = oWA.  iJ'hMj.RT  (3.7) 

was  used,  where 

a = evaporation  coefficient 
N,  = Avogardro's  number 
= partial  pressure  of  Se»(g) 

Mfy  = molecular  weight  ofSen(g) 

The  evaporation  coefficient  for  Se»(g)  were  studied  by  Huang  el  aL  [Hua  84)  for  molten 
selenium  at  500  K.  They  reported  the  vaporization  coefficients  of  about  0. 1 for  all  species 
Se»  (g)  (n=  2 to  8).  Their  data  were  re-evaluated  by  Drowart  el  al.  [Dro  84],  who 
determined  that  a=0.26  for  each  of  the  Se„(g)  species  Their  result  was  adopted  in  this 
calculation.  The  temperature  dependence  of  a is  neglected  in  the  calculation  since  no 
data  are  available  for  other  temperature.  The  calculation  results  are  given  in  Figure  3.16. 

3.6  Conclusions 

In  the  present  work,  the  thermodynamic  properties  for  condensed  and  gaseous 
selenium  have  been  critically  evaluated.  Entropy  and  heat  capacity  data  for  gaseous  phase 
species  Sen  (g)  were  calculated  from  molecular  constants.  The  obtained  Gibbs  energy 
functions  were  used  to  perform  third  law  evaluation  for  the  enthalpy  of  formation  value 
of  each  Se„  (g)  (n=3  to  8).  The  enthalpy  of  formation  value  of  Sej  is  critical  assessed 
using  the  selected  equilibrium  vapor  pressure  in  equilibrium  with  liquid  selenium  due  to 
the  disagreement  between  the  literature  data.  This  new  set  of  Gibbs  energy  function  show 


saustactorv  consistency  in  second  and  third  law  evaluations  for  standard  enthalpies  of 
formation  for  Se„(g).(n=3-S)  This  was  a problem  in  previous  evaluation  as  discussed  by 
Drowart  el  at.  [Dro  84],  The  vapor  pressures  obtained  from  equilibrium  calculation  arc 
well  within  the  experimental  results. 


Table  3.9  Assessed  heat  capacity  for  condensed  selenium. 


Phase 

Temperature 

Range (K) 

Coefficients  in  Cp(T)  expression 

C.  = a + bT  + cr 

Solid 

298.15  -494  3 

a = 19.140  b = 2 459  x 10*"  c = -l.606x  10’’ 

494.3  - 760 

a=  19  140  b = 2.459xlff'  c = -l606xlff’ 

760-1200 

a = 28.55 

Liquid 

298.15  -494.3 

a = 52.400  b = -4  985x10'-  c = 3.273  x Iff’ 

494.4  - 1000 

a = 52.400  b = -4  985  x 10"  c = 3.273  x 10'’ 

1000-1150 

a = -5.399  b=  7189x10-  c = -3.121  X Iff’ 

1150-1500 

a = 36.00 

Table  3.10  Assessed  heat  capacity  for  gaseous  selenium. 


Species 

Temperature 

Coefficients  in  Cp(T)  expression 

C,  = a + bT  + cT5  + d/T- 

Se 

298.15  -600 

a = 19  560  b=  1.327x10-  ' c 

3.871  x Iff" 

600  - 2000 

a = 23.230  b=  1.803  x 10°  c 

-7.953  x I0-' 

= -7  967  x 10s 

Se; 

a = 57.000  b=- 4 368x10''  c 

2.878  x I0-’ 

= -4  294  x 10' 

600  - 2000 

a = 34  140  b = 4 271  x10"  ' c 

- 4 885  x Iff' 

= 1.251  x 10“ 

Se, 

298  15  -2000 

a - 58  005  b = 2.221  x Iff4  c 

-6  686  x Iff* 

= -l  966x  10s 

Se, 

298  15  - 2000 

a = 82.578  b = 6.477  x Iff"  c 

-1  950  x Iff' 

= -5  353x  10’ 

298. 1 5 - 2000 

a =107.765  b = 3 668  x Iff'  c 

-1.106  x 10" 

= -4.877  x 10’ 

Sfc 

298  15-2000 

a =132.728  b = 3.425  x Iff"  c 

-1.033 x Iff' 

= -5  785  x 10’ 

Se? 

298  15  -2000 

a =157.585  b = 4.404  x Iff"  c = -1.328  x Iff  ' d = -7.082x!0’ 

Sc* 

298.15  -2000 

a =182  490  b = 4.836 x 10°  c 

-1.459  x 10'7 

= -7  845  x 10’ 

Table  3. 1 1 Gibbs  energy  relative  to  Han  of  condensed  selenium 


CHAPTER  4 

A CRITICAL  ASSESSMENT  OF  THERMODYNAMIC  AND  PHASE  DIAGRAM 
DATA  FOR  THE  CU-SE  BINARY  SYSTEM 


4.1  Introduction 

Cu-Se  system  is  a member  of  the  l-Vl  group  containing  compounds  with  wide 
variety  of  properties.  Cuj.,Se  is  a mixed  conductor,  which  exhibits  both  ionic  and 
electronic  conduction.  It  is  often  present  in  CulnSei-based  solar  cells  as  undesirable 
secondary  phase.  In  contrast,  the  CuSe2  phase  is  a low  temperature  superconductor. 
Knowledge  of  the  thermochemistry  and  phase  diagram  of  the  Cu-Se  system  is  essential 
for  understanding  the  lull  Cu-In-So  system  and  its  subsequent  use  in  developing  new 
process.  A critical  assessment  of  the  Cu-Se  system,  however,  was  lacking  in  the 
literature.  In  this  chapter,  a critical  assessment  of  the  literature  phase  diagram  and 
thermodynamic  data  for  the  Cu-Se  system  was  performed.  The  association  model  was 
used  to  model  the  liquid  phase,  and  a three-sublattice  compound  energy  formalism  was 
used  for  the  Cuj.,Se  phase.  The  remaining  intermediate  solid  phases  (CujSer,  CuSe, 
CuSe-)  were  modeled  as  line  compounds.  A self-consistent  set  of  phase  diagram  and 
thermodynamic  data  was  obtained  through  this  assessment. 

4.2  Phase  Diagram  Data 

The  Cu-Se  liquid  phase  exhibits  two  miscibility  gaps,  one  for  Cu-rich 
composition  and  the  other  for  Se-rich  composition.  Four  intermediate  phases  were 
reported  in  the  literature,  including  C’U'Se  (Cu2.*Se),  CujSej,  CuSe,  and  CuSe2.  The 
CujSe  (Cu;.,Se)  phase  is  the  most  stable  phase  in  this  system  with  a broad  homogeneity 
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range  exlending  towards  ihc  Se  side  to  form  a defect  compound,  Cu;.,Se  has  two 
polytypes:  a-Cu:-,Se  and  p-Cui„Se.  CuSe  has  three  polytypes:  a-CuSe,  P-CuSe,  and  y- 
CuSc.  Extensive  studies  have  been  carried  out  for  the  Cu-Sc  phase  diagram  using  DTA 
(Differential  Thermal  Analysis),  XRD  (X-Ray  Diffraction),  microscopy,  and 
microhardness  techniques.  Earlier  work  covered  the  copper-rich  side  at  high  temperature 
as  reviewed  by  [Han  58],  Later,  [Cha  81]  critically  reviewed  the  phase  diagram  and 
thermodynamic  data  for  the  system,  but  did  not  perform  an  optimized  assessment. 
Subsequent  to  their  review,  [Kor  84]  and  [Abr  84]  further  studied  the  Cu2.»Se  phase  using 
DTA  and  DSC  (Differential  Scanning  Calorimetry)  carried  out  more  studies  for  the  Cui. 
,Se  phase,  [Gla  91]  studied  the  miscibility  gap  on  the  Cu-rich  side  by  the  acoustic 
method.  All  the  phase  diagram  data  are  shown  in  Figures  4.2  a and  b.  The  invariant 
equilibria  are  listed  in  Table  4. 1 1 . 


The  standard  enthalpy  of  formation,  AH°r,29s.tsK,  value  for  a-CuiSe  was  measured 
by  several  authors  using  a variety  of  techniques  as  summarized  in  Table  4.1.  Gattow  and 
Schneider  [Gat  56]  reported  a value  of  AH°i;»*  iss=  -59.3  (kj/mol),  by  direct  synthesis 
calorimeter.  Several  investigators  have  used  solid  state  galvanic  cells  to  measure  the 
Gibbs  energy  of  formation  of  Cu:Se  phase  at  various  temperatures.  Valverde  [Val  68] 
reported  AG'Vi  iik“  -85,0  (kJ/mol)  using  the  following  cell; 


Pt  I C I Pb  I PbCla+KCl  I CuiStH-PbStwCu  I C I Pt 


The  standard  enthalpy  of  formation,  AITV.2ss.isk  =-69.9  (kj/mol),  was  evaluated  by  Mills 
[Mil  74]  from  the  available  literature  data  using  the  third  law  method.  Rau  and  Rabenau 


[Rau  70]  measured  the  equilibrium  Se2  vapor  pressure  in  the  Cu+Cu2Se  two-phase 
domain  and  estimated  the  enthalpy  and  entropy  values  using  the  heat  capacity  data 
reported  in  Kubaschewski'  book  [Kub  58],  Mills  re-evaluated  the  data  using  a third  law 
analysis,  AHVjm  isk  —65.3  (kJ/mol).  and  the  more  accurate  heat  capacity  and  phase 
transition  data  reported  by  Kubascheswki  and  Nulling  [Kub  72].  Askerov  el  al  [Ask  76] 
reported  a value,  AHcimmsk  =-65.7  (kj/ntol)  and  AS°tj9«.isK=28.9  (j/mol),  from  their 
EMF  measurements.  The  value  derived  from  Skeoch  and  Heyding’s  liquid  bismuth 
calorimetric  measurements  is  substantially  lower  than  other  reported  values.  Part  of  this 
discrepancy  is  likely  related  to  the  assumption  that  the  dissolved  liquid  mixture  is  an  ideal 
solution  The  entropy  values,  SWisk.  were  reported  by  Rau  and  Rabenau  [Rau  70]  from 
vapor  pressure  measurement  and  Askerov  el  al.  [Ask  76]  from  EMF  measurements  Their 
data  together  with  the  third  law  re-evaluation  by  Mills  [Mil  74]  are  listed  in  Table  4 7. 
The  heat  capacity  measured  by  Kubascheswki  and  Nolting  [Kub  72]  is  listed  in  Table  4.8. 


4.3.2  B-Cm-.Sc 


The  values  of  the  standard  enthalpy  of  formation  for  p-Cu:Se  as  determined  by 


different  authors  arc  listed  in  Table  4.2.  Gattow  and  Schneider  [Gat  56]  reported  a value, 
AFT fM»i5K  =-55.4  (kJ/mol),  by  direct  synthesis  calorimetry.  Again  the  values  derived 
from  Skeoch  and  Heyding’s  liquid  bismuth  calorimetric  measurements  are  too  low.  A 
value  of  the  standard  entropy  of  p-Cuj..Se  value,  S”  » was  calculated  by  [Rau  70)  from 
vapor  pressure  measurements  Their  data  together  with  [Mil  74]'s  evaluation  value  are 
listed  in  Table  4.7.  The  heat  capacity  measured  by  Kubascheswki  and  Nolting  [Kub  72]  is 
listed  in  Table  4.8,  while  the  values  of  AH””  (o-»P)  from  various  authors  are  listed  in 
Table  4.3 


4.3.3  QjjSsj 


Reported  values  of  the  standard  enthalpy  of  formation,  AHVm  ijk.  of  CujSe:  are 
listed  in  Table  4.4.  Gattow  and  Schneider  [Gat  56]  reported  the  value,  AH'V;,*  1?K  =-98.9 
(kJ/mol),  by  direct  synthesis  calorimetry.  Askerov  el  ai  [Ask  76]  reported  the  value, 
AH”r.K*»K  =-94.6  (kJ/mol),  from  their  EMF  measurements.  The  enthalpy  change  of  the 
peritectoid  reaction,  0.75CujSei(s)->Cui  7sSe(s)+0.5CuSe(s).  was  reported  by  several 
authors  [Hey  66,  Ber  72,  Mur  75]  and  the  values  are  listed  in  Table  4.4.  The  absolute 
entropy  value,  S^Wisk.  was  reported  by  [Ask  76]  from  EMF  measurement.  Their  data 
together  with  [Mil  74]'s  evaluation  value  are  listed  in  Table  4.7 


Reported  values  of  the  standard  enthalpy  of  formation,  Alf  r.29*  isk,  for  CuSe  are 
listed  in  Table  4.5.  Value  of  the  absolute  entropy,  S'WisK.was  reported  by  [Rau  70]  from 
their  vapor  pressure  measurements  and  by  [Ask  76]  from  their  EMF  measurements.  Their 
data  together  with  [Mil  74]'s  evaluated  value  are  listed  in  Table  4.7. 


Author 

dard  heat  of  formation  of  CuSe 

Method 

(kj/mol  > 

[Gat  56| 

Solution  calorimetrv 

[Val  681 • 

EMF 

46  0 

IRau  70|" 

Vapor  pressure 

406 

lAsk  76| 

EMF 

32.6 

[Hey  66]* 

from  reaction  1 75CuSe(s)-> 
Cui.7sSe(s)+0.75Se(l) 

41.8 

[Ber  72]“ 

from  reaction  1 75CuSe(s)  » 
Cui7>Se(s)H).75Se{l) 

42.3 

[Mil  741 

Evaluation 

‘Calculated  in  [Mil  74).  “calculated  in  this  work. 


4.3.5  CuSfcj 

Reported  values  of  the  AFrV.2wt.1sK  from  various  investigators 
4.6,  while  the  reported  values  of  S'ms.isk  are  listed  in  Table  4.7. 


Table  4.6  Standard  heat  of  formation  ofCuSe2  (s)  at  298.15  K. 


•Calculated  in  [Mil  74],  “calculated  in  this  work 


Table  4.7  Standard  eniropy  of  solid  Cu-Se  compounds  al  298  15  K 


Tabic  4.8  Heat  capacity  data  of  CujSe. 


84.2783  - 2.0290 x IQ'1  T 


4.4  Gibbs  Enrev  Dala 

4.4.1  EMF  measurements 

[Mos  89]  invesligaied  ihe  Cu!.,Se  phase,  using  the  solid  slate  galvanic  cell, 
Pl/Cu/CuBr/Cu;.,Se/graphitc,  in  the  temperature  range  350°C  to  450"C  The 
electromotive  force  was  measured  as  (unction  of  the  stoichiometry  of  the  selcnide 
compound  and  the  temperature.  The  composition  of  the  solid  solution  (Cuj.,Se)  was 
varied  by  coulometric  titration.  The  chemical  potential  of  copper,  pen.  in  Cu^Se  is  given 
by. 

Mcu-  M°cu  + EF 
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where  F is  the  Faraday  constant,  E is  the  open  circuit  potential,  and  u'Y-u  is  the  chemical 
potential  of  pure  Cu.  The  measured  values  are  plotted  in  Figure  4.4. 

4.4.2  Vapor  Pressure  Measurement 

The  following  gas  species  were  reported  to  exist  in  the  Cu-Se  system,  Se„  (n=l- 
8),  Cu.  Cui,  CuSe,  and  CujSe,  with  Se;  the  most  prominent  one  Rau  and  Rabenau  [Rau 
70]  studied  the  Sc.  equilibrium  pressures  for  Cu-Se  system  as  shown  in  Figure  4.2  The 
partial  pressure  data  have  been  converted  into  partial  Gibbs  energy  data  and  used  in  the 
assessment.  Blachnik  and  Bolle  [Bla  78]  determined  the  activities  of  selenium  in  molten 
mixture  of  copper  and  selenium  by  a modified  dew-point  method  at  1373  K.  Azakami 
and  Yazawa  |Aza  76]  measured  the  selenium  activities  in  copper-selenium  melt  at 
I200°C  by  the  transportation  method. 


4.5  Thermodynamic  Models 


The  Gibbs  energy  of  Cu  is  from  [Din  91],  and  that  of  Se  is  from  [Cha  99],  The 
three-term  equation  given  below  is  used  to  represent  the  temperature  dependence  of  the 
Gibbs  energies  of  stoichiometric  compound  phases, 


°G=a+bT+cTlnT 


Where  °G  is  the  standard  Gibbs  energy,  T is  the  absolute  temperature,  and  a,  b and  c . 
constants  whose  values  need  to  be  determined  from  experimental  data. 


A limited  solubility  of  Se  in  Cu  was  measured  by  [Sma  46,  Tay  76],  while  negligible 
solubility  ofCu  in  Se  was  proposed  by  [Cha  81],  In  this  work,  the  solubility  of  Cu  in  Se 


was  also  considered  negligible,  whereas  the  solubility  of  Se  in  Cu  is  modeled.  The  Gibbs 


energy  (Grcc)  of  the  terminal  solid  solution  based  on  face-centered  cubic  (fee)  Cu  is 
expressed  as 

(4.2) 

",G*'=xCm-G&+xs.'G£  ,4.3) 

‘'Gj6'  = «r(xt„lnxc,tx&lnxs,)  (4.4) 

tGft;  is  the  excess  Gibbs  energy,  and  expressed  by  the  Redlich-Kistcr  polynomial  (Red 
48]: 


CG/*=*c. At,  I '/£&(*■»-• 


where  L is  the  binary  interaction  parameter  to  be  optimized  in  the  present  work.  The 
temperature  dependence  of  L may  be  expressed  as 


aai 


(4.6) 


4.5.3  Liquid  Phase 


The  association  model  developed  by  [Som  82]  was  used  to  describe  the  Gibbs 
energy  of  the  liquid  phase.  Cu2.,Se  is  the  most  stable  compound  in  this  system,  which  is 
stable  to  1403  K and  congruently  melts  at  that  temperature,  while  the  other  compounds 
peritcctically  decompose  below  650  K It  is  thus  reasonable  to  select  CujSe  as  an 
associate  species  for  the  modeling  of  the  liquid  phase 

(4.7) 


(Cu,  CujSe.  Se)i 


The  Gibbs  energy  of  this  phase  (per  mole  of  atoms)  can  also  be  expressed  by 
equation  (4.2),  where  fee  should  be  replaced  by  I (Miquid).  The  term  "G1  is  equal  to 

L (4.8) 

where  yj  refers  to  the  site  fractions  of  the  species  i,  >’c.  * .V'cw,  * > : 1 The  terms  “c/ 
represent  the  Gibbs  energies  of  the  pure  liquid  phase  of  species  i. 

The  ideal  mixing  Gibbs  energy  "'G1  is  equal  to 

“G*  =«T(yol  lnyc.  + In  Iny,..).  ^ 

and  the  excess  Gibbs  energy  tGl  is  by  Redlich-Kister 

(4. 10) 

where  the  four  L terms  represent  the  interactions  between  the  indicated  species.  They  can 
be  expressed  as  the  junction  of  temperature. 

The  high  temperature  modification  of  Cuj.*Se,  |)-Cuj.,Se,  is  related  to  the  B4  (ZnO) 
or  CI(CaFj)  fee  structure.  [Ral  36,  Bor  45,  Ste  71]  suggested  that  p-Cua.,Se  is  the  ZnO- 
type  structure  with  the  space  group  /•'43m . In  this  structure,  four  Se  atoms  occupy  the 
4(a)  sites,  four  Cu  atoms  occupy  the  4(c)  sites,  and  the  other  Cu  atoms  are  statistically 
distributed  in  the  interstitial  sites  (four  tetrahedral  (d),  four  octahedral  (b),  and  sixteen 
trigonal  (e)  positions).  [Hey  76]  and  the  later  investigators  [Ton  8],  Oli  88.  Sak  89,  Yam 
91]  suggested  another  structure  model  based  on  the  space  group  Fm3m.  Each  of  these 
investigators  agreed  that  the  Se  atoms  occupy  the  4(a)  sites  and  form  a fee  sublattice 
While  [Oli  88]  proposed  that  all  Cu  atoms  occupy  the  trigonal  32(f)  sites,  and 


investigators  [Ton  81,  Sak  89,  Yam  91)  proposed  that  Cu  atoms  lie  on  two  different  sites. 
Based  on  the  structural  data,  the  ordered  nonstoichiometric  p-Cu;.,Se  phase  was  divided 
into  three  sublattices  one  is  occupied  by  Se  atom,  and  the  other  two  occupied  by  Cu 
atoms.  The  structural  data  suggested  the  nonstoichiometry  is  due  to  the  formation  of 
vacancies  on  the  Cu  sublattice.  The  energetics  of  the  two  Cu-sublattices  are  different.  The 
Cu  atoms  can  lie  more  easily  on  one  sublattice  with  lower  energy  than  the  other  one. 
Therefore,  vacancies  form  more  easily  on  the  sublattice  with  higher  energy,  and  thus 
assume  one  of  the  two  Cu  sublattices  is  always  fully  occupied  by  Cu  atoms.  As  a result, 
the  P-Cuj.«Se  phase  was  described  using  the  sublattice  model  developed  by  [Hil  70,  Sun 
81]  with  three  sublattices  after  the  formula: 

(Cu.Va)t(Se,Va)i(Cu)i  (4  11) 

where  Va  represents  vacancy.  Including  the  possibility  of  vacancies  in  Se  sublattice 
introduces  flexibility  in  the  modeling  of  the  Gibbs  energy,  but  does  not  imply  that  the  Se 
sublattice  contains  a high  concentration  of  vacancies  From  the  assessment  results  given 
later.  It  is  found  find  that  the  vacancy  concentration  in  the  Sc  sublattice  is  very  small, 
which  is  in  good  agreement  with  the  structural  information. 

The  low-temperature  modification  of  Cu:  ,Se.  a-Cu:  ,Se,  exhibits  a complex 
diffraction  pattern  and  makes  it  difficult  to  identify.  The  structure  of  a-Cu;.,Se  has  been 
described  as  cubic  [Ral  36,  Sor  74,  Yam  91 J,  pseudo-cubic  [Med  67],  tetragonal  [Bor  45, 
Asa  72],  orthorhombic  [Mar  69.  Ste  69,  Ste  71],  monoclinic  [Vuc  81,  Mil  87,  Fra  91], 
and  pseudo-monodinic  [Kas  88],  Each  investigator,  however,  agrees  that  the  Se  atoms 
form  a fee  sublattice  and  the  Cu  atoms  distribute  over  two  or  more  various  sites  (e  g , 
tetrahedral,  octahedral,  and  trigonal  sites).  For  modeling  the  homogeneity  range,  the 


ordered  nonstoichiomeiric  a-Cuj.,Se  phase  was  also  described  with  three  subiattices  after 


formula  (4.7).  The  Gibbs  energy  of  such  a phase  4>  (4=a-Cu2.*Se  or  P-Cu2-*Se)  can  also  be 
expressed  by  equation  (4.2)  (where  fee  is  replaced  by  $)  with. 

* ycJv. ‘G*^ 

yyn/s.'G^sM,  (4.12) 

UG*  »/W|Cv'cs  In  Vo  ♦ Vlv  In  v.. ) 

In /*,+>•  V.  In  vr, )|  m 13) 

SG'=J'a)’r.[vL(,’4ww»+,4j^a’a -*■'  ))*  y'v.  ("^erw-rx,  ♦'4xW«sC>ti.  ->V„>8 
+ y'uiy. [fa, +‘4wwws, - Vi't. ))+  y'v. (’4«ww* +'4r*(k& 0's,  - /r. )) 

(4.14) 

where  y,  and  y,  refer  to  the  site  fractions  of  the  component  i in  the  first  and  second 
sublattice,  respectively.  ’G*,*,.,  is  the  standard  Gibbs  energy  of  stoichiometric  a-Cu2Se 
or  p-CuzSe,  which  will  be  assessed  in  this  work.  ‘G*,  and  ._  can  be  expressed 

'<%*■*.  =rG£+°,+b,T  (4)5) 

(4I6) 

where,  is  the  standard  Gibbs  energy  of  pure  copper  in  the  fee  phase,  and  a,  and  bi 
are  the  optimized  parameters.  The  value  of  'G*JilC_  can  be  expressed  by  the  following 
relation  [Ans  85,  Zhu  99]: 

i’GC.iS.t.  +‘G*«  J',;C.-,C&;K.:CV  (4.17) 

To  reduce  the  number  of  parameters,  the  following  relations  arc  used  in  the  optimization 
[Ans  85,  Zhu  99]: 


(4.18) 
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'^-r.!j.<.ru  -"ihyv-c* 


’ 7?-*  t*~  l$asty*e*=  L 


A selected  set  of  thermodynamic  and  phase  diagram  data  (listed  in  Table  4.9)  were 
used  for  the  optimization  of  thermodynamic  model  parameters  of  all  phases  in  this 
system.  The  critical  review  of(Cha  81]  was  used  as  a guide  for  selecting  the  experimental 
phase  diagram  data.  New  experimental  phase  diagram  data  [Abr  84,  Kor  84,  Gla  91,  Pan 
99)  were  also  used  in  the  optimization.  The  optimization  was  performed  by  using  the 
PARROT  module  of  the  Thermo-Calc  program  package  [Sun  85]  and  split  in  six  steps 
First,  the  Gibbs  energy  coefficients  of  the  stoichiometric  compounds,  a-Cu-Se  and 
p-CujSe  and  the  stoichiometric  liquid  CujSe  were  estimated  from  enthalpy  of  formation, 
standard  entropy,  heat  capacity,  transition  enthalpy,  transition  temperature,  and  melting 
enthalpy  and  temperature  Then,  assuming  that  p-Cu2„Se  is  a stoichiometric  phase  and 
there  is  no  interaction  between  the  pure  elements  in  the  liquid  phase,  a preliminary 
optimization  of  the  liquid  coefficients  was  performed  based  on  the  activities  of  selenium, 
miscibility  gap,  and  three-phase  equilibrium.  Third,  fixing  the  coefficients  of  the  liquid 
phase,  the  coefficients  of  the  ordered  nonstoichiometric  P-Cu!.«Se  phase  were  roughly 
optimized  based  on  the  measurements  of  the  activities  of  copper  [Mos  99],  two-phase 
equilibria  between  liquid  and  p-Cu2,,Se,  and  the  monotcctic  and  eutectic  equilibria  data 
Fourth,  fixing  the  coefficients  of  the  liquid  and  P-Cui.»Se  phases,  the  coefficients  of  the 


stoichiometric  CujSej,  oc-CuSe,  P-CuSe  (y-CuSe),  and  CuSej  phases  were  optimized 
based  on  their  enthalpies  of  formation,  enthalpies  of  transition,  standard  entropies,  and 
equilibria  data  that  included  these  phases.  As  no  discontinuous  change  in  enthalpy  was 
detected  in  the  P-CuSes-*y-CuSe  transformation  by  DTA  method  (Mur  75].  the 
transformation  may  be  considered  to  be  continuous,  i.e.,  higher  order  [Cha  81],  In  this 
work,  the  P-CuSe  and  y-CuSe  phases  are  assumed  and  modeled  as  one  phase.  Fifth, 
fixing  the  coefficients  of  other  phases,  the  coefficients  of  the  ct-Cu!.,Se  and  terminal 
solid  solution  fcc(Cu)  phases  were  optimized  based  on  available  experimental  data. 
Finally,  all  parameters  of  the  phases  were  optimized  to  fit  the  experimental  data  listed  in 
Table  4 9. 


4.7  Result  and  Discussion 

The  optimized  parameters  of  the  stable  phases  in  the  Cu-Se  system  were  listed  in 
Table  4.10.  The  parameters  taken  from  [Din  91]  and  [Cha  99]  were  not  reproduced  in 
Table  4. 10  to  conserve  space.  The  phase  diagram  and  thermodynamic  properties  of  this 
system  have  been  calculated  by  using  the  optimized  parameters,  which  were  performed 
with  the  Poly-3  module  of  the  Thermo-Calc  program  package. 

The  calculated  phase  diagram  is  shown  in  Figure  4.1.  Tile  dash  line,  which  denotes 
the  [i-CuSewy-CuSe  - transformation  is  not  the  calculated  result.  It  is  based  on  the  XRD 
measurement  of  [Mur  75].  Figures  4.2a-c  compared  the  calculated  phase  diagram  and  the 
measured  data.  The  calculated  phase  diagram  agrees  well  with  the  data  used  in  the 
optimization.  It  does  not  agree  with  some  of  the  data  which  was  not  used  in  the 
optimization,  including,  the  Cu-rich  miscibility  gap  data  measured  by  [Bur  74]  (Figure 


4.2a),  and  the  a-Cui.*Se/p-Cu2*Se  solvus  line  and  related  eutectoid  measured  by  [Ogo 
69,  Ogo  72]  (Figure  4.2b).  These  data  are  inconsistent  with  the  new  measured  results. 
Specifically,  the  Cu-rich  miscibility  gap  data  measured  by  [Bur  74]  arc  disagreed  with 
those  measured  by  [Aza  76,  Gla  91];  while  the  a-Cu;.,Se/[l-Cu;.,Se  equilibria  and  related 
eutectoid  measured  by  [Ogo  69,  Ogo  72]  are  inconsistent  with  those  measured  by  [Abr 
84], 

Table  4. 1 1 presents  the  experimental  and  calculated  temperatures  and  compositions 
of  the  invariant  reactions  in  this  system.  The  calculated  values  agree  well  with  the  data 
summarized  by  [Cha  81].  The  calculated  critical  temperature  of  the  Cu-rich  liquid 
miscibility  gap  agrees  well  with  the  experimental  data  of  [Gla  91],  whereas  the  calculated 
critical  composition  has  lower  ae  concentration  than  the  experimental  one  [Gla  91], 
However,  the  calculated  compositions  of  the  whole  Cu-rich  miscibility  gap  agree  well 
with  the  experimental  data  of  [Gla  91]  as  shown  in  Figure  4.2a.  In  the  calculated  phase 
diagram,  the  temperatures  and  compositions  of  the  Sc-rich  liquid  miscibility  gap  have 
been  predicted,  which  need  to  be  verified  by  furture  experiments. 

Figure  4.3  shows  the  comparison  between  the  calculated  and  measured  chemical 
potential  of  selenium  in  liquid  phase  The  calculated  values  agree  with  the  measured  ones 
[Bla  78]  when  the  mole  fraction  of  selenium  is  larger  than  0.333.  When  the  mole  fraction 
of  selenium  is  less  than  0.333,  the  calculated  values  lie  between  the  measured  results  of 
[Aza  76]  and  [Bla  78].  The  calculated  and  measured  chemical  potential  of  copper  in  the 
P-Cu!.,Se  phase  is  shown  in  Figure  4.4.  The  calculated  curves  essentially  agree  with  the 
experimental  data.  The  calculated  curves,  however,  are  not  as  concave  as  the  measured 
ones  Figure  4.5  shows  the  comparison  between  the  calculated  and  measured  Se2  partial 
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pressure  along  the  liquids  and  in  equilibrium  with  the  compounds  in  this  system.  To 
calculate  the  Se-  partial  pressure,  an  ideal  gas  phase  containing  12  species  is  added  to  the 
optimized  coefficient  set.  The  parameters  of  the  gas  species  Se,  S(.-:.  Sei,  Sea,  Se-.  See, 
Se7.  and  Sex  are  given  in  chapter  3,  whereas  the  parameters  of  the  gas  species  Cu,  Cu:, 
CuSe,  and  Cu;Se  are  taken  from  [Sun  94],  The  calculated  Se2  partial  pressure  values 
agree  well  with  the  measured  data  of  [Rau  70]  except  those  at  the  low  pressure  fcc(Cu)/[l- 
Cu2.,Sc  region.  Attempt  to  fit  will  not  allow  other  experimental  data  to  be  fitted, 
including  the  experimental  Se:  partial  pressures  measured  by  the  same  authors  in  [Rau 
70], 

The  calculated  vacancy  site  fractions  in  the  fi-Cu:.,Se  phase  are  shown  in  Figures 
4.6  and  4.7.  In  the  (l-Cu:.,Se  single-phase  region,  the  Cu  vacancy  concentration  in  the 
first  sublattice.  yVa , increases  with  the  increase  of  selenium  concentration;  whereas  the  Se 
vacancy  concentration,  yVa , decreases  with  the  increase  of  selenium  concentration.  The 
horizontal  lines  from  the  mole  fraction  of  selenium  at  0.338  to  0.342  at  380  1C  in  Figures 
4.6  and  4.7  correspond  to  the  a-Cu:.»Se+p-Cu2.,Se  two-phase  region.  The  calculated  Cu 
vacancy  site  fraction  in  the  first  sublattice  is  between  0.008  to  0.271  in  Figure  4.6, 
whereas  the  Se  vacancy  site  fraction  is  between  2.07 *I0*17  to  2.23  xIO*4  in  Figure  4.7.  No 
experimental  data  are  available  for  comparison  to  these  calculation  results.  However,  the 
very  limited  Se  vacancy  and  a large  Cu  vacancy  concentration  are  in  reasonable  accord  to 
the  structural  research  results  [Ste  71,  Hey  76,  Yam  91], 

The  experimental  and  calculated  standard  enthalpies  of  formation  of  the 
intermediate  compounds  at  298. 1 5 K in  the  Cu-Se  system  are  presented  in  Table  4. 1 2 The 
corresponding  standard  entropies  are  presented  in  Table  4.13.  The  reference  state  for 
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these  properties  is  the  enthalpy  of  the  pure  stable  element  at  298. 1 5 K.  The  calculated 
enthalpies  of  formation  in  the  present  work  lie  between  the  data  measured  by  [Gat  56] 
and  those  evaluated  by  [Mil  74],  The  calculated  standard  entropies  agree  with  those 
evaluated  by  [Mil  74],  except  that  ofCuSej,  which  lies  between  the  value  of  [Rau  70]  and 
that  evaluated  by  [Mil  74],  Figure  4.8  shows  the  calculated  heat  capacity  for  Cu;Se  along 
with  the  experimental  result  [Kub  73],  The  calculated  values  agree  well  with  the 
measured  ones  The  calculated  enthalpy  of  transition  of  the  a-Cu;Se->[i-Cu;Se  transition 
at  396  K is  6.83  kj/mole,  which  is  in  good  agreement  with  that  measured  by  [Kub  73], 
The  calculated  enthalpy  of  the  a-CuSe->P-CuSe  transition  at  324  K is  1.38  kj/mole, 
which  is  in  good  agreement  with  that  measured  by  [Hey  66],  As  no  experimental  heat 
capacity  data  for  CuSe  are  available,  those  data  are  optimized  from  the  phase  diagram 
data  after  assuming  no  difference  between  the  heat  capacities  of  a-CuSe  and  p-CuSe. 

4.8  Conclusions 

A thermodynamic  description  of  the  Cu-Se  binary  system  was  obtained  by 
optimization  of  the  available  phase  equilibrium  and  thermodynamic  data  The  Redlich- 
Kister  polynomial,  associate  solution  model,  and  sublattice  model  were  used  to  represent 
the  Gibbs  energy  of  the  fcc(Cu),  liquid,  and  both  Cu;.,Sc  phases,  respectively.  Other 
intermediate  phases  are  modeled  as  line  compounds.  The  phase  diagram  and 
thermodynamic  properties  of  this  system  have  been  calculated  by  using  the  optimized 
model  parameters  There  is  reasonable  agreement  between  the  model-calculated  values 
and  selected  phase  equilibrium  and  thermodynamic  data  available  in  the  literature  for  all 
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phases.  Further  experimental  work  is  suggested  for  determining  the  shape  and  critical 
point  of  the  Se-rich  liquid  miscibility  gap. 


Table  4. 10  Optimiz 
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Phase  or 

*arameters 

Linuid 

aCu..,Se 

a.='cnT  +90000 

'OJja  = 40(1011  * GCu2Se  _ A-'G[“ 

“'•Lav.  =11180+10/’ 

I$.r.-c  =-5678 

)Cu.,Se 

C&J.O.  =OCu2Se_B 

GiL-M.-G£  +80000+  36/' 

liana.  -20OO0 

ll'Sc’: 

3 '"=-25590+62  7S2UT-I5 ,66455TlnT 

xC'nSe 

3"“'=-25859+83.0760T-l8  24325TlnT 

jf‘+’c=-2447S+78X145T.|8  24325TlnT 

IlSe. 

J - "=-2 1489*95  3 1 S»T- 1 9 977  lOTInT 

ec(Cu) 

EST*  ">  'lus-m 

3cu2Se_A=-!(U2Ji34+28S  l6728T-59.0572TlnT-0  0375096T’  (298<T<395) 

/=-9S588.35+664.34fi7IT-1200866TlnT+0.03785T: 

— 6:9635 * 1 (r^T'+ 1 0 1 9900T ' (395<T<800) 

3cu2Sc  B=GCu2Sc  A+6830-I7.29I  I4T 

3Cu2Sc  l=GCu2Se  B+I6000-1 1 422T 

P-CuSc  and  CuSc;  are  in  J/mole  of  atoms  The  Gibbs  energies  ofa-tu.  ,Sc  and  p-Cu - ,Se  arc  in 
J/molc  of  (Cu.Va),(Sc.Va)i(Cu)i.  The  symbol  • indicales  Cu,  Sc  or  Va.  and  lrl  represents 
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reactions  in  the  Cu-Se  system.  P P 
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Table  4. 12  Experimental  and  calculated  standard  enthalpies  of  formation  (AH°  um.uk)  of 


Table  4. 1 3 Experimental  and  calculated  standard  entropies  <"S»s)  of  the  intermediate 
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Figure  4 IThe  calculated  Cu-Sc  binary  phase  diagram  based  on  the  optimized  parameters 
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Figure  4.2a  Comparison  belween  the  calculated  Cu-Se  phase  diagram  and  various 
experimental  data. 
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Figure  4.2b  Comparison  between  the  calculated  Cu-Se  phase  diagram  and  various 
experimental  data  available  in  the  literature  in  the  vicinity  of  u-Cu;.,Se  and 
P-Cu2.*Se  equilibria  and  the  related  culccloid. 


Mole  fraction  of  Se 


Figure  4.2c  Comparison  between  the  calculated  Cu-Sc  phase  diagram  and  experimental 
data  available  in  the  literature  for  the  Se  solubility  in  Cu  terminal. 
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Figure  4.3  Calculated  chemical  potential  of  Se  in  the  liquid  phase  along  with  the 

experimental  data  [76Aza,  78Bla],  The  reference  state  is  the  liquid  phase  of  pure 
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Figure  4.4  Calculated  chemical  potential  of  Cu  in  the  (l-Cu;.,Sc  phase  with  the 
experimental  data  [89Mos],  The  reference  state  is  the  fee  phase  of  pure  a 
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Figure  4.5  Calculated  Sej  partial  pressures  in  the  Cu-Se  system  (atm)  along  with  the 
experimental  data  [70Rau]. 
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Figure  4.6  The  calculated  site  fractions  of  Cu  vacancy  on  the  first  sublattice  for  the 
0-Cu:.1Se  phase. 
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Figure  4.7  The  calculated  site  fractions  of  Se  vacancy  on  the  second  sublattice  for  the 
[S-Cu..,Se  phase. 
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Figure  4 8 Calculated  heal  capacity  for  a-Cu2Se  and  0-CujSe  along  with  Ihe 
experimental  values  [73Kub], 


CHAPTERS 

A CRITICAL  ASSESSMENT  OF  THERMODYNAMIC  AND  PHASE  DIAGRAM 
DATA  FOR  THE  IN-SE  BINARY  SYSTEM 


5.1  Introduction 

The  In-Se  system  is  a member  of  the  III-VI  group  containing  semiconducting 
compounds  (InsSes,  InSe,  IncSe?,  IihtSes)  with  possible  applications  in  electronic  devices. 
Knowledge  of  the  thermochemistry  and  phase  diagram  In-Se  system  is  necessary  for  this 
study  of  the  Cu-In-Se  system.  A critical  assessment  of  the  In-Se  system,  however,  is 
lacking  in  the  literature.  In  this  work  an  assessment  of  the  In-Se  system  using  the  phase 
diagram  and  thermochemical  data  in  the  literature,  the  SGTE  unary  assessment  [Din  91) 
for  indium,  and  the  one  for  selenium  given  in  chapter  3.  Nine  intermediate  solid  phases 
(ImSej,  InSe,  lifeSe,,  InsSeM,  ImSer,  a-ln:Se>,  p-  InjSej,  y-  InjSe, , and  5-In;Se,)  are 
modeled  as  line  compounds,  The  associate  model  [Som  82)  was  used  to  describe  the 
liquid  phase.  By  coupling  the  thermochemical  and  the  phase  diagram  data,  optimized  G- 
Hssa  expressions  for  the  solid  compounds  and  the  liquid  phase  are  obtained  together  with 
the  major  vapor  phase  species  (Se„(n=I  to  8),  In,  InSe,  InjSe,  InjSej)  in  In-Se  system. 
The  T-x  and  P-T  projections  of  the  P-T-x  diagram  were  assessed. 


5.2  Phase  Diagram  Dala 


The  phase  diagram  of  In-Se  system  was  studied  by  ['Sla  63,  !Sln  63,  *Sla  63,  Gul 
65,  Lik  74.  Ima  81.  Gla  89,  Did  96.  Dao  98,  and  G6d  98J.  The  early  studies  by  ['Sla  63, 
JSla  63,  'Sla  63,  and  Gul  65J  were  performed  by  DTA  and  X-ray  analysis.  Four 
intermediate  compounds.  ImSe(s),  InSe(s),  In-Sc^fs),  and  IniSe-(s)  were  reported.  InSe(s) 
and  InjSeifs)  were  identified  as  congruent  melting  species.  One  miscibility  gap  was 
proposed  in  their  diagram.  Likforman  and  Guittard  [Lik  74]  investigated  the  In-Se  system 
from  10  to  95  at.%  Se  using  DTA  and  X-ray  analysis.  They  reported  five  compounds 
liuSes  (s)  (instead  of  ln:Se).  InSe  (s),  liv.Se?  (s)  (instead  of  In-Se,).  In-Se-  (s),  and  In-Se- 
ts) Only  In-Se-  (s)  was  identified  as  a congruent  melting  species  and  two  misicibility 
gaps  were  proposed  in  their  diagram.  Imai  el  al  [Ima  8 1 ] investigated  the  phase  diagram 
from  30  to  56  at.%  by  DTA.  They  also  carried  out  the  composition  analysis  of  the 
compounds  grown  from  the  stoichiometric  and  non-stoichiometric  melts  using  He'  ion 
Rutherford  backscniienns  techniques  (RBS).  Their  diagram  was  in  qualitative  agreement 
with  the  diagram  reported  by  Likforman  el  al,  [Lik  74].  The  main  differences  are  the 
compositions  of  the  peritetic  reactions  and  the  shape  of  the  liquidus  lines. 

Glazov  el  al  [Gla  89]  measured  the  ultrasound  propagation  rate  at  different 
composition  and  temperature  of  liquid  In-Se  alloys  to  study  the  In-rich  misicibility  gap 
Okada  and  Ohno  [Oka  93]  investigated  the  electronic  properties  of  liquid  In-Se  alloys, 
including  electrical  conductivity,  thermopower,  and  magnetic  susceptibility.  They  also 
determined  the  melting  temperature  of  In-Se  alloys  at  different  compositions  from  the 
temperature  dependence  of  the  electronic  properties.  The  In-Se  phase  diagram  was 
reviewed  by  Okamoto  [Oka  96]  based  on  the  work  of  ['"  Sla  63,  Lik  74,  and  Ima  81], 


Didoukh  [Did  96]  measured  the  electroconductivuy  for  liquid  immiscible  Ini.Se,  alloys 
and  determined  the  miscibility  gap  on  the  In-rich  side  of  the  phase  diagram.  Daouchi  cl 
at.  [Dao  98)  reinvestigated  the  In-Se  phase  diagram  using  DTA,  DSC,  and  XRD.  Their 
diagram  is  in  qualitative  agreement  with  [Oka  96],  Gildecke  el  al  [God  98]  re- 
determined the  stable  and  metastable  In-Se  phase  diagram  using  DTA,  XRD,  optical 
microscopy,  SEM,  and  TEM.  They  reported  ln*>Sen  and  InsSer  are  stable  phases  at 
stoichiometric  composition  and  P-InjSe*  was  observed  at  59.6  at.  % contrary  to  previous 
studies,  which  suggested  p-ln2Sej  was  one  of  the  polymorphic  stoichiometric  ImSej 
phases.  They  also  suggested  a-InjSe.1  phase  is  slightly  Se-rich  compared  to 
stoichiometric  y-InjSc)  phase,  however,  the  exact  composition  was  not  reported  The 
formation  of  InsSe?  and  InySeu  was  suppressed  while  applying  cooling  rate  between  2 to 
10  K/min.  They  attributed  the  difference  between  their  Ending  and  the  literature  data  to 
the  different  experimental  procedures.  The  phase  diagram  data  from  various  authors  are 
shown  in  Figure  5.2.  The  invariant  equilibria  are  listed  in  Table  5.1. 


5.3,1  Unary  dala 

Gibbs  energy  data  of  condensed  phase  indium  was  taken  from  the  SGTE  database 
[Din  91 J.  The  thermodynamic  properties  of  pure  selenium  were  assessed  in  this  work  as 
discussed  in  chapter  3 The  data  are  given  in  Tables  5 4 and  5.5  in  the  form  of  (G-Hsua) 


Table  5.4  Gibbs  energy  relative  to  Hum  of  condensed  indium. 


Phase 

Range  IK) 

Coefficients  in  (G-H-j., ) expression 

G-Hser  = A * BT  + CTInT  + DT'+Er+FfT+GT'  * HT" 

Hexagonal 

298.15-429.75 

A ■-6978  89  B-  92.33381 15  C ■ - 21 .8.386  D - - 5 .72566  x 10  ' 

429.75  - 3800. 

A ■ -7033.516  B ■ 124.476588  C = - 27.4562  D - -0.54607  x It)  ' 

Uqutd 

298.  U -429.75 

A = - 3696.798  B=  84.701255  C> -21.8386  D ■- 5.72566  x 10  ’ 

494.4-3800 

A - - 3749.81  B - 1 16  835784  C - - 27  4562  D ■ 0.54607  x 10  ' 

E ■ - 0 08367  x 10 *‘  F ■ -211708  G • 0 H - u 

Note  *-  superheated  solid.  *•-  supercooled  liquid 


5.3.2  !bjSsj(s) 

In;Se»(s)  was  reporled  lo  exist  in  several  polymorphic  phases,  however,  the 
literature  contains  contradictory  results  on  the  structure  and  phase  transition  temperature 
[Sla  63.  Med  65,  Lan  74,  Lik  78,  Pop  71,  Pop  79,  Lik  80,  Man  80,  Oka  91,  Dao  98,  and 
God  98],  The  uncertainty  about  the  structure  of  InjSej  is  likely  related  to  several  reasons, 
including:  I . Complexity  and  similarity  of  the  different  polymorphic  structures  that  exist 
in  the  vicinity  of  this  composition.  2.  Different  techniques  used  to  prepare  the  samples.  3. 
The  volatile  nature  of  taSe)  (s)  caused  the  composition  to  drill  during  the  heat  treatment 
involved  in  phase  transition  experiments.  4.  The  difficulty  in  performing  in-situ  structural 
invesugauons  without  losing  volatile  components  from  the  vapor  phase. 

Mills  [Mil  76)  reported  the  only  information  concerning  the  high  temperature  heat 
capacity  and  enthalpy  of  transition.  He  reported  the  heat  capacity  for  ln:Ser(s)  and 


IniSei.Js)  in  Ihc  temperature  range  T=298  to700  K and  the  iransilion  enthalpy  at  486  K 
The  XRD  analysis  on  the  quenched  in  In;Sei.,(s)  sample  after  the  DSC  measurements 
indicate  this  In-rich  InrSe,  is  (l-ln^Sev  This  result  is  consistent  with  [Clod  98]‘s  phase 
diagram.  Mill's  Cp  data  were  used  to  tit  the  heat  capacity,  Cp(T),  equations  for  at,  0,  and 
y-Iil;Sc)  The  £-In2Scr(s)  heat  capacity  was  estimated  by  [Pas  88]  from  In2Sj  heat 
capacity  data.  The  fitted  Cp(T)  results  were  given  in  Table  5.8. 

The  standard  enthalpy  of  formation,  AH'r,2wr  isk.  for  a-InrSej  has  been  studied  by 
several  authors  using  different  techniques,  including  combustion  calorimetry  [Han  53], 
synthesis  calorimetry  [Sha  67],  mass-spectrometry  [Ber  69],  vapor  pressure  (absorption 
spectra)  [Gri  72],  EMF  measurements  [Mus  74],  [Mai  99],  vapor  pressure  measurements 
(Knudsen  or  torsion)  [Ven  79,  Pas  88,  Bar  88],  and  analysis  of  the  reaction:  In2Sej(s)+ 
l!(g)=2Inl(g)+l.5Sej(g)by  [Gri  75]. 

Mills  [Mil  76]  reported  the  value  for  a-InjSei  to  fi-In-Sej  transition  enthalpy, 
AH«=3.85  kJ/mol,  using  DSC.  His  value  was  re-interpreted  as  ct-In-Sej  to  y-ln2Sej 
transition  based  on  the  recent  phase  diagram  [G9d  98],  The  standard  enthalpy  of 
formation,  AH°r298.isK=  -319.9  kJ/mol,  for  y-InjSer  was  estimated  from  the  standard 
enthalpy  of  formation,  AH°r.M«i>K=  -318  kJ/mol  [Sha  67],  for  a-lnjScj  using  this 
transition  enthalpy. 

The  vapor  pressure  measurements  by  [Gri  72,  Gri  75,  Ven  79,  Pas  88,  Bar  88] 
were  performed  in  high  temperature  range,  where  S-ln2Se2  and  y-ln2Sej  are  the  stable 
phases  The  vapor  pressure  data  reported  by  [Ven  79,  Pas  88,  Bar  88]  were  used  to  derive 
the  standard  enthalpy  of  formation  value  for  6-In2Sej.  The  results  are  summarized  in 


Table  5.6. 
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Srinivasa  and  Edwards  [Sri  87]  studied  the  vaporization  chemistty  by  the 
simultaneous  Knudsen  and  torsion  method  in  the  vicinity  of  stoichiometric  In-Sei  If  we 
re-interpret  their  data  based  on  the  new  diagram  by  [God  98).  The  data  they  interpreted  as 
the  a-lnaSe3+InsSe6  two-phase  equilibrium  vapor  pressure  could  be  re-interpreted  as  the 
T-P  projection  above  the  y-InjScj  and  (1-In.Se,  two-phase  domain.  Their  vapor  pressure 
data  were  used  to  estimate  the  value  for  p-InjSej  according  to  the  following  vaporization 

y-InjSej— ♦P-ItijSea.sj  0.02475Se2 

The  value.  AH°f.29*irKe  - 3 14.8  kJ/mol.  was  obtained  by  second  law  method. 

Tile  value  of  entropy.  S'Wu  k for  a-lnaSet(s)  was  reported  by  [Dem  75,  Rym  75. 
Kos  84]  from  measurements  of  the  low  temperature  heat  capacity,  and  by  [Mus  74]  from 
EMF  measurements.  The  data  are  summarized  in  Table  5.7.  The  entropy,  S°»s,  for  y- 
InjSej(s)  was  estimated  from  S*m.ts  n of  ct-lnzSej(s)  using  the  heat  capacity  data  of  a- 
InrSejfs)  and  y-In;Sei(s)  and  AS",,..,  transition  entropy  data  retained  in  this  work.  The 
value,  S°2j«  i!K-  210.68  J/K  mol,  was  obtained  for  y-lniSej(s). 

The  transition  enthalpy  for  y-lniSej  to  8-lnrSej  reaction  was  estimated  from  the 
values  of  standard  enthalpy  of  formation  and  the  heat  capacity  of  y-InjSei  and  8-IniSei 
The  transition  entropy  was  then  calculated  from  the  equation,  ASn^AHn/Tv  The  value, 
S"i«i«k=260  99,  was  obtained  for  6-InjSej(s)  from  the  value  of  y-IniSej(s)  and  the 
AS",-,n  transition  entropy  data  retained  in  tltis  work. 


5.3.3  InsSer 

No  therniochemica]  data  are  available  for  the  InsSe?  phase.  For  ionic  or  quasi- 
ionic system,  the  shape  of  the  AH°r,misK  polygon  vs.  composition  usually  lies  outside  the 
basic  triangle,  and  the  summit  is  often  the  compound  with  the  highest  melting  point  [Kub 
79],  Following  this  trend,  the  standard  heat  of  formation,  AH°c2w.isk,  for  ln$Se7  is 
estimated  to  be  -784.4  kj/mole  from  Figure  5.1.  The  value  of  standard  entropy,  S'Wts.  is 
estimated  to  be  469.9  (J/K  mol)  by  Latimer's  rule  [Lat  75],  The  heat  capacity  is 
estimated  by  Unitl's  rule  [Kub  77].  The  enthalpy,  entropy  and  heat  capacity  data  arc  listed 
in  Tables  5.9, 5.10  and  5.11,  respectively. 

5.3.4  IbsSsm 

As  with  In-Set  there  are  no  ilierinochemical  data  available  for  InqSen  phase.  The 
standard  heat  of  formation,  AH°r,:9tu5K,  for  ln,>Sen  is  estimated  to  be  -1282  kj/mole  using 
Figure  5. 1 . The  standard  entropy,  S11;.;,.  is  estimated  to  be  796.8  (J/K  mol)  by  Latimer's 
title  [Lat  75],  The  heat  capacity  is  estimated  by  Unal's  rule  [Kub  77).  The  enthalpy, 
entropy  and  heat  capacity  data  are  listed  in  Tables  5.9,  5. 10  and  5.1 1,  respectively. 


Table  5.6  Standard  heal  of  formation  oflnjSej(s)  at  298.15  K 


Table  5.7  Reported  S%s.is  k values  for  ln;Se, 


Experimental  method 


Adiabatic  Calonmclrv 


Adiabatic  Calorimetry 


* calculate  in  this  work 


Table  5.8  Assessed  heat  capacity  for  ImSesfs) 


Cr=  a + bT  (J/K*mol) 


107  64907*0.054328  T 


1 1678158+0026741  T 


1 13.41683+0  035989  T 


S.3.5  iDsSfi+ts) 

The  standard  enthalpy  of  formation.  AH^stsK.  and  entropy.  S'Wisk.  for 
irlrSe?<s)  have  been  measured  by  Mustafaev  el  al.  [Mus  74],  but  the  values  were  assigned 
to  lnjSer,(s).  Chatillon  [Cha  93]  has  recalculated  the  AG“r(T)  = -761578+290.849  T 
(J/mol).  The  value,  AH'V.wrijk  =-688.63  (kJ/mol),  is  obtained  by  recalculation  of  [Mus 
74]'s  EMF  data.  The  value.  S'Wisk~514.  14  (J/mol).  is  obtained  from  the  second  law 
AS°r.t9s.tsK  value  from  [Mus  74].  The  value,  520.85  (J/mol),  is  obtained  if  estimated  by 
Latimer's  rule  [Lat  75].  The  heat  capacity  was  determined  by  Mills  [Mil  76]  from  DSC 


measurement  The  enthalpy,  entropy  and  heat  capacity  data  are  listed  in  Table  S.9,  5,10 
and  5,11,  respectively. 

5.3.6  InSets) 

The  standard  enthalpy  of  formation,  AH'Y;?*  isk.  was  determined  by  Hahn  and 
Burrow  [Han  52]  using  combustion  calorimetry,  and  by  Mustafaev  cl  al.  [Mus  74]  from 
their  EMF  measurements.  Values  of  the  enthalpy  and  entropy  of  fusion  were  reported  by 
Mills  [Mil  76]  using  DSC.  A value  for  the  absolute  entropy,  S°nw uk,  data  were  reported 
by  Mamedov  er  al.  [Mam  67]  from  low  temperature  heat  capacity  measurements.  The 
value,  S°29sisk  = 77.44  (J/K  mol),  is  obtained  from  the  [Mus  74]'s  second  law  AS"r.misK 
value  of  EMF  measurement.  Mills  [Mil  76]  measured  the  heat  capacity  using  DSC.  The 
enthalpy,  entropy  and  heat  capacity  data  are  listed  in  Tables  5,9,  5.10  and  5 11, 
respectively. 

5.3.7  iDaSsjOU 

The  standard  enthalpy  of  formation,  AH'V.2«i>k,  and  entropy,  AS"kmi.i!k,  for 
ln2Se(s)  have  been  measured  by  Mustafaev  el  al.  [Mus  74],  which  should  be  assigned  to 
ImSeifs)  based  on  the  correct  phase  diagram.  The  values,  AH’Vm  ISk  =-332.24  (kj/mol) 
and  S°2ssuk=336.83  (J/K  mol),  were  oblained  by  recalculation  of  [Mus  74]’s  EMF  data. 
Chatillon  [Cha  93]  estimated  the  entropy  value  to  be  292.9  (J/K  mol)  using  Latimer’s  rule 
[Lat  75]  and  the  heat  capacity  data  by  Ungl,s  rule  [Kub  77],  The  enthalpy,  entropy  and 
heat  capacity  data  were  listed  in  Table  5.9,  5. 10 and  5.1 1 , respectively. 


Table  5.9  Standard  enthalpy  of  formation  values  at  298,15  K,  AH'Y.IWIJK.  of  In-Se 


Species 

Experimental  method 

Reference 

In,  Sc? 

Estimated 

This  work 

793.75 

This  work 

1215.9 

Estimated 

This  work 

1337.6 

Assessed 

This  work 

580.6*  16.3 

EMF 

[Mus  74] 

68863 

EMF 

[Mus  741* 

808  6 ± 0.5 

EMF 

[Mai  99) 

855  67 

Assessed 

This  work 

InSe 

1 17.0  ± 12.5 

Combustion 

Calorimetry 

[Han  52] 

109  1 ± 7 9 

EMF 

[Mus  74) 

1 18.0  ± 12.6 

Evaluated 

[Mil  74) 

98 

Evaluated 

[Cha  93] 

[Mai  99] 

Assessed 

This  work 

In-Se 

147  1 * 1 1 7 

EMF 

[Mus  74] 

332.24 

EMF 

[Mus  741* 

300.7 

Evaluated 

[Cha  931 

This  work 

Note:  * recalculated  in  this  work,  **  should  be  InjSe, 


Table  5.10  Entropy  values  al  298.15  K ofln-Se  compounds. 
ln?Se7 


lixpctimcnlal  method 


Adjabauc  Calorimetry 


e:  "recalculated  in  this  work,  **  should  be  IttaSej 


Table  5.11  Heat  capacity  data  for  In-Se  compounds. 


' Estimated  by  U rial's  rule.  [Kub  77] 


>us  molecules  In(g),  ImSei(g),  In.Selg),  and 
lnSc(g)  were  studied  by  [Col  68,  Gre  71,  Gre  73]  and  assessed  by  [Cha  93],  including 
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5.4  Gibbs  Energy  Data 


5.4.1  EMF  measurements 

Mustafacv  el  al.  [Mus  74]  measured  the  indium  activity  in  the  temperature  range 
510  to  710  K by  electromotive  force  in  the  two  phase  domains,  InjSer-Se,  (ImScr,)- 
ln-Se.t,  InSe-On-Se,.).  (ln,Se)-lnSe  In, Sc..  should  be  corrected  to  InsSer  and  ln:Se  should 
be  corrected  to  InjSej.  More  recently  Mallika  el  al.  [Mai  99]  measured  the  indium 
activity  in  solid  and  liquid  (indium  4 selenium)  solutions  to  SO  mass  per  cent  selenium  in 
the  solid-state  galvanic  cell. 


The  following  gas  species  were  reported  to  exist  in  In-Se  system,  Se„(n=l  to  8) 
(g),  ln(g),  In2Se2<g),  InjSe(g),  and  InSe(g).  Both  total  vapor  pressure  [Gri  72,  Ven  79,Gri 
71,  Gre  73,  Sri  85,  Sri  87]  and  partial  pressure  [Ber  69,  Col  68]  measurements  in  the  In- 
Se  system  are  presented  in  Table  5. 14. 

Colin  and  Drowart  [Col  68]  determined  the  partial  pressure  in  equilibrium  with 
liquid  InSe  by  Knudsen  cell  mass-spectrometry.  In  their  investigation,  the  condensed 
phase  composition  was  not  constant  due  to  fact  that  the  ion  intensity  ratio 
l(Se2~)/l(IniSe')  varied  appreciably  with  temperature.  Grinberg  el  al.  [Gri  71]  studied  the 
saturated  total  vapor  pressure  above  liquid  InSe  by  Bourdon  quartz  gauge  and  the  vapor 
constituents  by  electron  and  IR  absorption  spectroscopy. 


Sublimation  of  ln:Set(s)  has  been  studied  by  mass-spectrometriy  [Ber  69],  torsion 
or  Knudsen  effusion  [Ven  79.  Pas  88,  Bar  88,  Sri  87.  Sri  85],  and  membrane  manometry 
[Gri  72],  ln:Se,(s)  was  reported  to  sublime  congruently  according  to  the  reaction  [Ber  69, 
Pas  88.  Bar  88]: 

InjSej  (s)  = toiSe  (g)  + Sej  (g) 

Congruent  behavior  has  not  always  been  observed  and  the  occurrence  depended 
on  the  experimental  conditions.  One  reason  can  ascribe  to  the  fact  that  the  effusion 
process  may  change  the  condensed  phase  composition.  This  efffecl  was  shown  in 
Srinivasa  and  Edwards  investigations  (Sri  87.  Sri  88].  The  difference  between  congruent 
sublimation  and  congruent  effusion  was  discussed  in  (Edw  9]].  In  their  studies,  the 
condensed  phase  was  not  the  Iti;Sei(s)  one  phase  domain,  but  a liquid-ln:Se,  two-phase 
boundary.  Their  data  agreed  qualitatively  with  the  P-T  phase  diagram  from  Grinberg  el 
al.  [Gri  72],  although  the  pressure  measured  by  [Sri  85,  Sri  87]  is  much  lower.  Grinberg 
el  at  [Gri  71]  studied  the  vapor  phase  along  the  three-phase  boundary  “ ln’Se_i(s)-liquid- 
vapor ",  by  measuring  the  total  pressure  and  optical  absorption  spectra  of  the  vapor.  They 
constructed  the  P-T  projection  of  the  phase  diagram  along  the  Iil;Se\(.s)  liquidus.  The 
spectroscopic  analysis  indicated  three  major  components  in  the  saturated  vapor,  InjSe  (g). 
Sc-  (g),  and  In  (g).  A small  homogeneity  range  (-0.1  at.%)  was  found  in  InjSej  (s)  phase 
by  analyzing  the  vapor  composition  along  the  three  phase  boundary. 


I OS 


5 5 Thermodynamic 


ns  and  stoichiometric  compound  phases 


= y,.°G'„  +ys.’G's,  (5.2) 
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The  ideal  mixing  Gibbs  energy  "!G'  is  equal  lo 

"'G'  = RT(y,„  In  yh  + In  _phjS0  +j'x,lajrs,)  (5.3) 

The  excess  Gibbs  energy  tGl  is 

® = JVj  + /-ta.hjSVj  CP/a  ~ J'/r.jS.j  ))  (54) 

+ Tto2s^Vs.(%.!s», ,Sr+'4/^Sfl.S.CV/«jS<j  ■ J’s.)) 
where  the  four  terms  L represent  the  interactions  between  the  species.  They  can  be 
expressed  as  the  function  of  temperature. 


5,6  Optimization  Procedure 

A selected  set  of  thermodynamic  and  phase  diagram  data  were  used  for  the 
optimization  of  thermodynamic  model  parameters  of  all  phases  in  this  system.  The 
optimization  was  performed  by  using  the  PARROT  module  of  the  Thermo-Calc  program 
package  [Sun  85]  and  procedures  in  four  steps. 

First,  initial  estimates  of  the  Gibbs  energy  coefficients  of  the  stoichiometric 
compounds  InjSe,,  InSe,  Irv.Se?.  lri.;Sen.  InjSe?,  Ji-in.Sei.  v-ln-Se*.  5-IniSer,  and  o- 
InjSej,  were  deduced  from  the  reported  enthalpy  of  formation,  standard  entropy,  heat 
capacity.  No  enthalpy  of  melting  for  InjSej  was  reported  in  the  literature  The  Gibbs 
energy  coefficients  for  liquid  InjSej  were  deduced  from  the  melting  enthalpy  of  InSe 
reported  by  [Mil  76],  Second,  a rough  optimization  of  the  liquid  coefficients  was 
performed  based  on  the  activities  of  selenium,  miscibility  gap,  and  three-phase  equilibria 
(those  involved  the  liquid  phase).  Third,  fixing  the  coefficients  of  the  liquid  the 
coefficients  of  the  stoichiometric  IniSe.i,  InSe,  InriSe?.  In,, Sen.  InsStr?,  P-ln^Sei.  y-IniSej. 
5-ln;Se„  and  a-ln2Sej  phases  were  optimized  based  on  the  enthalpies  of  formation. 


5.7  Result  and  Pis 


The  optimized  parameters  of  the  stable  phases  in  the  In-Se  system  are  listed  in 
Table  5.15.  The  phase  diagram  and  thermodynamic  properties  of  this  system  were 
calculated  with  the  Poly-3  module  of  the  Thermo-Calc  program  package.  The  calculated 
phase  diagram  is  shown  in  Figure  5.2.  Comparisons  between  the  experimental  and 
calculated  phase  diagram  are  shown  in  Figures  5.3  to  5.6.  Table  5. 1 presents  the 
experimental  and  calculated  temperatures  and  compositions  of  the  invariant  reactions  in 
this  system. 

The  calculated  values  are  well  within  the  uncertainty  of  experimental  data.  The 
calculated  In-rich  liquid  miscibility  gap  agrees  well  with  the  experimental  data  of  [Gla 
91]  and  [God  98]  as  shown  in  Figure  5.4.  The  calculated  peritetic  reactions  and  liquidus 
lines  for  ltijSe;,  InSe,  InsSe?,  InsSe.i,  In.Se?  and  (l-ln?Se.i  are  shown  in  Figure  5.5  in 
detail  along  with  the  experimental  data.  In  general,  the  calculated  diagram  is  in 
agreement  with  the  experimental  data  The  calculated  phase  diagram  at  the  Se-rich  side  is 
given  in  Figure  5.4  in  comparison  with  the  experimental  data  Figure  5.5  shows  the 
comparison  between  the  calculated  and  measured  activity  of  indium.  The  calculated 
activity  deviated  more  positively  than  the  experimental  data. 

Figure  5.7  shows  the  comparison  between  the  calculated  and  measured  total  vapor 
pressures  in  this  system.  To  calculate  the  total  pressure,  an  ideal  gas  phase  containing  12 
species  is  added  to  the  optimized  coefficient  set.  The  parameters  of  the  gas  species  Se, 
Se2,  Sej,  Se4,  Sej,  Se6,  Se?.  and  Seiare  from  chapter  3,  whereas  the  parameters  of  the  gas 
species  In,  InSe,  ln:Se.  and  In:Se;  are  from  [Sun  94].  The  calculated  total  pressures  agree 
well  with  the  measured  data.  The  calculated  total  pressure  data  are  in  good  agreement 


wilh  the  measurements  of  [Gri  72).  The  calculated  T-p  data  along  the  0-In-Se,  liquidus 
agree  well  with  the  measurements  of  [Sri  87]  and  [Pas  89). 

The  experimental  and  calculated  standard  enthalpies  of  formation  of  the 
intermediate  compounds  at  298.15  K in  the  ln-Sc  system  are  presented  in  Tables  5.6  and 
5.9.  The  corresponding  standard  entropies  are  presented  in  Tables  5.7  and  5.10.  The 
reference  states  of  these  data  are  the  enthalpies  of  the  pure  stable  elements  at  298. 1 5 K. 
The  assessed  enthalpy  of  formation  for  a-ln:Sej  in  the  present  work  are  close  to  the 
synthesis  calorimetry  value  [Sha  67).  Their  value  was  obtained  by  direct  synthesis  of  a- 
ImSei  in  a bomb  calorimeter.  The  synthesized  compound  was  confirmed  by  X-ray 
difiraction.  This  technique  does  not  rely  on  other  thermodynamic  functions,  which 
introduces  additional  uncertainty.  The  assessed  standard  enthalpies  for  p-In;Se»  and  y- 
ln;Sei  are  close  to  the  values  derived  from  [Sri  87]  and  [Sha  67],  respectively.  The 
assessed  enthalpy  of  formation  for  InSe  is  well  within  the  experimental  uncertainty  of  the 
value  measured  by  [Han  52]  ffom  combustion  calorimetry.  It  was  found  that  the  enthalpy 
of  formation  value  for  InsSe?  is  much  more  negative  than  the  value  derived  ffom  EMF 
me  e ne  it  of  [Mus  74]  and  closer  to  the  value  by  [Mai  99].  It  was  also  found  that  the 
enthalpy  of  formation  value  for  ImSer  is  much  more  negative  Ilian  the  value  evaluated  by 
[Cha  93]  and  the  recalculated  value  ffom  [Mus  74).  A more  positive  value  of  the 
enthalpy  of  formation  resulted  in  the  disappearance  of  lnjSej  phase  in  the  calculated  T-x 
phase  diagram. 

The  assessed  entropy  values  are  well  within  the  available  experimental  values.  The 
calculated  heat  capacity  data  for  a-ln:.Se(,  y-ImSes,  |i-In,Sei.  ln,.Se?  and  InSe  along  with 


the  measured  dala  by  [Mil  76]  are  given  in  Figure  5.8,  5,9,  and  5,10,  respectively.  The 
calculation  values  agree  well  with  the  measured  ones. 


5,8  Conclusions 

A thermodynamic  description  of  the  In-Se  binary  system  is  obtained  by 
optimization  of  the  available  phase  equilibrium  and  thermodynamic  data  The  associate 
solution  model  was  used  to  represent  the  Gibbs  energy  of  the  liquid  phase.  The 
intermediate  phases  were  modeled  as  line  compounds.  The  phase  diagram  and 
thermodynamic  properties  of  this  system  have  been  calculated  by  using  the  optimized 
model  parameters.  There  is  reasonable  agreement  between  the  model-calculated  values 
and  the  phase  equilibrium  and  thermodynamic  data  available  in  the  literature  for  all 


Figure  5 2 Calculated  In-Se  phase  diagram  based  on  the  optimized  parameters. 


: 5.3  Calculated  In-Se  phase  diagram  along  with  the  experimental  i 


Figure  5.4  Comparison  between  the  calculated  In-Se  phase  diagram  and 
experimental  data  in  the  vicinity  of  the  In-rich  miscibility  gap. 


Figure  5.5  Comparison  between  the  calculated  In-Se  phase  diagram  and  ■ 
experimental  data  in  the  vicinity  of  Se-rich  miscibility  gap. 
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Figure  5.6  Comparison  between  the  calculated  In-Se  phase  diagram  and  • 
experimental  data  from  0.3  to  0.65  mole  fraction  of  Se. 


Figure  5.7  Calculated  activity  of  In  with  the  experimental  data  [Mai  99]  with  the  pure 
liquid  In  phase  as  the  reference  state. 


Figure  5.8  Comparison  of  the  lotal  vapor  pressure  in  the  In-Se  system  with  the 
experimental  data. 
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Figure  5.9  Assessed  and  measured  heat  capacity  data  for  a-ImSe>,  p-ln.Se,.  and  y-InjSej. 
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TEMPERATURE_KELVIN 


Figure  5. 1 1 Assessed  and  measured  heat  capacity  data  for  InSe. 


CHAPTER  6 

CulnSe;  THIN  FILM  FORMATION  BY  RAPID  THERMAL  PROCESSING 


6.1  Introduction 

In  conventional  absorber  film  formation  processes,  CIS  is  not  directly  formed 
during  the  thin  film  deposition  process.  Rather,  the  requirement  for  a large  area,  high 
throughput,  and  low  cost  absorber  film  formation  process  has  motivated  the  use  of 
physical  deposition  techniques  using  elemental  sources  to  first  deposit  a precursor  film 
of  an  appropriate  overall  composition.  Subsequent  annealing  then  leads  to  the 
synthesis  of  large  grained  CIS.  Selenization  and  Stacked  Elemental  Layer  (SEL) 
techniques  are  two  promising  processes.  Selenization  is  a two-stage  technique.  Firstly, 
sequential  layers  of  copper  and  indium  are  deposited  on  a molybdenum-coated 
substrate.  CIS  is  then  obtained  by  reacting  the  metal  precursor  with  gas  phase  Se. 
either  H,Se  or  Se  vapor.  The  SEL  technique  is  a process  involving  the  sequential 
deposition  of  Cu,  In  and  Sc  layers  on  a substrate  which  is  then  annealed  in  vacuum  or 
a Se  environment  to  promote  reaction  of  the  elements  to  form  CIS.  Among  the 
annealing  techniques,  RTP  is  attractive  with  the  advantages  of  low  thermal  budget  and 


better  control  of  the  process  kinetics. 


6.2  Previous  Work  on  RTP.  CIS 


Several  precursor  structures  have  been  utilized  in  RTF  synthesis  of  CulnSe:  thin 
films  Basically,  three  types  of  precursor  structures  have  been  reported  in  the  literature 
One  type  of  precursor  structure  stacked  the  elemental  components  to  form 
sandwich-like  structure.  A variety  of  precursors  could  be  fabricated  by  alternating  the 
layer  sequence  and  changing  the  number  of  layers.  Compositions  of  the  films  were 
controlled  by  the  relative  thickness  of  the  individual  layers.  It  often  incorporated  excess  Se 
due  to  the  high  volatile  nature  of  Sc  gas  species.  A number  of  groups  have  reported 
successful  CuInSej  thin  film  formation  using  this  approach,  including  [Oum  89,  Kar  93, 
Bau  94],  or  Cu(ln,.,Ga,)Se;  [Pro  94,  Pro  96]. 

The  important  processing  parameters  in  this  approach  include  stacking 
arrangement  of  the  layers,  Cu/In  ratio,  RTP  ramping  rate,  and  annealing  temperature. 

The  first  device  quality  CIS  produced  by  RTP  was  synthesized  by  Karg  e/ «/.  [Kar 
93],  It  consists  of  Mo,  Cu,  In , and  Se  sequentially  deposited  by  DC-magnetron  sputtering 
(Mo,  Cu.  In)  and/or  evaporation  (Cu.  In.  Se).  The  Cu/In  ratio  varied  from  0.85  to  0.95, 
while  Se  was  typically  4054  in  excess  of  the  stoichiometric  amount  A furnace  type  (low 
ramp  rate,  f'C/s)  processing  was  also  performed  in  their  investigation.  Two  critical 
process  steps  were  found  from  this  type  of  process.  The  first  one  is  the  severe  dewetting 
of  the  amorphous  Se  layer  starting  at  about  I00°C  which  cause  a laterally  inhomogeneous 
nuolcation  and  grain  growth  of  the  CIS  layer.  This  laterally  inhomogeneous  growth 
process  resulted  in  voids  between  the  CIS  and  Mo  back  contact,  which  produced  poor 
adhesion  and  low  cell  performance.  The  second  critical  process  period  is  when  the  sample 
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in  the  temperature  range  of  100  to  250”C  where  the  formation  of  CuiSe  is  favored  as 
compared  to  CulnSej,  a detrimental  secondary  phase.  Their  results  indicated  the  Cu/ln 
ratio  is  the  most  significant  parameter.  At  a Cu/ln  precursor  of  0.89  small  grained  material 
is  obtained  in  contrast  to  the  coarse  grained  film  at  Cu/ln  of  0.92  with  a grain  size  in  the 
range  of  I '3  pm.  TEM  results  revealed  that  copper-rich  and  indium-rich  films  were 
distinct  in  their  defect  structure.  Copper-rich  films  had  a significant  number  of  dislocations 
and  stacking  faults.  In  contrast  to  copper-rich  films,  there  were  a significant  number  of 
twins  in  In-rich  films.  These  features  revealed  different  reaction  mechanisms  during 
synthesis  for  these  two  Cu/ln  ratios. 

Using  the  same  process  strategy,  [Pro  94]  fabricated  the  Cu(ln,Ga)SeJ  thin  film 
cells  by  RTP.  They  also  developed  two  ways  to  incorporate  sodium  into  CIS  to  quantify 
its  impact  on  grain  morphology  and  cell  performance.  They  varied  the  outdiffusion  of 
sodium  from  the  soda  lime  glass  by  employing  a partially  permeable  Mo-backclectrode  in 
the  first  process.  Sodium  selenide  was  added  to  the  precursor  in  the  second  process.  The 
optimization  of  the  sodium  amount  in  RTP-CIGS  films  on  soda  lime  glass  resulted  in  cells 
with  averaging  efficiency  12%  and  a peak  efficiency  of  13.2%. 

More  recently  Probst  el  at.  [Pro  96]  developed  a new  stacked  elemental  layer 
process  to  produce  Cu(In.Ga)Sc,  and  Cu(In,  GaXSe,  S)a  cell.  The  key  process  steps  are 
I (deposit  a barrier  coating  on  the  soda  lime  glass  substrate  combined  and  add  a sodium 


thin  layer  ofSiO.,  AljOj  or  Si,N,  was  used  as  the  barrier  material.  They  suggested  a liquid 
assisted  growth  mechanism  for  the  stacked  elemental  layer  process  has  two  main 
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prerequisites:  I)  the  processing  conditions  must  favor  the  formation  of  the  CuSe  phase,  2) 
a dynamic  heating  process  to  exceed  the  melting  temperature  of  CuSe  before  it  is 
consumed  in  the  reaction.  Cell  peak  efficiencies  of  13.3%  for  Cu(In,Ga)Se;  and  of  14.6% 
for  Cu(ln,Ga)(ScS>2  devices  were  obtained  with  this  technique 

The  second  precursor  structure  used  co-evaporating  elemental  Cu,  In  and  Se  [Moo 
91,  Gab  92],  This  processing  approach  had  bubbling  and  adhesion  problems,  especially 
with  Mo-coated  glass  substrates.  The  bubbling  was  worse  with  higher  Cu  contents  and 
faster  ramp  rates.  The  device  made  from  this  approach  had  an  efficiency  < 1%.  after  2hrs 
200"C  air  annealing  an  efficiency  3,5%  was  obtained. 

The  third  approach  used  RTP  to  re-crystallized  CulnSe:.  Albin  el  al.  [Alb  91] 
rapid  thermal  annealed  small-grained  (SI  pm)  CulnSe-  to  producing  large-grained 
CulnSe2  (-100  pm).  This  approach  suffered  from  the  severe  oxygen  incorporation  and 
CulnSe2  decomposition  at  high  temperature,  and  other  defects  such  as  excessive  voids  and 
pinholes  in  the  Dims. 

It  is  learned  from  this  review  that  the  precursor  film  structure  and  the 
subsequent  reaction  pathway  during  annealing  can  dramatically  affect  the  eventual 
device  performance.  For  example,  studies  by  [Mat  95]  showed  a different  chemical 
reaction  pathway  when  the  precursor  was  formed  by  mixing  powders  of  pure  Cu,  In 
and  Sc  in  the  1:1:2  stoichiometry  as  compared  to  mixing  powders  of  the  intermetallic 
compound  Culn  with  Se  in  the  same  stoichiometry.  The  latter  approach  led  to  CIS 
formation  at  a much  lower  temperature. 

This  is  consistent  with  the  observation  that  a co-deposited  Cu+ln+Se 
amorphous  mixture  [Gab  92,  Moo  91]  precursor  structure  does  not  lead  to  device 
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quality  CIS  films,  in  contrast  to  the  successful  SEL  process.  In  the  SEL  process  the 


elements  are  sequentially  deposited  on  a substrate,  usually  in  the  sequence 


substrate/Cu/In/Se  [Pro  96],  This  precursor  structure  usually  leads  to  the  formation  of 
a Cu,ln,  intcrmetallic  phase  during  an  intermediate  reaction  step  [Kar  93).  Based  on 
these  results,  it  is  anticipated  that  alternative  precursor  structures  might  yield  improved 
CIS  thin  film  growth. 


Phase  diagram  information  is  essential  in  designing  the  precursor  structures 
Figure  6. 1 shows  the  Cu-In-Se  ternary  composition  diagram  along  with  the  accepted 
compounds  and  several  tie-lines  connecting  pairs  of  them.  Inspection  of  this  diagram 
suggests  several  reaction  pathways  to  synthesize  CulnSe2  from  binary  compounds 
using  RTP.  Although  previous  approaches  have  started  with  the  pure  elements  to  form 
a physical  mixture  or  a slacked  layer  structure,  certain  advantages  exist  when  starting 
with  binary  compounds  Among  these  compounds  those  within  the  Cu-Se  and  In-Se 
couples  are  attractive  due  to  existence  of  low  melting  temperature  phases  [Sta  97). 
There  are  several  combinations  of  these  two  couples  that  could  lead  to  CIS  formation. 
The  following  guidelines  were  used  to  suggest  suitable  Cu-Sc  and  In-Se  compounds  to 
serve  as  precursors  for  ex-situ  RTP  synthesis  of  CIS. 

• Avoid  high  melting  temperature  phases  such  as  In-Sej  and  Cu-Se, 
especially  the  conducting  Cu:Se  secondary  phase. 

• Utilize  an  intermediate  liquid  phase  to  assist  grain  growth. 

• Minimize  evaporation  losses  during  RTP 


Select  highly  reactive  pairs. 


At  first  glance,  one  choice  would  be  to  deposit  a bilayer  structure  of  CujSe  and 
ln:Se-.  and  react  them  to  form  CIS  according  to: 

Cu;Se*TniSej-*2CuInSe; 

Examining  the  evaluated  Cu^Sc-ln^Se.  pseudobinary  T-x  pltase  diagram  (Figure  6.2) 
reveals  that  the  liquid  phase  that  exists  at  the  lowest  temperature  in  this  section  is  the 
-870"C  eutetic  alloy  lying  in  the  In;Se;-rich  composition  field.  This  temperature  is 
much  too  high  for  the  glass  substrate. 

An  examination  of  Cu-Se  (Figure  6.3)  and  In-Se  (Figure  6.4)  binary  T-x 
diagrams  from  chapters  4 and  5 suggests  that  a more  suitable  combination  would  be 
selenium-rich  copper  compounds  (CuSe,  CuSc;)  and  selenium-rich  indium  compounds 
(InSe.IntSei)  When  these  binary  compounds  are  heated  the  copper  compounds  will 
form  a liquid  phase  and  the  indium  compounds  will  form  a solid-liquid  mixture  in  the 
temperature  range  523  to  620*0.  To  reach  the  liquid  phase,  rapid  thermal  processing 
might  be  required  to  avoid  solid-solid  phase  transformations.  The  binary  diagrams 
alone,  however,  do  not  indicate  if  a liquid  phase  exists  after  the  reacting  precursor  film 
enters  the  ternary  domain. 
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Figure  6. 1 Cu-ln-Se  ternary  composition  diagram. 
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Figure  6.2  CurSe-IihSej  pseudobinary  phase  diagram. 


Figure  6.3  Cu-Se  binary  T-x  phase  diagram. 
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Figure  6.4  In-Se  binary  T-x  phase  diagram. 


Se  should  easily  volatilize  to  leave  single  phase  CIS.  An  estimation  of  the 
ternary  Cu-ln-Se  phase  diagram  has  been  developed  and  a projection  of  the  liquidus 
surface  is  shown  in  Figure  6.6.  With  this  diagram  possible  compound  combinations  can 
be  forther  evaluated.  In  particular,  the  eutetic  valley  between  the  binary  eutetic 
compositions  Ei  and  E,  should  be  nearly  degenerate  and  thus  exist  at  slightly  below 
220°C.  The  dash  line  in  Figure  6.6  connecting  the  pure  Se  comer  to  the  CulnSe- 
compound  indicates  it  should  be  possible  to  form  CIS  in  the  presence  of  an 
intermediate  low-temperature  liquid  phase  by  the  reaction: 

Liquid  -*  CulnSe..  + Se 

The  relative  amounts  of  the  binary  compounds  should  be  adjusted  to  lie  on  this  line 
The  above  approach  suggested  a precursor  structure  that  consists  of  an  In-Se 
layer  (50  to  55  Se  at.%)  (see  dashed  box  in  In-Se  phase  diagram  in  Figure  6.4) 
deposited  on  a Cu-Se  (50  to  55  Se  al.%)  (see  dashed  box  in  Cu-Se  phase  diagram  in 
Figure  6 3)  layer,  which  is  deposited  on  a Mo  coated  soda-lime  glass  substrate  A 
schematic  diagram  of  this  structure  is  given  in  Figure  6.5.  The  composition  ranges  for 
the  binary  pairs  and  the  overall  precursor  films  are  indicated  in  Figure  6.6  as  heavy 
lines  along  the  In-Se  and  Cu-Se  limits  near  the  CuSe  and  InSe  compounds.  With  this 
composition  a Cu-Se  liquid  phase  should  form  to  assist  growth  and  the  small  amount 
of  remaining  excess  Se  should  be  volatile.  The  presence  of  solid  In-Se  layer  on  top  will 
minimize  the  evaporation  loss  during  the  RTP  process. 


Figure  6,5  Proposed  precursor  f 


Figure  6.6  Liquidus  projection  for  Cu-In-Se  ternary  system. 


136 

6.4  Precursor  Growth 

A low  lemperaiure  precursor  growth  is  desirable  for  higher  manufacturing 
equipment  throughput,  lower  equipment  cost,  lower  thermal  budget,  and  thus  a lower 
capital  cost  for  modules.  A rotating-disc  migration-enhanced  deposition  system  was 
used  for  precursor  deposition,  which  included  two  Langmuir  evaporation  sources  for 
copper  and  indium  and  a double-oven  thermal  cracking  source  for  selenium.  A 
schematic  drawing  of  the  reactor  is  shown  in  Figure  6.7.  The  growth  conditions 
(substrate  temperature  and  (lux  ratio)  of  In-Se  binary  system  have  been  thoroughly 
studied  by  [Eme  92],  In  their  investigations  they  utilized  a Langmuir  evaporation 
source  and  suggested  the  low  selenium  sticking  coefficient  is  most  likely  due  to  the 
formation  of  larger  polymer  molecules  (e.g.,  Sea).  The  use  of  a cracking  source 
allowed  the  Se  decomposition  temperature  to  be  adjusted  independent  of  source 
temperature  (flux  control)  and  thus  control  the  selenium  flux  species  distribution. 

The  precursor  films  used  in  this  study  were  deposited  at  a 200°C  substrate 
temperature.  To  fabricate  the  desired  precursor  structure,  the  In-Se  and  Cu-Se  growth 
conditions  were  studied  separately  by  fixing  the  selenium  flux  constant  and  adjusting 
the  selenium  to  copper  and  selenium  to  indium  molar  flux  ratios.  The  molar  fluxes 


t cross-calibrated  by  EIES  and  quartz-crystal  i 


137 


Figure  6.7  Sc 


ition-enhanced  < 
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6.5  Rapid  Thermal  Processing 

Ex-situ  RTP  was  performed  on  samples  placed  in  a closed  rectangular  quartz 
susceptor  A flat  sheet  of  SiC  coated  graphite  was  placed  at  the  bottom  of  the 
susceptor  to  improve  the  temperature  uniformity.  The  susceptor  and  sample  were 
loaded  into  a customized  RTP  furnace  heated  by  a bank  of  12  tungsten-halogen  quartz 
lamps  from  both  sides  under  flowing  NV  The  schematic  diagram  of  the  RTP  reactor  is 
given  in  Figure  6.9.  The  processing  temperature  was  controlled  by  a thermocouple  in 
direct  contact  with  the  glass  substrate.  The  RTP  annealing  recipes  was  input  and 
controlled  by  a PID  controller.  Due  to  the  significant  difference  in  the  absorption 
characteristics  between  the  precursor  films  and  the  soda-lime  glass  substrate  It  is 
anticipated  a non-negligible  temperature  gradient  between  the  thermocouple  reading 
and  surface  film  temperature.  This  temperature  difference  was  estimated  by  visualizing 
the  eutetic  reaction  at  575°C  of  a Mg/Au  layer  deposited  prior  to  the  ln-Se  compound 
(see  Figure  6.8).  The  In-Se  layer  is  about  1000  A.  The  temperature  difference  is 
estimated  to  be  around  200"C. 


Figure  6.  8 Temperature  calibration  sample. 
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rf  the  rapid 


6.6  Film  Characterization 


Electron  Microscope  (SEM).  Compositional  analysis  was  performed  by  Inductively 
Coupled  Plasma  (ICP)  emission  spectroinetty  for  the  atomic  ratios  of  CIS.  The  samples 
were  dissolved  in  5 % (by  volume)  nitric  acid.  The  standard  solution  is  also  based  on  S% 
nitric  acid.  Composition  depth  profiles  were  measured  by  sputter  Auger  Electron 
Spectroscopy  (AES)  using  a Perkin-Elmer  PH1660  AES/SAM  microscope,  equipped  with 
a CMA  spectrometer  and  an  Ar-ion  gun.  All  spectra  were  obtained  at  SKeV. 

The  structural  properties  of  the  films,  such  as  phase  constitution,  crystallinity,  and 
texture  were  determined  by  Powder  X-ray  Diffraction  (PXRD)  and  Grazing  Incidence  X- 
ray  Diffraction  (GIXRD)  at  various  angles  The  GIXRD  measurements  were  carried  out 
using  a Siemens  diffractometer  equipped  with  an  area  detector  using  Cu  Ka  radiation. 
This  non-destructive  technique  allows  determining  the  depth  distribution  of  the  phases  in 
the  films.  The  utilization  of  an  area  detector  provided  the  opportunity  to  assess  the 
texturing  of  the  films  at  the  same  time.  This  feature  is  essential  when  diffraction  peak 
intensities  are  needed  to  identify  the  phase  constitution. 

The  unpolarized  Raman  scattering  signals  were  measured  using  a double 
monochromator,  an  Ar-  laser  excitation  at  two  wavelength  (A,=514,5nm  and  477  nm).  and 
standard  backscattering  sample  geometry.  The  longer  wavelength  excitation  laser  will 
provide  information  with  deeper  penetration  in  the  film. 


6 7 Results  and  Discussion-Precursor  Structure-! 


Two  precursor  smiclures  were  fabricated.  The  first  one  possesses  a three  layered 
structure,  ln1Se.i/CuInSc,/Cu;Se.  The  second  one  possesses  a two  layered  structure, 
InSe/CuSe 

6.7.1  Precursor  structure 

It's  a challenge  to  distinguish  CulnSe?  from  Cu:Se  using  XRD  spectra  since  most 
of  the  peaks  overlap.  A major  difference  in  the  powder  spectra  between  these  two  phases 
is  the  peak  intensity  ratio  of  the  (1 1 1)/(022)  = 0.5  for  CuiSe  and  (II 2)/(024)  = 2 for 
CuinSc;  phase.  Simulated  powder  XRD  patterns  for  CulnSe.,  and  Cu^Se  were  given  in 
Figure  6.10(a)  and  (b),  respectively.  Unfortunately,  the  possibility  of  preferred  orientation 
renders  this  difference  useless  in  a conventional  diffractometer  equipped  with  a point 
detector.  With  the  aid  of  an  area  detector,  the  difference  can  be  exploited.  An  area 
detector  allows  assessment  of  the  effects  of  preferred  orientation.  It  is  clearly  seen  in 
Figure  6. 1 1 that  only  those  diffraction  rings  belonging  to  the  Mo  substrate  indicate  signs 
of  texturing.  The  XRD  spectra  of  the  precursor  film  taken  at  various  X-ray  incidence 
angles  (0.5  to  6")  clearly  shown  a change  of  the  relative  intensity  of  the  peak  around  20  = 
26.5"  [(filers. . (I I2)ciiwu]  and  the  peak  around  20  = 44.5"  [(022)c„k„  (024)0*5*]. 
The  relative  intensity  for  peak  at  44.5"  is  increasing  with  the  X-ray  incident  angle,  as 
shown  in  Figure  6.12  (a)  and  (b)  All  the  spectra  show  the  presence  of  IniSe>  phase  in 
addition  to  CurSe  and  CulnSe:  These  data  suggested  a three  layered  structure, 
Mo/Cu;Sc''CLilnSe;/ln.,Se,1  was  formed  in  the  deposition  process. 
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Figure  6 10  Simulated  powder  Cu-Ka  XRD  patterns  for  (a)  CulnSe?.  and  (b) 
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Figure  6.11.  X-ray  diffraction  pattern  at  6"  incident  angle  of  the  precursor  film. 
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Figure  6. 12  XRD  spectra  of  precursor  trim  measured  ar  (a)  0.5"  incident  angle  and  (b)  6" 
incident  angle. 
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The  Raman  speclra  shown  Figure  6. 1 3(a)  and  (b)  provide  further  evidence  for  ihis 
structure.  A much  more  intense  CulnSe;  phonon  peak  is  evident  in  the  spectrum 
measured  at  514  nm  as  compared  to  the  one  measured  at  477  nm. 


Figure  6.13  Raman  spectra  at  two  excitation  wavelengths  (a)  \ = 477  nm  and  (b)  - 514 
nm  measured  on  the  precursor  film. 

RTP  290 

The  majority  of  the  ItuSe?  phase  disappeared  and  reacted  to  form  CuinSc;  after 
290“C  (TC  temperature)  RTP  annealing  for  70s  as  indicated  by  the  PXRD  spectrum 
shown  in  Figure  6.14.  The  appearance  of  the  (101)  peak  and  the  increasing  (1 12)  peak 
intensity  of  CulnSej  are  good  indications  of  productions  of  additional  CIS.  The  G1XRD 
patterns  vary  only  slightly  with  various  X-ray  incident  angles.  This  indicated  the  binary 


precursor  layers  were  reacting  and  formed  more  CulnSe;  during  the  RTP.  GIXRD  spectra 
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, however,  still  shows  the  presence  of  liuSer  peaks.  This  indicated  the  reaction  is  not 
complete  at  this  processing  condition.  The  XRD  spectra  at  6"  incident  angle  are  given  in 
Figure  6.  IS. 


Figure  6. 14  PXRD  spectrum  of  the  film  after  RTP  at  290°C  (TC  temperature)  for  70 : 
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Figure  6. 1 5 (a)  X-ray  diffraction  pattern  (b)  integrated  spectra  at  6°  incident  angle  of 
the  film  after  RTP  at  290”C  (TC  temperature)  for  70  s. 


Figure  6. 16  shows  (he  XRD  pattern  of  the  film  after  annealing  at  500°C  for  70  s. 
The  spotty  pattern  instead  of  the  ring  pattern  for  the  precursor  film  indicated  a dramatic 
grain  growth  process  at  this  annealing  temperature.  The  reason  for  this  enhanced  grain 
growth  is  the  formation  of  an  In-Sc  liquid  by  ln4Scj  undergoing  a monotetic  reaction  (see 
Figure  5. 1).  Figure  6.17  shows  the  XRD  diffractogram  measured  at  0.5°  incident  angle 
indicated  formation  of  Cu  or  Cunlm  intermetallic  and  high  temperature  Cu2.xSe  phases  at 
the  near  surface  region.  Figure  6.18  shows  an  integrated  XRD  spectrum  measured  at  6" 
incident  angle. 


Figure  6. 1 6 X-ray  diffraction  pattern  at  6“  incident  angle  of  the  film  after  RTP  at  500“C 
(TC  temperature)  for  70  s. 


149 


Figure  6. 1 7 XRD  spectra  at  0.5°  incident  angle  of  thin  81m 
after  RTP  for  70  s at  500°C  (TC). 


10  30  50  70 


2:CuInSe- 

3:Cu:Se 


2©  (degree) 


Figure  6. 18  X-ray  diffraction  spectra  at  6°  incident  angle  of 
the  him  after  RTP  at  500°C  (TC  temperature)  for  70  s, 
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Morphology 

Plane-view  (a,  c,  e)  and  cross-sectional  (b.  d.  and  f)  SEM  photographs  of  the 
precursor  and  RTP  films  are  given  in  Figure  6.19.  The  plane-view  SEM  photograph 
shows  no  change  in  the  surface  morphology  after  the  290°C  (TC  temperature).  70  s heat 
treatment  (Figure  6.19(c)).  A significant  change  was  observed  in  the  cross-sectional 
photograph  in  Figure  6.19(d).  These  results  suggested  an  interface  reaction  between  the 
precursor  layers,  and  the  reaction  did  not  extend  to  the  top  surface  region.  Dramatic 
changes  were  observed  after  the  S00°C  (TC  temperature),  70  seconds  heat  treatment  in 
both  plane-view  and  cross-section  photographs  as  shown  in  Figures  6.19(e)  and  (f).  The 
surface  morphology  suggests  a liquid  phase  formed  during  the  rapid  thermal  process. 
This  could  be  attributed  to  the  melting  of  lmiSej  phase,  which  undergoes  a binary 
monotetic  reaction  at  S24°C  (sec  Chapter  S).  This  is  an  expected  result  since  we  already 
know  the  temperature  on  the  surface  is  about  200°C  higher  than  the  TC  temperature  at 
this  ramp  rate  as  indicated  by  our  temperature  calibration  studies. 


Figure  6.1 
films. 


ction  SEM  photographs  of  the  precursor  and  RTP 
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The  overall  composition  of  Ihe  samples  was  analyzed  by  (ICP)  spectrometry  The 
estimated  uncertainty  for  ICP  analysis  is  ±0.5%  relative.  Major  sources  for  uncertainty 
are  variability  in  the  plasma,  variability  in  the  pump  rate,  and  variability  in  the  nebulizer 
efficiency  The  precursor  has  an  overall  Composition  ofCu:  In:  Se  = 35.6:  25.2:  39.1.  The 
AES  survey  spectrum  for  the  as  received  precursor  is  given  in  Figure  6.20.  The  intense 
oxygen  signal  indicated  the  precursor  surface  is  mainly  covered  by  native  oxide.  The 
surface  composition  estimated  by  AES  is  given  in  Table  6.2.  The  peak  intensity  depth 
profile  data  are  shown  in  Figure  6.21.  The  oxygen  signal  decreased  drastically  after 
sputtering  for  a shot  time  The  bottom  layer  composition  as  estimated  by  AES  is  given  in 
Table  6.3.  The  depth  profile  is  consistent  with  the  XRD  data,  which  suggested  the 
formation  of  a three  layer,  litiSej/CuInSej/Cu^Se,  structure. 

Table  6.1.  ICP  chemical  analysis  of  the  investigated  samples 


RTP  290 

The  overall  composition  of  the  films  change  a little  after  RTP  at  290°C  (TC)  for 
70  s as  shown  by  the  ICP  data  in  Table  6. 1 , However,  the  AES  data  shown  in  Figure  6.23 
did  show  a dramatic  change  after  RTP.  The  depth  profile  clearly  showed  a deep  Cu 
diffusion  from  the  backside  of  the  film  into  the  top  ln.Sc)  layer.  This  indicated  the 
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reaction  front  between  CuiSe  and  IiuScj  layer  was  moving  forward  and  forming  CulnSe2 
during  the  RTP.  This  is  consistent  with  the  XRD  data,  which  showed  the  disappearance 
of  XRD  iil-iSei  peaks  and  increasing  CuInSe?  peak  intensity.  However,  the  reaction 
mechanism  is  mainly  through  solid  state  dittusion. 

RTP  500 

Again  the  overall  composition  of  the  film  change  a little  after  RTP  at  500°C  (TC) 
for  70  s as  shown  by  the  ICP  data  in  Table  6. 1 . This  clearly  demonstrates  an  advantage  of 
using  binary  layers  as  precursors  as  compared  to  slacked  elemental  layer  structures.  The 
AES  survey  spectrum  on  the  as  received  sample  is  given  in  Figure  6.24.  The  strong  Cu 
signal  indicated  the  reaction  had  reached  the  surface  of  the  film  during  the  RTP  The  AES 
depth  profile  given  in  Figure  6.25  showed  a uniform  distribution  of  Cu,  In  and  Se  atomic 
species.  Such  fast  mass  transport  is  believe  to  be  a result  of  liquid  phase  formation  as 
IniSei  underwent  a monotetic  reaction,  as  supported  by  the  XRD  and  SEM  dala. 
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received  precursor  film. 


Figure  6.20  AES  survey  scan  on  as  i 


Table  6.3  Precursor  film  composition  measured 
by  AES  after  100  min  sputtering 


Figure  6.21  AES  survey  scan  on  precursor  film  after  100  min  sputtering  etching 
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Figure  6.22  AES  depth  profile  of  the  precursor  film. 
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Figure  6.23  AES  depth  profile  of  the  film  after  RTP  at  290°C. 
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Figure  6.24  AES  survey  scan  of  the  film  surface  after  RTP  at  S00°C. 


1H013H  W3d 

Figure  6.25  AES  depth  profile  of  the  film  after  RTP  at  500”C. 
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6 8 Results  and  Discussion-Precursor  Siniclurc-ll 


The  composition  analysis  for  the  precursor  and  the  reference  binary  layers  is 
shown  in  Table  6.4.  The  WDS  analysis  is  in  agreement  with  the  ICP  data  for  the  binary 
samples.  For  the  bi-layer  precursor  structure,  WDS  is  not  able  to  provide  the  correct 
overall  composition.  The  ICP  data  show  the  overall  precursor  structure  is  a slightly  Cu- 
rich.  The  overall  composition  should  be  adjusted  by  increasing  the  ln-Se  binary  layer 
thickness. 

The  XRD  spectra  are  shown  in  Figures  6 26  and  6.27  for  the  ln-Se  and  Cu-Se 
reference  samples  respectively.  The  ln-Se  layer  primarily  consists  of  small  crystallites, 
which  showed  a broad  peak  in  the  XRD  spectrum.  The  Cu-Se  layer  mainly  consists  of  the 
a-CuSe  phase.  The  XRD  spectrum  for  the  bilayer  precursor  is  given  in  Figure  6.28.  The 
increased  peak  height  at  20=26.7  value  and  the  small  extra  peak  20=44.2  indicates  a 
small  amount  of  pre-reaction  to  form  CulnSc2  during  the  deposition  process.  A schematic 
drawing  for  the  precursor  structure  is  shown  in  Figure  6.29.  The  AES  depth  profile  of  the 
precursor  film  given  in  Figure  6.30  is  consistent  with  this  structure. 


Table  6.4.  Composition  analysis. 


InSc/ViFC  Scbl  :r  Cu.-26.78  | 1 


Figure  6.26  XRD  spectrum  for  In-Se  binary  reference  sample. 


Figure  6.27  XRD  specirum  for  Cu-Se  binary  reference  sample. 


Figure  6.28  XRD  spectrum  for  the  bilayer  precursor. 


Figure  6.29  Schematic  drawing  of  the  bilayer  precursor. 
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Figure  6.30.  AES  depth  profile  of  the  precursor  structure-11 


A sel  of  experiment  was  performed  lo  understand  the  influence  of  ramp  rate. 
Three  ramp  rales  1 , 2 , and  10  °C/s  were  used  in  the  study.  The  slow  ramp  rate  resulted  in 
bubble  formation  and  de-lamination  of  the  films.  A plane-view  optical  micro-graph  is 
shown  in  Figure  6.3 1 . The  mechanism  can  be  easily  understood  by  examining  the  Cu-Se 
binary  phase  diagram  shown  in  Figure  6.3.  The  CuSe  phase  decomposes  to  from  Cli;,Se 
phase  and  Sej  vapor  during  the  heating  process  before  forming  the  CuSe  liquid  phase 
The  production  of  Se2  vapor  causes  bubble  formation  and  thus  de-lamination  of  the  film. 

At  higher  ramp  rate,  it  is  expected  that  the  single  liquid  phase  domain  will  appear 
more  quickly.  Although  the  Se  vapor  pressure  is  very  high  (Figure  4.5)  in  this  domain  the 
mass  transport  and  reactivity  arc  also  much  higher.  On  the  other  hand,  at  the  slower  ramp 
rate  most  of  the  a-CuSe  phase  has  time  lo  transform  to  Cu;.,Se  before  the  temperature 
can  reach  the  single-phase  domain  Thus,  the  higher  ramp  rate  produced  more  CulnSei 
with  better  crystallinity  and  reduced  de-laminating.  This  is  clearly  shown  in  the  optical 
plane-view  image  of  the  RTP  film  processed  at  a ramp  rate  of  10  'C/s  (Figure  6.32).  The 
XRD  spectra  for  the  10  °C/s,  and  I °C/s  cases  are  shown  in  Figure  6.33  and  6.34, 
respectively.  The  XRD  spectra  indicate  the  peaks  for  the  10  ”C/s  ramp  rate  RTP  sample 
have  higher  intensity  than  the  corresponding  XRD  peaks  processed  at  a ramp  rate  of 
1'C/s. 

A reference  a-CuSc  binary  film  was  subjected  to  RTP  treatment  at  a ramp  rate  I 
'C/s  for  70  s to  further  confirm  this  conclusion.  WDS  analysis  indicated  the  film 
composition  changed  from  Cu:50.33,  Se:49.46  to  Cu:60,78,  Se:  39,22.  XRD  analysis 
indicated  the  film  has  changed  to  the  Cu:.xSe  phase  (Figure  6.35).  These  results  clearly 
show  the  Se*  vapor  evaporated  from  tx-CuSe  during  the  RTP  process  to  produce  Cuj.xSe. 
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The  ICP  analysis  of  the  bi-layer  samples  after  RTP  treatment  is  given  in  Table  6.5.  In 
contrast  to  the  result  for  the  RTP  experiment  of  reference  a-CuSc  sample,  the  ICP  data 
indicated  the  overall  film  composition  did  not  change  much  after  the  bi-layer  precursor 
was  subjected  to  RTP  treatment. 

These  results  are  encouraging  that  our  precursor  structure  design  is  feasible. 
Putting  the  InSe  layer  on  top  of  the  CuSe  layer  to  minimize  the  Se„  evaporation  loss 
appears  to  be  effective  Figure  6.36  shows  the  AES  depth  profiles  of  the  film  after  RTP 
treatment  with  a ramp  rate  of  IO°C/s  and  annealed  at  470°C  (TC  temperature)  for  10  s A 
Uniform  composition  profile  was  achieved  throughout  the  film.  The  high  mass-transport 
and  reaction  rate  is  believed  to  be  a result  of  liquid  phase  formation  in  the  CuSe  layer. 
Plane-view  SEM  photographs  of  the  precursor  film  and  the  film  after  RTP  annealing 
(10°C/s)  are  shown  in  Figure  6.37a  and  b,  respectively.  The  SEM  images  show  the 
formation  of  crystals  with  grain  size  on  the  order  of  several  microns. 

Table  6.5  Composition  analysis  after  RTP  treatment. 


CuSc/InSc  bi-laycr  after  RTP  ( 1"C  / s)  ICP  Cu:2T6l>  I ln:23.54  I Sc:4M6 


Figure  6.31  Plane-view  optical  micro-graph  at  the  ramp  rate  (I  “Os). 


Figure  6.32  Plane-view  optical  micro-graph  at  ramp  rale  (to"C/s) 


Figure  6.33  XRD  spectrum  for  RTP  CuInSe2  film  with  ramp  rate  10 "( 
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Figure  6.34  XRD  spectrum  lor  CuSe/lnSe  bi-layer  film  after  RTP  with  ramp  rate  I T/S 


Figure  6.35  XRD  spectrum  for  a-CuSe  film  after  RTP  with  ramp  rate  1 "C/s. 
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Figure  6.36  AES  depth  profile  of  the  precursor  structure  (II) 


Figure  6.37  The  plane-view  SEM  photographs  of  the  precursor  (a)  and  RTP  (b)  films. 


Three  different  CulnSe,  layer  growth  mechanisms  are  proposed  in  Figure  6.38, 
6.39,  and  6.40.  The  first  growth  mechanism  is  solid  state  inter-diffiision  between  IrpSe, 
and  CujSe  binary  layers  (Figure  6.38).  The  second  growth  mechanism  is  a liquid-solid 
reaction  due  to  the  melting  of  ln,Sej  Coarse-grained  CulnSe;  crystals  grow  from  the 
liquid-solid  mixture.  Uniform  mass-transport  was  achieved  in  short  time  (70  s)  as  shown 
in  the  AES  depth  profile  (Figure  6.25).  The  main  phases  to  appear  in  the  film  after  500°C 
RTP  are  CulnSe,.  Cu,Se  and  a small  amount  of  a Cu-ln  inter-metallic  phase  The  third 
mechanism  is  a liquid-solid  reaction  due  to  the  melting  of  ct-CuSe  in  the  bottom  layer 
Coarse-grained  single  phase  CulnSe,  crystals  were  synthesized  from  the  liquid-solid 
mixture  as  suggested  in  the  reaction  pathway  analysis  (section  6.3).  Uniform  mass- 
transpon  was  achieved  in  short  time  (70  s)  as  shown  in  the  AES  depth  profile  (Figure 
6.36). 

In  this  study  an  alternative  strategy  for  low  temperature  precursor  fabrication  of 
CIS  films  has  been  proposed  and  tested  using  ex-situ  rapid  thermal  processing.  The 
experiments  showed  that  large  grain  CIS  could  be  synthesized  through  this  approach 
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Figure  6.38  Illustration  of  the  solid-state  diffusion  growth  model. 
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Figure  6.39  Illustration  of  the  liquid  assisted  growth  mechanism  for  precursor 
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CHAPTER  7 

X-RAY  ABSORPTION  FINE  STRUCTURE  (XAFS)  INVESTIGATIONS  OF  COPPER 
INDIUM  SELENIDE  SEMICONDUCTOR  ALLOYS 

CIS-based  solar  cells  have  been  thought  to  be  p-n  structures  (p-CIS,  n-CdS). 
However,  such  a model  does  not  predict  the  failure  of  high-Ga-content  devices.  A novel 
model  has  been  proposed  by  Schmid  el  at.  [Sch  93]  that  the  growth  kinetics  of  PVD- 
grown  CulnSet  result  in  the  in-situ  formation  of  a surface  layer  of  CulnjSes  and  that  this 
CuInSet/CulnjSej  helerojunction  is  the  true  heterojunction  partner  in  CIS  devices.  This 
surface  material  is  reported  to  show  compositions  near  CuIntSe;  and  has  been  referred  to 
as  an  "Ordered  Vacancy  Compound"  (OVC).  This  significant  finding  resembles  earlier 
observation  of  a buried  junction  elucidated  by  electron-beam  induced  current  (EBIC) 
[Mat  86]  measurements  However,  the  crystal  structure  of  CuhtjSej  and  its  phase  relation 
with  CulnSe:  is  still  not  dear  at  this  point  In  this  chapter  a detailed  study  of  the  local 
structure  of  CulnSe:  and  CulniSej  by  X-ray  Absorption  Fine  Structure  was  reported 


7,1  Introduction 

X-ray  absorption  fine  structure  (XAFS)  spectroscopy  refers  to  the  measurement  of 
the  modulation  of  the  X-ray  absorption  coefficient,  p.  as  a function  of  photon  energy,  E, 
at  and  above  the  X-ray  absorption  edge.  One  of  the  major  advantages  of  XAFS  is  its 
atomic  selectivity,  which  is  capable  of  probing  the  local  structure  of  each  atomic 
constituent  of  a sample.  Conventionally,  EXAFS  spectra  refer  to  the  region  40  to  1000 
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eV  above  the  absorption  edge.  The  pre-edge  region  contains  bonding  information  such  as 
the  electronic  configuration,  and  the  site  symmetry.  In  between  the  pre-edge  and  the 
EXAFS  regions  is  the  X-ray  Absorption  Near-Edge  Structure  (XANES)f  which  contains 
information  about  the  local  electronic  structure  and  geometrical  configuration  of  the 
materials.  Figure  7 1 shows  a typical  X-ray  absorption  spectrum  around  an  absorption 


Figure  7. 1 A typical  XAFS  spectrum. 


7.2  Theory  of  EXAFS 

If  A,  and  / represent  (Ire  imonsiiy  of  the  incident  X-ray  beam  and  transmitted  X-ray 
beam  respectively,  the  relation  between  h,  and  / is  given  by 

/ = /«  «'>',  (7.|) 

In  this  equation  p(£)  is  the  X-ray  absorption  coefficient,  E is  the  photon  energy,  and  x is 
the  sample  thickness.  The  term  p(£)  is  a monotonically  decreasing  function  of  E in  most 
energy  ranges.  When  the  photon  energy  equals  the  binding  energy  of  a core  electron,  a 
sudden  increment  on  p due  to  the  photoelectric  absorption.  This  energy  is  denoted  as  £>: 
the  threshold  energy.  When  the  electrons  are  in  the  n =1  shell  the  absorption  edge  is 
called  the  A'-edge.  For  the  n =2  shell,  the  corresponding  edges  arc  called  the  /.-edges  In 
the  photoelectric  process,  the  energy  of  an  X-ray  photon  is  absorbed  by  a core  electron 
and  promoted  to  an  unoccupied  valence  level  or  to  continuum.  The  momentum  of  the 
emitted  photoelectron  can  be  determined  by  the  relation 


where  p is  the  momentum  of  the  photoelectron,  h is  Planck's  constant,  m is  the  electron 
mass,  and  v is  the  X-ray  photon  frequency. 

For  an  isolated  atom  the  photoelcctron  electron  can  be  represented  as  an  outgoing 
wave  as  shown  in  Figure  7.2  by  the  solid  lines.  When  the  absorbing  atom  is  surrounded 
by  other  atoms,  as  in  any  molecule  or  condensed  phase,  the  outgoing  photoelectron  wave 
is  backscattered  by  the  surrounding  atoms  as  indicated  by  the  dashed  lines  in  Figure  7.2. 
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The  final  slate  is  the  superposition  of  the  outgoing  and  the  scattered  waves.  The 
backscattcred  waves  will  add  or  subtract  from  the  outgoing  wave  at  the  center  depending 
on  their  relative  phase.  The  total  amplitude  of  the  electron  wave  function  will  be 
enhanced  or  reduced,  thus  modifying  the  probability  of  absorption  of  the  X-ray.  The 
absorption  coefficient  p is  proportional  to  the  atomic  absorption  cross  section  o.  If  we 
denote  Pas  the  probability  of  an  atom  absorbing  a photon  of  energy  hv,  from  the  Fermi's 
Golden  rule  within  the  dipole  approximation  gives: 

P = 4itah<oZ\(i\er\f'f{s(E, -E,  -hoi)  (7.3) 

where  l/>  and  \f>  represent  the  initial  and  final  states  with  E,  and  E/  as  the 
corresponding  energies,  respectively,  a = — is  the  hypcrfine  structure  constant,  ft  is  the 
Planck  constant  divided  by  2it,  e is  the  electronic  charge,  c is  the  speed  of  light,  a is  the 
electric  field  polarization  vector  of  the  X-ray  photons,  r is  the  position  vector  of  the 
scatter  with  respect  to  the  absorber  atom,  and  to  is  the  angular  frequency  of  the  photon. 

In  EXAFS  we  are  only  interested  in  the  oscillatory  part  of  the  total  absorption.  A 
simple  way  of  deriving  the  EXAFS  equation  is  by  considering  one-electron  single 
scattering  in  plane  wave  approximation.  With  these  approximations  the  oscillatory  part 
X(*)  of  the  total  absorption  coefficient,  p(*),  can  be  described  by  the  summation  of  sine 
waves  originating  from  all  neighboring  atoms,  and  the  EXAFS  equation  is  defined  as 
[Say  87], 

*<*)  = I^So^,(*)exp(-2*V(!)exp(--^)sin[2*R,  +£,(*)]  (7.4) 


Here y refers  10  the  yih  coordination  shell,  /(,  is  the  average  interatomic  distance,  8 ft)  is 
the  phase  shift  experienced  by  the  electron  in  the  scattering  process.  N,  is  the  average 
coordination  number  in  the  jth  shell  and  a,  is  a Debye-Waller  type  of  term  which 
represents  the  root  mean  square  fluctuation  in  R,  caused  by  thermal  motion  of  the  atoms 
and  the  structural  disorder,  h',(k)  is  the  back-scattering  amplitude,  is  an  amplitude 
reduction  factor  which  accounts  for  the  photoelectron  energy  loss  due  to  many  body 
effects  and  shake-up/shake-off  processes,  and  the  term  accounts  for  the  finite 
lifetime  of  the  excited  state. 


Figure  7.2  Outgoing  and  backscattered  portions  of  the  final  state  wave  function 
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Effects  of  Thermal  Vibration  and  Sialic  Disorder 

The  Debye-Waller  factor  o2  in  the  EXAFS  equation  has  two  contributions:  the 
thermal  vibration,  c a,,  and  static  disorder,  a\  It  is  important  to  note  that  the  o in  the 
EXFAS  equation  (7  4)  is  different  than  the  Debye-Waller  term,  p2.  in  diffraction  The  p2 
term  in  diffraction  accounts  for  the  mean  square  displacement  of  the  atom  from  its  lattice 
site,  while  for  EXAFS,  o2,  is  the  relative  deviation  between  the  absorbing  and 
backscartering  atoms.  Thus  in  EXAFS  the  vibration  modes  causing  the  atomic  pairs  of 
interest  moving  in  the  same  phase  and  amplitude  have  no  contribution  to  o2.  <r*  could  be 
calculated  from  first  principle  [Dim  98]  or  estimate  using  the  Einstein  or  Debye 
approximations.  The  Einstein  model  assumes  the  single  optical  vibration  frequency,  be, 
and  is  given  by 


(7.5) 


where  M,  is  the  reduced  mass  of  absorbing  atom  and  backscattering  atom  pair. 
In  the  Debye  model,  the  a2n,(T)  has  a simple  form  given  by 


Debye  temperature,  N is  the  number  of  atoms,  and  V is  the  crystal 


The  structural-disorder  Debye  Waller  factor.  cr„  measures  the  degree  of  deviation  of  the 
equilibrium  inter-atomic  distance  of  each  absorber  and  backscatter  from  the  average  bond 
length.  R.  is  given  by. 


(7.8) 


where R,  is  the  equilibrium  bond  length  of yth  atom  in  that  shell. 


7 3 Experimental  XAFS  Apparatus 

An  X-ray  source.  X-ray  optics,  and  the  X-ray  detectors  are  the  primary  experimental 
system  for  XAFS  experiments.  A schematic  diagram  for  the  experimental  setup  is  shown 
in  Figure  7.3. 


Figure  7.3  Schematic  diagram  for  XAFS  experiment 
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7,3.1  X-ray  source 

A high  intensity  X-ray  source  with  a continuous  energy  spectrum  is  necessary  for 
XAFS  experiments.  The  principal  X-ray  source  used  for  XAFS  is  that  from  an 
accelerating  charge  particle  emits  radiation.  Synchrotron  radiation  is  the  most  widely 
used  X-ray  source  for  XAFS  experiments.  All  the  measurements  in  this  study  were  done 
at  the  Materials  Research  Collaborative  Access  Team  (MRCAT)  beam  line,  Advanced 
Photon  Source,  at  the  Argonne  National  Laboratory. 

The  beam  acceleration  and  storage  process  begins  at  the  electron  gun,  which 
emits  electrons  that  exit  the  gun  at  100  keV.  Electrons  are  raised  to  an  energy  level  of  200 
MeV  by  a series  of  accelerating  structures  in  the  first-stage  linear  accelerator  (linac). 
Next,  the  480W  beam  of  electrons  strikes  a 7-ram-thick  tungsten  positron  conversion 
target,  which  produce  electron-positron  pairs.  Positrons  are  accelerated  to  4S0  MeV  by 
the  second-stage  linac.  The  368-m-tong  booster  synchrotron  raises  positron  energies  at  a 
rate  of  32  keV  per  turn,  in  0.2S  s positrons  orbit  the  booster  200000  limes  as  their  energy 
climbs  to  7 GeV.  The  positrons  are  then  injected  into  the  1014-m-circumference  storage 
ring.  The  positron  beam  orbits  the  storage  ring  more  than  271000  times  per  second, 
steered  and  focused  by  1097  powerful  electromagnets  as  it  circulates  The  beam 
decelerates  at  a rate  of  about  6 MeV  per  turn  as  it  emits  synchrotron  radiation  The 
energy  loss  is  replaced  by  the  352-MHz  storage  ring  rf  systems. 

Every  10  to  20  hours,  the  storage  ring  is  refilled  with  positrons.  A storage  ring  is 
actually  a set  of  curves  connected  by  straight  sections.  Linear  arrays  of  north-south 
permanent  magnets  with  alternating  polarity  are  inserted  into  the  straight  section  (hence 
the  name  Insertion  Device  (ID)),  one  array  above  the  beam  path,  the  other  below.  When 
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where  E is  the  X-ray  energy,  p,  is  the  eleciron  density  The  fundamental  will  be  totally 
reflected  while  the  harmonics  will  be  absorbed  by  the  mirror  by  setting  the  mirror  angle 
between  0C  of  the  fundamental  and  9,.  of  the  harmonics. 

7.3.3  X-ray  detector 

The  most  common  way  of  measuring  X-ray  intensity  in  XAFS  is  by  using  a gas 
ionization  detector'  or  an  "ion  chamber"  The  basic  structure  of  an  ion  chamber  is  a pair 
of  parallel  metal  plates.  A certain  type  of  gas,  such  as  He,  Ar,  N;  or  Kr,  or  a mixture  of 
these  gases  fills  the  chamber.  The  choice  of  gas  depends  on  the  energy  range  of  interest. 
A DC  voltage  is  applied  between  the  plates.  The  X-ray  will  ionize  the  gas  molecules  in 
the  chamber  The  ions  (electrons)  will  transport  to  the  negative  (positive)  collecting  plate 
generating  a current  in  a closed  circuit.  Figure  7 4 shows  a photograph  of  the 
experimental  set-up. 
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7 4 Composition  Effects  on  the  Local  Structure  of  CulnSe; 


The  crystal  structure  of  a-CulnSej  has  been  determined  and  refined  by  several 
authors.  It  belongs  to  A!BmC'1  family  with  the  chalcopyrite  (Eli)  structure  (space  group 
1 4 2d),  with  Cu  in  site  4a  at  (0  0 0),  In  in  site  4b  at  (0  0 1/2),  and  Se  in  site  8d  at  (x  1/4 
1/8)  [Spi  74].  It's  a ternary  extension  of  the  Grimm-Sommcrfcld  compounds,  which  are 
characterized  by  tetrahedral  bonding.  Each  anion  C is  coordinated  by  two  A and  two  B 
cations,  whereas  each  cation  is  tetrahedrally  coordinated  by  four  anions  C.  The  anion  C 
usually  adopts  an  equilibrium  position  closer  to  one  pair  of  cations  than  to  the  other, 
Rac*Rbc  (bond  alternation).  The  nearest  neighbor  bond  lengths  for  a perfect  crystal  can 
be  derived  from  the  structural  parameters  and  given  by 

RAC  =[*!  + (l+l/2)/16),,Io  (7.11) 

Rk  =[(x-i)’+(l  + 7:)/l«r!a  (7.12) 

’!  m c 2a  (tetragonal  deformation)  is  the  ratio  between  the  lattice  parameters  c and  a . A 
schematic  illustration  of  chalcopyrite  CulnSe;  unit  cell  is  shown  in  Figure  7,5. 

An  interesting  question  was  raised  by  recent  studies  by  Merino  el  ai  [Mer  96, 
Mer  98].  They  reported  the  Se  atomic  coordinate,  v(Se).  thus  the  Cu-Se  and  In-Se  inter- 
atomic distances,  depend  strongly  on  the  Cu  occupation  number  based  on  X-ray  powder 
diffraction  and  Rietveld  refinement.  Using  their  results  they  suggested  that  the  spread  in 
reported  energy  gap  value.  E„  in  the  literature  could  be  due  to  the  variation  of  the  x 
values,  predicted  by  Jaffe  and  Zunger's  model  [Jaf  84],  A direct  measurement  of  the  band 
gap  of  their  samples  (Mer  98],  however,  show  that  there  is  no  correlation  between  the 


measured  gap  and  ihe  suggcsied  x values  In  fact,  for  some  samples  E,  decreases  as  x 
increase,  in  contrast  to  the  Jaffe  and  Zunger's  model  [Jaf  84).  In  this  work.  Cu  K-edge 
X-ray  Absorption  Fine  Structure  (XAFS)  measurements  were  performed  on  Cu,In,Se, 
crystals  with  different  compositions  close  to  the  CulnSej  stoichiometry  and  CulnsSex. 


The  samples  were  prepared  from  powder  CuiSe  and  InjScj  in  boron  nitride  (BN) 
coated  evacuated  quartz  ampoules.  The  BN  coating  was  used  to  prevent  copper 
devitriiication  of  the  quartz  at  high  temperature  [Shu  93).  These  mixtures  were  then 
slowly  heated  to  the  melting  temperature  of  CujSe  then  gradually  cooled  to  700°C  and 
were  kept  in  the  furnace  at  700°C  for  one  week  to  assure  equilibrium.  One  sample  was 
then  quenched  in  liquid  nitrogen  and  another  one  is  slowly  cooled  to  room  temperature. 
The  third  sample  was  heated  above  the  order-disorder  temperature  (8I0°C)  to  900°C,  and 
held  there  for  a week  then  quenched  in  liquid  nitrogen.  Electron  Probe  Micro-Analysis 
(EPMA)  (JEOL  Superprobe  733  with  wavelength  dispersive  spectrometer)  was 
performed  to  assure  the  composition  uniformity  within  of  the  samples.  Powder  X-ray 
Diffraction  (XRD)  was  used  to  verify  the  formation  of  chalcopyrite  CuInSe;  and  the 
absence  of  other  phases.  The  XAFS  specimens  were  prepared  by  grinding  the 
polycrystalline  material  The  fine  powder  was  uniformly  spread  over  Scotch  tape.  The 
XAFS  spectra  were  measured  in  the  fluorescence  mode  at  room  temperature 

The  samples  were  analyzed  by  ICP  emission  spectrometry  for  the  atomic  ratios  of 
CIS 
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Figure  7 5 Schematic  diagram  of  chalcopyrite  structure  for  CuInSe; 
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7.4.2  Results 

The  ICP  analysis  of  thecomposition  of  three  CulnSc-  samples  show  that  each 
sample  has  a different  overall  composition  as  given  in  Table  7.1.  including  sample  III 
with  nearly  ideal  1:1:2  Cu/ln/Se  ratio  The  estimated  uncertainty  for  ICP  analysis  is 
±0.5%  relative. 

Table  7. 1 Chemical  analysis  of  the  investigated  samples 

I Cu  at.%  I In  at.%  I Se  at.% 

CulnSc- 1 I 26.4  1 25,8  I 47.8 

CulnSe;  II  26,3  26  3 j 47  4 

CulnSe;  in  I 25  2 I 25.1  I 49.8 


The  X-ray  powder  diffraction  was  used  to  verily  the  average  long-range  order 
structure  The  resulting  spectra  for  all  three  CulnSe:  (I.  11.  Ill)  samples  showed  typical 
chalcopyrite  CulnSe:  powder  XRD  patterns  [Mer  96], 

The  three  samples  were  taken  to  Argonne  national  laboratory  and  X-ray 
absorption  spectra  were  taken  using  synchrotron  radiation  from  the  Advance  Photon 
Source.  The  XAFS  data  were  analyzed  by  using  the  WinXAFS  package  [Res  97],  The 
collected  data  were  processed  by  standard  methods  [Say  87].  The  experimental  intensity 
data  were  first  converted  to  an  absorption  coefficient  based  on/r(E)x  = ln(l„/l)  The  Cu 
K-edge  spectrum  of  CulnSe,  sample-1  is  shown  in  Figure  7.6(a).  The  pre-edge 
backgiound  was  then  removed  form  the  spectra.  The  next  step  was  to  convert  the  photon 
energy  E to  the  photoelectron  wave  vector  k using  the  relation  k = ^2m(E-E„)/h: , 
where  m is  the  electron  mass.  E is  the  X-ray  photon  energy,  and  Eo  is  the  K-edge  energy. 


The  normalized  XAFS  spectra  were  obiained  by  subtracting  the  background  fio(k)  from 
the  measured  absorption  coefficient  fi(k)  and  were  normalized  by  the  edge  jump  Ago(O): 


(7  13) 


The  resulting  xW  oscillations,  shown  in  Figure  7 6(b),  were  then  Fourier  transformed 
into  r space  with  a*-’  weighting  and  k-space  ranges  of  3-12  A'1.  The  Fourier-  transformed 
Cu-K  edge  r-space  data  are  shown  in  Figure  7.6(c).  To  isolate  the  nearest  neighbor 
information,  data  were  inverse  Fourier  transformed  with  a square  window  range  of  1 .6  to 
2.6  A The  resulting  first  shell  EXAFS  are  shown  in  Figure  7.6(d).  The  same  <t-space  and 
r-space  windows  were  applied  to  the  data  analysis. 

Figure  7.7  shows  the  Cu-K  edge  r-space  data  for  samples  I and  III  The  first  sharp 
peaks  appear  at  about  2.1  A and  are  due  to  Cu’s  first  nearest  neighbors  (Sc).  The  second 
broader  peaks  are  at  around  3.9  A,  and  are  produced  by  scattering  from  the  second 
nearest  neighbors  (Cu,  In).  It  can  be  clearly  see  the  two  spectra  are  almost  identical  This 
indicates  the  local  structures  (i.e,,  the  bond  length,  dcu.se,  and  the  number  of  nearest 
neighbors.  N)  around  Cu  atom  in  these  two  samples  are  nearly  identical. 


Figure  7 6 (a)  Absorption  as  a (unction  of  photon  energy  at  the  Cu-K  edge  in  CulnSe;  (I) 
at  298  K.  (b)  Cu-K  edge  EXAFS  oscillations,  /r%(k),  as  a (unction  of  k after  background 
removal,  (c)  Fourier  transform  of  (b)  to  real  space.  The  transform  window  is  3-12  A'1 
(d)  Fourier  filtered  first  shell  EXAFS 


Figure  7 7 Cu-K  edge  EXAFS  in  real  space  transformed  using  a window  k=3  to  12  A for 
CuInSei  samples  1 and  III 
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In  Figure  7.8,  the  Fourier  tillered  EXAFS  on  lhe  Cu-X  edge  For  the  first-neighbor 
peaks  of  samples  1, 11,  111  are  shown.  The  first  shell  EXAFS  from  these  three  samples  are 
identical  This  indicates  the  first  nearest  neighbors  around  Cu  atoms  are  identical  in  these 
three  samples  To  obtain  quantitative  information,  the  standard  EXAFS  equation  (7.4) 
was  used  (Say  87]  in  the  least  square  analysis. 

To  fit  the  data  to  this  equation,  the  structural  parameters  of  chalcopyrite  CulnSe: 
are  needed  for  obtaining  the  backscattering  amplitude.  Ffk),  and  phase  shift,  Sfflt),  from 
the  CulnSe2  EXAFS  spectrum  A review  of  the  reported  structural  data  for  chalcopyrite 
CulnSe2  is  listed  in  Table  7.2,  The  agreement  in  cell  parameters  is  good,  but  a 
considerable  difference  is  observed  in  the  reported  Se  atomic  coordinate,  x(Se),  and  the 
resulting  bond  lengths,  do,*,  and,  di„.S[  The  values  from  Spiess  el  al.  [Spi  74]  dd-s.- 
2.424  A.  and  di„*,=  2.598  A are  accepted  here  tor  chalcopyrite.CuInSe2  structure. 


Table  7.3  Structural  parameters  of  CulnSe2  determined  by  X-ray  diffraction 
Type  a (A)  c lAl  X (Se)  d .(A)  d ..  JAI  Reference 


Single  Crystal  5.784  11.616  0.224 


The  backsc 


nplitude.  Ft(k),  and  phase  shift,  &lk).  for  Cu  K-edge 


extracied  from  sample  III  here  using  Cu-Se  inter-aiomic  dislance  equal  to  2.424  A These 
functions  were  then  used  to  calculate  the  Cu-Se  inter-atomic  distances  of  sample  I and  II. 
The  least  square  fitting  results  are  given  in  Table  7.3.  The  values  of  Cu-Se  inter-atomic 
distance  are  within  0.01  A of  each  other. 


Table  7.3  Least-square  fitting  results 


CulnSer-l 

CulnSe.-III 


Coordination  number  I Cu-Se  inter-atomic  distance  A 

4 2.424  *-0.002 

4 2.425+- 0.002 

4 2424 


The  samples  examined  in  this  study  did  not  give  evidence  for  a large  density  of  Cu 
vacancies  The  question  whether  the  inter-atomic  distance  is  a function  of  the  Cu 
occupancy  number  remains  unanswered  experimentally,  even  though  a large  density  of 
Cu  vacancies  reported  in  [Met  96]  are  questionable.  This  question,  however,  could  be 
addressed  by  X-ray  absorption  measurements  on  In-rich  samples  such  as  CuInjSes.  The 
CuInjSes  compound  lies  on  the  CuiSe-IniSer  tie  line.  In  general  the  structure  of  this 
compound  belongs  to  defect  tetrahedral  structures.  They  can  be  described  as  normal 
tetrahedral  structures  with  a certain  fixed  number  of  unoccupied  structure  sites.  The 
EXAFS  data  on  this  compound  indicated  that  the  Cu-Se  and  in-Se  inter-atomic  distances 
are  the  same  in  CulnScj,  even  though  20%  of  Cu  sites  are  vacant  in  this  structure.  The 
detailed  EXAFS  studies  of  CuIn3SeS  are  described  in  section  7.5. 


k(A') 


Figure  7.8  Fourier  filtered  EXAFS  on  the  Cu  K edge  for  the  first-neighbor  peaks  of  three 
different  stoichiometric  CuInSei  samples  I,  If,  III 
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Considerable  work  has  been  carried  out  to  determine  the  structure  of  this 
compound.  A summary  of  the  reported  structure  of  this  compound  is  given  in  Table  7.3. 
Palatnik  and  Rogacheva  [Pal  67]  first  reported  the  structure  for  CuInjSe*  and  assigned  the 
crystal  structure  to  the  space  group  of  either  / 4 or  /42m  based  on  the  observed 
extinction  rule  from  X-ray  powder  diffraction  data.  Later,  a refined  structure  solution  was 
given  by  HOnle  et  al.  [Hon  88]  using  single  crystal  X-ray  diffraction  data  They  determined 
the  space  group  was  PA2c  and  proposed  the  name  P-chalcopyrite  for  this  phase.  Tseng 
and  Wert  [Tse  89]  proposed  an  ordered  vacancy  compound  (OVC)  structure  model,  with 
/ 4 symmetry,  from  their  transmission  electron  diffraction  results.  The  proposed  structure 

model,  however,  has  stoichiometry  of  Cu:ln:Se=l:5:8  rather  than  1:3:5.  Xiao  el  al.  [Xia 
94]  studied  epitaxial  CulmSc,  thin  films  by  transmission  electron  microscope.  They 
constructed  a structure  model  by  incorporating  ordered  point  defects  into  the  chalcopyrite 
CuInSe.  structure  to  account  for  the  extra  diffraction  spots.  They  further  suggested 
CuInjSc,  is  an  ordered  defect  chalcopyrite  structure.  Nomura  and  Endo  [Norn  96]  also 
constructed  a structure  model  for  CuImSes  from  the  same  principle.  Their  model, 
however,  is  not  an  Ordered  Vacancy  Compound  (OVC),  although  they  referred  to  it  as  an 
OVC  in  the  reference.  Hanada  el  a/.[Han  96]  investigated  the  crystal  structure  of  CuInjSet 
by  a combination  of  electron  and  X-ray  diffraction.  They  limited  the  space  group  to  / 4 or 
1 42 in  based  on  the  observed  extinction  rule  followed  by  powder  X-ray  data.  They  further 
determined  the  space  group  to  be  /42m  from  a Convergent  Beam  Electron  Diffraction 
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(CBED)  pattern  along  the  [OOI]  zone  axis.  They  proposed  a structure  model  ABIn:Se4  for 
CuImSes,  which  is  represented  by  (Cuoi  ll:)(ln,i1-«i6)inj$e1  (□:  denotes  structure 
vacancy).  They  concluded  that  CuImSe*  is  neither  an  OVC  nor  a defect  chalcopyrite 

Zhang  et  ai  [Zha  97, 98]  proposed  an  ordered  vacancy  compound  structure  model 
from  first  principle  calculation.  They  followed  two  rules:  (!)  minimum  Madelung  energy 
and  (2)  minimum  deviations  from  the  octet  rule.  The  model  suggested  the  OVC  is  a 
weighted  distribution  of  three  types  of  local  tetrahedral  around  each  Se  to  minimize  the 
deviation  from  the  octect  rule:  2Cu+21n  (called  k=  8),  V„+Cu+2In  (called  *= 7),  and 
V„+3In  (called  k= 9).  The  quantity  k denotes  the  sum  of  the  valence  electrons  from  the  Se 
nearest  neighbor  cations.  The  k=7  and  k=9  clusters  must  occur  in  equal  numbers  to  ensure 
charge  neutrality  and  likely  reside  to  each  other  from  mutual  Coulomb  attraction. 
Chalcopyrite  CulnSej  consists  of  100%  of  the  *=8  clusters.  CuInjSej  has  20%  of  the  k=i 
cluster,  and  40%  each  of  the  k=l  and  k= 9 clusters.  They  also  suggested  the  formation  of 
the  defect  pair  ( 2VC,  + In^, ) stabilizes  the  formation  of  an  OVC. 

Clearly  the  structure  of  CulniSe-  needs  further  investigation.  The  disagreement 
among  these  structural  models  arises  from  the  sensitivity  of  the  different  characterization 
techniques.  The  two  questions  in  single  crystal  X-ray  structure  determination  are  whether 
a true  single  crystal  was  used  and  the  limitation  in  unique  determination  of  the  space 
group.  Transmission  Electron  Microcopy  is  capable  of  selecting  a single  crystal  for  study. 
The  intensity  data,  however,  are  not  easy  to  utilize.  Although  the  observed  Selected  Area 
Electron  Diffraction  (SAED)  patterns  in  their  TEM  samples  [Tse  89][Xia  94]  can  be 
explained  qualitatively  by  their  proposed  ordered  vacancy  structure  model,  their  structure 


model  is  not  unique  to  the  observed  SAED  patterns  and  the  intensity  data  were  not  used 
quantitatively  "Does  the  vacancy  in  Culn)Ses  order?  * is  still  an  open  area  for  research  In 
this  work,  the  X-ray  absorption  Fine  Structure  (XAFS)  measurements  were  performed  on 
the  Cu.  In  and  Se  K-cdges  to  study  the  local  structure  of  CuIn,Sei.  The  first-principles 
total  energy  calculated  results  are  compared  and  discussed 


The  samples  were  prepared  from  powder  Cu.Se  and  InjSej  in  boron  nitride  (BN) 
coated  evacuated  quartz  ampoules.  The  BN  coating  was  used  to  prevent  copper 
devitrification  of  the  quartz  at  high  temperature  (Shu  93]  These  mixtures  were  then 
slowly  heated  to  the  melting  temperature  then  gradually  cooled  to  700°C,  and  were  kept 
in  furnace  at  700°C  for  a week  to  assure  equilibrium  then  quenched  in  liquid  nitrogen 
Powder  X-ray  Diffraction  (XRD)  and  Wavelength  Dispersive  X-ray  Spectroscopy  (WDS) 


were  used  to  verify  the  formation  of  desired  homogeneous  single-phase  material  The 
XAFS  specimens  were  prepared  by  grinding  the  polycrystalline  material  then  spread  fine 
powders  uniformly  over  Scotch  tape. 

The  XAFS  measurements  were  performed  on  the  Materials  Research  Collaborative 
Access  Team  (MRCAT)  beamline  of  the  Advanced  Photon  Source  at  Argonne  National 
Laboratory.  The  beamline  optics  incorporates  a Si(lll)  double  crystal  monochromator 
and  a harmonics-rejection  mirror  to  minimize  the  higher  order  harmonic  content  in  the  X- 
ray  beam  The  X-ray  intensities  were  monitored  using  ionization  chambers.  The  XAFS 
spectra  were  measured  in  the  transmission  mode. 

7.S.2  Data  Analysis 

The  XAFS  data  was  processed  in  the  standard  procedures  as  described  in  section 
7.4.2.  The  backscattering  amplitude,  F,(kJ.  and  phase  shift,  S/k),  for  Cu  and  In  K-edgc 
were  extracted  from  the  CuInSe?  spectra  using  the  known  structural  data.  The  ah  inno 
multiple-scattering  code  FEFF7  [Ank  95]  was  used  to  calculate  F/k),  5,(k).  and  Mk)  for 
Se  K-edgc  data  analysis. 


Figure  7.9  (a)  Absorpiion  as  a function  of  photon  energy  at  the  Cu-K  edge  in  CulnSe,  at 
298  K.  (b)  Cu-K  edge  EXAFS  oscillations,  k* yfk),  as  a function  of  k after  background 
removal,  (c)  Fourier  transform  of  (b)  to  real  space.  The  transform  window  is  3 to  12  A'1. 


The  Cu-K  edge  EXAFS  spectra  and  FT-R  space  data  for  CulniSes  are  shown  in 
Figure  7.9,  The  Erst  FT  peak  appears  at  about  2.1  A is  due  to  the  Erst  nearest  neighbors 
(Se).  The  second  broader  peak  around  3.9  A arises  from  interaction  with  second  nearest 
neighbors  (Cu,  In).  The  isolated  first  shell  spectra  for  CuInSes  and  CuhuSe,  are  given  in 
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Figure  7.1 1.  li  can  be  clearly  seen  lhai  ihe  two  spectra  arc  almost  identical,  except  the  one 
for  CulniSe.  has  a slightly  larger  damping  coefficient.  This  indicates  the  local  structures 
(i.e.,  the  bond  length,  dcvsc,  and  the  number  of  nearest  neighbors,  N)  around  the  Cu  atoms 
in  these  two  semiconductor  alloys  are  nearly  the  same.  To  obtain  quantitative  data  a least- 
squares  fit  to  the  standard  EXAFS  equation  was  performed 

The  least  square  fits  of  the  k'-weighted  data  at  298  K are  shown  in  Figure  7 12. 
The  backscattering  amplitude,  Fjfk).  and  phase  shift,  S>(k),  for  Cu  K-edge  were  extracted 
from  the  CulnSer  spectrum  using  a Cu-Se  inter-atomic  distance  equal  to  2.424  A.  The 
obtained  structural  parameters  are  given  in  Table  7.4.  The  fitting  results  confirmed  that  the 
Cu-Se  inter-atomic  distance  is  conserved  between  CulnSe2  and  CulmSej.  The  nearest 
neighbor  coordination  number  for  Cu  is  fixed  at  4 for  the  least  square  fitting  in  both 
phases  A schematic  drawing  of  the  Cu  first  nearest  neighbor  local  structure  is  given  in 
figure  7. 13. 
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Figure  7.10  Cu-K  edge  (a)  Icyfk)  EXAFS  spectra,  and  (b)  magnitudes  of  the  Fourier 
transform  for  CubbSe* 
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Figure  7.1 1.  Fourier  filtered  EXAFS  on  the  Cu  K edge  for  the  first  shell  peaks  in  CulnSc, 
and  CuInjSej. 
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Figure  7 12  Least-square  fits  to  the  Fourier-filtered  first  shell  Cu  A'-edge  /(-weighted  x(*) 
data  of  CuItiiSej.  (The  data  for  CulmSe-  and  reference  CulnSe?  were  filtered  over  an 
identical  range  in  r space,  from  1.58  to  2 58  A.) 


Table  7 4 EXAFS  analysis  results 
N:  4 +-  0.2 
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Figure  7. 1 3 A schematic  drawing  of  the  Cu  first  nearest  neighbor  local  structure. 
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7 .5.3.2  In  local  structure 

The  In-K  edge  FT-R  space  data  for  CulniSei  are  given  in  Figure  7. 14.  The  first 
FT-peak  arises  at  2.26  A and  is  attributed  to  the  first  nearest  neighbors  (Se).  The  isolated 
first  shell  spectra  for  CulnSe2  and  Culn.Se-  are  given  in  Figure  7.15.  It  can  be  clearly  seen 
that  the  two  spectra  are  almost  identical.  This  indicates  the  first  nearest  neighbor 
environment  (i.e.,  bond  length,  dt^s*  and  number  of  nearest  neighbors,  N)  around  In 
atoms  in  these  two  semiconductor  alloys  are  identical.  The  backscattering  amplitude,  F,(k), 
and  phase  shift,  Sfli).  required  for  the  analysis  of  the  In-K  edge  CuInjSe;  spectrum  were 
extracted  from  the  In  K-edge  CuInSejEXAFS  spectrum. 

The  least  square  fits  of  the  ^-weighted  data  at  298  K are  shown  in  Figure  7. 16. 
The  estimated  structural  parameters  arc  given  in  Table  7.5.  The  fitting  results  confirmed 
that  the  In-Se  inter-atomic  distance  is  conserved  between  CuInSei  and  CuIntSe*.  The 
nearest  neighbor  coordination  number  is  fixed  at  4 for  the  least  square  fitting.  A schematic 
drawing  of  the  In  first-nearest-neighbor  local  structure  is  given  in  Figure  7 1 7. 
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Figure  7 14  Magnitudes  of  the  Fourier  transform  of  the  In-/f  edge  EXAFS,  ^x(*)  for 
CuInjScj. 
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Figure  7 15  Fourier  filtered  EXAFS  on  the  In  K edge  for  the  first  shell  peaks  in  CulnSej 
and  CulnrSes. 


Figure  7 17  Least-square  fit  to  the  Fourier-filtered  first  shell  In  A'-edge  ^'-weighted  z(*) 
data  of  CuIniSei  The  data  for  CuIniSe;  and  reference  CuinSe,  were  filtered  over  an 
identical  range  in  r space,  from  17  to  2.7  A. 


Table  7 5 EXAFS  analysis  results 
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Figure  7. 1 8 A schematic  drawing  of  the  In  first-n 
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7. 5. 3. 3 Se  local  structure 


The  Se-K  edge  raw  X-ray  absorption  data,  k'xik)  spectra,  and  FT-R  space  data  for 
CuIn.Se,  are  shown  in  Figure  7. 19.  Unlike  the  Cu-K  edge  data,  the  first  peak  shows 
interference  because  the  first  shell  consists  of  both  Cu  and  In  atoms.  The  Se-center 
tetrahedron  of  chalcopyrite  CuInSei  is  given  in  Figure  7.23(a),  which  consists  of  two  Cu 
and  two  in  atoms  around  the  central  Se  atom.  CulnrSe*  belongs  to  defect-tetrahedral 
structure  (Par  64],  which  can  be  described  as  normal  tetrahedral  structure  with  20% 
unoccupied  tetrahedral  sites.  The  first  shell  Se-K  edge  EXAFS  data  should  reveal  the 
number  of  nearest  neighbors  in  the  ratio  of  Cu:  In:  □ = 1 : 3:  1,  if  CulmSe.  did  belong  to 
defect-tetrahedral  structural  family. 

The  isolated  first  shell  spectra  for  CuInSei  and  Culn.iScj  are  given  in  Figure  7.19. 
It  is  clear  that  the  two  spectra  are  very  different.  This  indicates  the  first  nearest  neighbor 
environment  around  Se  atoms  in  these  two  semiconductor  alloys  is  different.  The  ab  uilio 
multiple-scattering  code  FEFF7  [Res  97]  was  used  to  calculate  F/k),  S>(k),  and  X(k)  to 
obtain  quantitative  estimates  of  the  number  of  nearest  neighbors  and  their  distances.  The 
isolated  first  shell  spectrum  for  CulnScj  and  the  least  squares  fit  are  given  in  Figure  7.20. 
The  fitting  was  performed  by  assuming  the  number  of  nearest  neighbors  (NN),  Cu=2, 
ln=2.  and  their  distances  do*se=  2.424  A,  and  di„.^, : 2.598  A.  Fitted  results  showed  the 
Se-K  edge  data  are  consistent  with  the  Cu  and  In  K-edge.  The  resulting  reduction  factor, 
So,  was  retained  and  used  in  the  analysis  of  Se  K-edge  spectrum  for  CuInjScs.  The  fitting 
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was  performed  using  NN  Cu=0.8,  NN  In=2.4.  dc»s.=  2.424  A,  dh.s.=  2.598  A.  The  FT- 
filler  first  shell  spectrum  for  CulmSer  with  the  best  least  square  lit  is  shown  in  Figure 
7.21  The  fitting  results  confirm  that  CulmSe*  does  belong  to  defect-tetrahedral  structure, 
which  is  characterized  by  a vacant  tetrahedral  site  in  the  Se-centered  tetrahedron  The  data 
are  consistent  with  the  ideas  proposed  by  [Zha  98J,  which  suggested  that  the  structure  of 
CuIn.Sc,  consists  of  three  types  of  local  tetrahedral  cationic  clusters  around  each  Se: 

3+Cu+2In  (k— 7),  2Cu+21n  (k=8),  and  0+3In  (k=9).  (□  denotes  a vacancy  and  k denotes 
the  sum  of  the  cation  valence  electron.  A schematic  drawing  of  these  local  structure 
models  is  given  in  Figure  7.22.  How  these  different  types  of  tetrahedron  arrange 
themselves  and  form  the  long-range  order  structure  (e.g.,  symmetry)  need  lurther  study  by 
other  techniques. 


210 


Figure  7,  IS  Se-K  edge  for  CuInjSc,  (a)  raw  X-ray  absorplion  data  (b)  EXAFS 
spectra,  and  (c)  magnitudes  of  the  Fourier  transform 


Figure  7 19  Fourier  filtered  EXAFS  on  the  Se-K  edge  for  the  first  shell  peaks  in  CulnSe. 
and  CulmSe. 


Figure  7.20  Least-square  fits  to  the  Fourier-filtered  first  shell  Se-K  edge  ^'-weighted  *(*) 
data  ofCulnSe2. 


Table  7.6  EXAFS  analysis  results  for 
CulnScj  Se-K  edge  data. 
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Figure  7 21  Least-square  fits  to  the  Fourier-filtered  first  shell  Se-K  edge  A1- weighted  /(A) 
dataofCuInjSes. 
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Figure  7.22  Schematic  drawing  of  the  three  types  of  local  tetrahedral  cation 
around  each  Se:  (a)  2Cu+2ln  (k=8),  (b)  [H-Cu+21n  (k=7).  and  (c)  0+3In  (k=9) 


To  further  investigate  this  problem,  a collaborative  research  effort  was  established 
with  the  solid  state  theory  group  at  the  National  Renewable  Energy  Laboratory  (NREL). 
The  group  at  NREL  has  calculated  the  Cu-Se  and  In-Se  bond  lengths  in  stoichiometric 
CulnSe;  and  Cu-poor  Culn.Se.  and  CuIniSc*  ordered  vacancy  compounds  using  first- 
principles,  self-consistent  electronic  structure  method.  The  total  energy  calculation  was 
performed  within  the  local  density  approximation  (LDA)  [Cep  81]  as  implemented  by  the 
general  potential  linearized  augmented  plane  wave  (LAPW)  method  [Wei  85] 

The  calculated  lattice  constants  a = 5.768  A,  t]  = 1.008  and  x = 0.217  for 
chalcopyrite  CulnSej  are  in  good  agreement  with  the  experimental  values  of  a = 5.784  A , 
tl  = 1 004  and  and  x = 0.224  [Spi  74],  To  minimize  the  total  energy,  our  calculation 
considers  only  structures  that  have  minimal  deviations  from  the  octet  rule,  i.e.,  in  these 
structures  the  sum  (denoted  as  k)  of  the  valence  electrons  of  the  four  cations  (including 
the  possibility  of  a vacancy  I ) that  surrounds  each  Se  atom  should  be  close  to  8.  For 
CulnSc;,  only  k=8  clusters  exist,  that  is  each  Se  is  surrounded  by  2Cu  + 21n.  Some  of  the 
crystal  structures  that  satisfy  this  condition  are  show  in  Figure  7 23.  Their  recent  study, 
[Wei  99]  show  that  all  the  k=8  structures  can  be  described  as  zinc-blende-based  polytypes, 
i.e.,  they  can  be  obtained  by  shifting  part  of  the  (001)  planes  of  the  chalcopyrite  (CH) 
structure  (Figure  7 23(a))  along  the  [1 10]  direction  by  r = (|}0)a.  For  example,  the 
CuAu  structure  (Figure  7.23(b)  with  c axis  in  the  [010]  direction)  can  be  obtained  from 
the  chalcopyrite  structure  by  periodically  shifting  two  out  of  every  four  planes,  while  the 
primitive  chalcopyrite  structure  (Figure  7.23(c))  can  be  obtained  by  shift  periodically  one 
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out  of  every  four  planes,  etc.  Notice  that  the  shift  along  the  [I  I0J  direction  by  r = ( r tO)  a 
is  equivalent  to  the  permutation  of  all  the  Cu  atom  with  all  the  In  atoms  in  that  plane.  For 
Culn.Se,  the  octet  rule  requirement  dictates  that  only  an  equal  number  of  k-7  (V  + Cu  + 
2ln)  and  k=9  (V  + 31n)  Se-centercd  clusters  can  exist.  These  structures  can  be  obtained 
from  CuInSe;  in  Figure  7.5  by  having  one  (2VC.  + Inc)  defect  pair  in  every  16-atom  (4 
CulnSei)  unit  cell  [Zha  98],  For  CuImSes,  it  can  be  treated  as  a mixture  of  CulnScj  and 
CulmSe,,  i.e„  they  contain  20%  of  the  k = 8 clusters  and  40%  each  of  the  k = 7 and  k = 9 

The  NREL  group’s  calculations  show  that  for  stoichiometric  CulnSei  the 
chalcopyrite  structure  (Figure  7.23(a))  has  the  lowest  energy,  in  agreement  with 
experiment  [Pal  67],  Thus  chalcopyrite  is  the  ground  state  structure  for  CulnSei.  The 
calculations  also  show  that  for  all  the  k = 8 polytypes  the  total  energy  difference  between 
any  of  the  two  polytype  stiuctures  is  very  small,  less  than  2 meV/atom.  This  result 
indicates  that  although  the  k=S  local  environment  is  well  conserved  below  the  order- 
disorder  transition  temperature  (Te  =1083  K)  [Pal  67],  the  long  range  chalcopyrite 
ordering  is  not  perfect  at  finite  temperature  [Wei  99].  Many  of  the  chalcopyrite  (001) 
layers  are  shifted  along  the  [1 10]  direction  where  Cu  occupies  the  nominally  In  site  and 
vise  versa.  The  actual  degree  of  the  long-range  chalcopyrite  order  of  CulnSe;  may  depend 
sensitively  on  the  growth  kinetics,  history  of  annealing,  and  the  configurational  entropies 
Indeed,  the  coexistence  of  chalcopyrite  and  CuAu-like  phases  in  CulnS;  has  been  reported 
[Sud  99],  This  imperfectness  of  long  range  chalcopyrite  ordering  will  have  large  effect  on 
the  XRD  analysis  of  the  crystal  stiuctures  (see  below)  Although  CuAu-like  phases  had 
not  been  reported  in  CulnSei,  The  observed  diffuse  scattering  in  the  electron  diffraction 
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studies  of  CulnSe,  by  Manolikas  and  Spyridclis  [Man  SI]  and  Kiely  et  al  [Kie  91]  clearly 
show  the  disruption  of  the  long  range  order  and  a strong  tendency  for  short  range 
ordering. 

This  result  is  consistent  with  the  classical  work  of  Mickkelsen  and  Boyce  [Mic  82] 
and  previous  theoretical  work  on  semiconductor  alloys  [Mar  84]  which  show  that  bond 
lengths  in  semiconductor  alloys  are  more  closer  to  their  ideal  values  than  to  the  alloy 
averaged  values.  It  is  also  consistent  with  the  fact  that  bond  length  conservation  is  more 
complete  for  alloys  with  large  ionic  character  (eg,,  II- VI  alloys  and  CulnSe,)  than  for 
more  covalent  compounds  (eg.,  III-V  alloys),  since  the  bond  bending  force  is  much 
weaker  for  ionic  alloys  [Mar  84]. 

Unlike  a full-scale  order-disorder  transition,  which  can  lead  to  a change  of  band 
gap  by  as  much  as  0.35  eV  [Wei  92],  the  formation  of  polytypes  leads  to  only  small 
changes  in  the  band  gap,  less  than  0.03  eV  (Wei  99],  The  band  gap  in  the  defect 
tetrahedral  compounds  is  larger  than  the  band  gap  for  CulnSe,  (Zha  98,  Sch  95],  The 
reason  for  this  is  as  follows:  In  the  defect  tetrahedral  compounds  some  of  the  Cu  atoms 
are  replaced  by  Cu  vacancies  Vc„  and  ltv„  antisites,  thus,  the  coupling  between  Se  p and 
Cur/ is  reduced  in  the  defect  tetrahedral  compounds  as  compared  to  CulnSe,.  This  weaker 
P d repulsion  lowers  the  Valence  Band  Maximum  (VBM)  of  the  ordered  vacancy 
compounds  and  increases  their  band  gaps 

Both  EXAFS  and  first  principle  calculation  results  show  that  the  first-nearest- 
neighbor  local  structures  (i.e.,  the  bond  length,  d,-„.s„  and  the  number  of  nearest 
neighbors,  N)  around  Cu  atoms  in  CulnSe,  is  not  a function  of  stoichiometry  nor  is 
vacancy  This  is  contrary  to  the  conclusions  of  Merino  M al.  [Mer  96]  from  X-ray  powder 
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diffraction  (XRD)  and  Rietveld  refinement.  The  reason  behind  this  discrepancy  is  found  in 
the  differences  between  these  two  characterization  techniques  XRD  is  not  suitable  for 
distinguishing  individual  bond  length  when  two  different  type  of  atoms  share  the  same 
crystallographic  site.  On  the  other  hand  XAFS  is  capable  of  measuring  local  structure 
around  specific  type  of  atoms,  including  the  bond  lengths. 

It  is  believe  that  the  large  dependence  of  the  Cu-Se  and  In-Se  bond  lengths  on  Cu 
occupation  number  reported  in  [Mer  96]  is  an  artifact  of  the  fitting  procedures.  Part  of  the 
difficulty  is  related  to  the  determination  of  the  Cu  occupation  number  [Mer  96]  used  three 
schemes  were  used  in  the  Rietveld  refinement  process  In  the  first  scheme,  the  Se  content 
Is  fixed  at  2,  the  In  content  at  1,  and  the  Cu  content  was  refined.  In  the  second,  only  Se  is 
fixed  at  2,  the  In  is  as  observed  by  EDAX,  and  the  Cu  is  refined  In  the  third,  both  In  and 
Se  are  from  EDAX  and  the  Cu  is  refined.  The  three  different  schemes  yielded  similar  Cu 
occupation  number,  despite  the  5 76%  difference  between  the  In  content  in  scheme  1 and 
scheme  2.  This  observation  indicates  that  the  refinement  is  not  very  sensitive  to  the 
occupancy  number.  For  example,  it  would  be  interesting  to  see  how  much  the  R factor 
changes  if  one  fixed  the  Cu  content  using  the  EDX  data.  Another  interesting  scheme  is 
allowing  all  three  occupancy  numbers  to  vary  in  the  refinement. 

Second,  as  show  by  Jaff  and  Zunger  [Jae  84]  and  more  recently  Wei  el  at.  [Wei 
98],  large  change  in  x(Se),  and  thus,  the  Cu-Se  and  In-Se  bond  lengths,  will  lead  to  a large 
change  in  the  band  gap.  A change  from  0 224  to  0.234  in  x(Se)  is  expected  to  cause  an 
increase  of  band  gap  by  0. 125  eV  [Wei  98],  This  large  band  gap  increase  is  not  supported 
by  the  measurement  of  Merino  et  al.  [Mer  98]  for  the  same  samples  In  fact,  some  of  the 
samples  with  large  derived  x(Se)  actually  have  smaller  band  gaps. 


It  is  likely  thai  the  samples  used  in  the  measurement  by  Merino  el  al.  show  some 
degree  of  partial  disordering  of  Cu  and  In  cations  or  the  formation  of  CuInSe,  polytypes  in 
short  range  order.  By  assuming  the  samples  have  perfect  chalcopyrite  ordering  in  their 
Rietveld  analysis,  some  of  the  In-Se  bond  is  treated  as  Cu-Se  bonds  and  vise  versa.  This 
will  overestimate  the  smaller  Cu-Se  bond  and  underestimate  the  longer  In-Se  bond,  thus 
give  a false  impression  that  J<Se)  increases  The  fact  that  formation  ofCulnSe,  polytypes 
in  the  samples  is  also  consistent  with  the  observed  small  band  gap  changes  in  the  samples, 
and  the  fact  that  sample  B I which  has  "largest"  .t(Se)  actually  has  smallest  band  gap 

In  addition  to  the  two  reasons  staled  above,  another  source  of  complication  could 
come  from  formation  of  the  CujSe  phase  in  the  samples,  especially  for  the  Cu-rich  samples 
used  in  [Mer  96].  It  is  well  known  that  the  powder  XRD  pattern  for  CujSe  and  CulnSei 
are  not  easy  to  resolve  if  a mixed  phase  was  examined  Simulated  powder  XRD  patterns 
for  a single-phase  chalcopyrite  CuInSe;  and  a mixed  phase  (10%  CujSe  and  90% 
CulnSe,)  are  shown  in  Figure  7.24(a)  and  (b).  respectively  The  only  difference  is  a small 
change  in  the  relative  peak  intensity  and  slight  shift  in  peak  position.  These  can 
not  be  easily  distinguished  from  the  preferred  orientation  effect  and  lattice  parameter 


Figure  7.23  Crystal  structures  of  CulnSe,  in  (a)  chalcopyrite,  (b)  CuAu-like,  and  (c) 
primitive  phases  For  clarity  the  Se  atoms  are  not  shown  The  lines  with  arrow  give  the 
primitive  unit  cell  lattice  vectors  for  these  structures, 
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Figure  7 24  (a)  Calculated  Cu-Ka  powder  XRD  pattern  for  chalcopyrite  CulnSe, 
structure  (Spi  74]  (b)  Calculated  Cu-Ka  powder  XRD  pattern  for  a mixed  chalcopyrite 
CulnSej  and  Cu.Se  [Hey  66]  sample  (10  % CujSe  and  90%  CulnSe;),  DBWS  [You  95] 
program  was  used  for  the  simulation. 
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[Mcr  96,  98]  extended  their  results  to  explain  the  spread  in  the  energy  gap  values 
reported  in  the  literature  based  on  Jaffc  and  Zunger's  model  [Jaf  84]  [Mer  96,  98] 
compared  their  experimental  Sd^.JSt  and  &,.„/&  values  with  that  derived  from 
[Mer  98]  and  find  good  agreement  between  them.  However,  the  value  of  Sd^.JSt  and 
K..*  ^ can  be  directly  derived  from  the  structure  model  and  given  as, 

= *(\  + Tf)n6t'a  (7  |4) 

+(l  + <r)/|6]1  '«  (7  15) 

A plot  of  Sd^JSc  and  Sd,._„lm  vs.  x (Se)  is  given  in  Figure  7.24  The  agreement 
between  the  values  of  and  «&/,„.„/&  in  [Mer  96]  and  [Jaf  84]  are  a direct 

consequence  of  using  the  same  structural  model  It  could  not  be  used  as  a measure  of  how 
well  the  experimental  data  agree  with  the  theoretical  calculations. 
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Figure  7 26  (a) SJC, and(b)  vs.  x (Se) (a=5.782  A,  ifl.005). 
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7.7  Conclusions 

XAFS  experiments  were  performed  to  study  the  local  structure  of  various  copper 


Fluorescence  Cu  K-edge  X-ray  Absorption  Fine  Structure  (XAFS)  measurements 
were  performed  on  Cu,ln,Se,  crystals  with  difterent  compositions  close  to  the 
stoichiometry  CulnSej.  The  result  indicates  the  local  structures  (i.e.  bond  length,  d(vse, 
and  number  of  nearest  neighbors,  N)  around  Cu  atoms  in  these  semiconductor  alloys  are 
the  same.  Contrary  to  other  X-ray  powder  diffraction  (XRD)  and  Rietveld  refinement 
work,  which  suggested  the  Cu-Se  and  In-Se  bond  length  values  depend  strongly  on  Cu 
occupation  number,  these  studies  indicate  the  environment  around  Cu  is  not  affected  by 
small  changes  in  chemical  composition.  Total  energy  calculations  were  also  performed, 
and  show  that  the  Cu-Se  bond  length  is  independent  of  the  Cu  occupation  number  while 
In-Se  bond  lengths  decrease  only  slightly  when  the  Cu  occupation  decreases.  The 
calculation  also  suggested  the  formation  of  other  polytype  structures  in  CulnSej  in 
addition  to  the  chalcopyrite  structure.  It  is  likely  the  samples  used  in  the  measurement  by 
[Mcr  96)  showed  some  degree  of  partial  disordering  of  Cu  and  In  cations  or  the  formation 
CulnSej  polytype  in  short  range. 

A detail  EXAFS  study  of  CulnSej  and  CuIn)Se>  on  Cu-K,  In-K,  and  Se-K  edges 
was  performed  It  was  found  the  Cu  and  In  first  nearest  neighbor  local  structures  in 
Culn.Ses  arc  almost  identical  to  those  in  CulnSej.  The  Se  first  nearest  neighbor  local 
structures,  however,  arc  quite  different  The  least-square  fitting  of  Se-K  edge  EXAFS 
spectra  indicated  CulnSej  consists  of  Se-centercd  tetrahedron  with  average  2 Cu  and  2 In 
as  nearest  neighbors  On  the  other  hand,  Culn,Ses  consists  of  Se-centered  tetrahedron 
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with  average  0.8  Cu  and  2.4  In  as  nearest  neighbors  This  result  directly  proved  Culn,Ses 
is  indeed  a defect  tetrahedral  structure. 


CHAPTER  8 

CONCLUSIONS  AND  FUTURE  WORK 
This  dissertation  has  focused  on  developing  a detailed  thermodynamic  description 
of  the  Cu-ln-Sc  system  for  belter  understanding  of  the  processing  chemistry  of  CulnSe; 
absorber  layers.  Knowledge  of  thermochemistry  and  phase  diagram  was  applied  to 
reaction  engineering  of  the  precursor  structure  for  RTP  synthesis  of  CuIilSe:  absorber 
layers.  Another  focus  of  this  work  was  a detailed  study  of  the  local  structure  of  copper 
indium  sclenide  semiconductor  alloys  by  EXAFS. 

8. 1 Conclusions 

8.  I I Thermodynamic  and  phase  diagram  assessment 

In  Chapter  3.  the  assessment  of  the  thermodynamic  properties  of  elemental  Sefc, 
hexagonal).  Se(l).  and  Se„(g).  (n=I  to  8)  resulted  in  suggested  values  of  AH'rjM.isK, 
S'Wuk.  and  Cp(T).  The  thermodynamic  (unctions  of  each  gas  phase  species  was  assessed 
from  the  available  literature  data  or  calculated  from  molecular  constants  assuming 
harmonic  oscillator,  rigid  rotor,  and  ideal  gas  behavior.  The  values  of  AH°un  ,!K  for 
Ses(g)  and  Ses(g)  were  optimized  using  the  vapor  pressure  data  in  equilibrium  with  Se(l) 
and  Se(c,  hexagonal),  respectively  Gas  phase  nonidcalitics  were  taken  into  account  by 
using  the  Peng-Robinson  equation  of  sate.  The  result  of  this  assessment  are  expressions 
for  the  Gibbs  energy  of  each  of  the  ten  species  Se(c,  hexagonal),  Se(l),  and  Sty(g).  {n=l 
to  8}  as  a function  of  temperature.  These  expressions  were 


then  used  to  calculate 
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equilibria  in  ihe  Se  system.  The  vapor  pressures  obtained  from  equilibrium  calculations 
agree  well  with  the  experimental  results.  The  equilibrium  calculations  showed  that  the 
predominant  vapor  species  at  a temperature  above  approximately  1000  K is  the  dimer, 
while  at  a lower  temperature  I 600  to  1000  K),  Se.  is  the  dominant  species.  Below  600 
K the  major  species  is  Se«(g)  During  growth  using  a Se  evaporation  source,  a change  in 
the  deposition  chemistry  is  expected  as  the  temperature  is  changed  through  one  of  these 
transitions.  The  results  were  used  to  estimate  the  Langmuir  evaporation  flux  for  each 
species  and  is  useful  in  characterizing  the  performance  of  an  evaporative  source. 

In  Chapter  4,  a complete  database  for  ihe  Cu-Se  system  was  assembled  and  a 
critical  assessment  was  performed.  A self-consistent  set  of  phase  diagram  and 
thermochemical  data  was  selected.  The  selected  database  was  combined  with  the  unary 
data  for  copper  provided  by  SGTE  and  the  assessed  Se  data  in  Chapter  3.  The  association 
model  was  used  to  describe  the  liquid  phase.  A three-sublatticc  compound  energy 
formalism  was  used  for  the  Cu;.*Se  phase.  The  Redlich-Kister  polynomial  was  used  to 
represent  the  Gibbs  energy  of  the  terminal  solid  solution  fcc(Cu).  The  remaining 
intermediate  solid  phases  (CujSea,  CuSe,  CuSe.-)  were  modeled  as  line  compounds. 

The  phase  diagram  and  thermodynamic  properties  of  this  system  have  been 
calculated  using  the  optimized  parameters.  The  calculated  phase  diagram  agrees  well  with 
the  data  used  in  the  optimization.  In  the  calculated  phase  diagram,  the  temperatures  and 
compositions  of  the  Se-rich  liquid  miscibility  gap  have  been  predicted.  The  site  tractions  of 
vacancy  in  the  fl-Cu;„Se  phase  were  calculated.  In  the  p-Cu;.,Se  single-phase  region,  the 
Cu  vacancy  concentration  in  the  first  sublattice,  , increases  with  an  increase  of 
selenium  concentration;  whereas  the  Se  vacancy  concentration,  .iy, , decreases  with  an 
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increase  in  the  selenium  concentration.  The  calculated  Cu  vacancy  concentration  in  the 
first  sublaltice  is  between  0.008  to  0.271.  whereas  the  Se  vacancy  site  fraction  is  between 
2.07xl0',T  to  2.23x  lO-1. 

The  calculated  enthalpies  of  formation  in  the  present  work  are  between  the  data 
measured  by  [Gat  56]  and  those  evaluated  by  [Mil  74],  The  calculated  standard  entropies 
agree  with  those  evaluated  by  [Mil  74]  except  that  of  CuSej,  which  lies  in  between  the 
value  of  [Rau  70]  and  that  evaluated  by  [Mil  74]  The  calculated  heat  capacity  values 
agree  well  with  the  measured  ones  [Kub  73].  The  transition  of  a-CujSe-»p-Cu;Se  is 
assessed  to  occur  at  (306  K with  an  enthalpy  change  of  6830  J/mole,  which  is  in  good 
agreement  with  that  measured  by  [Kub  73],  The  calculated  enthalpy  of  the  transition  of  a- 
CuSe— sP-CuSe  at  324  K is  1381  J/mole.  which  is  in  good  agreement  with  that  measured 
by  [Hey  66].  As  no  experimental  heat  capacity  data  for  CuSe  are  available,  those  data 
were  optimized  from  the  phase  diagram  data  after  assuming  that  the  heat  capacities  of  a- 
CuSe  and  p-CuSe  are  identical. 

In  a similar  fashion,  a complete  database  for  the  ln-Se  system  was  assembled  and 
presented  in  Chapter  5.  The  phase  diagram  and  thermochemistty  of  the  In-Sc  system  has 
been  experimentally  studied  by  several  authors  and  the  results  for  selected  properties  are 
inconsistent  Uarly  studies  reported  four  intermediate  compounds  (InjSe,  InSe.  InjSej,  and 
ln?Ser).  while  InSe  and  In-Se.  were  identified  as  congruent  melting  compounds.  A 
miscibility  gap  on  the  In-rich  side  was  also  suggested.  Later  work  reported  the  existence 
of  the  four  compounds  ImSej  (instead  of  ImSe).  InSe.  ln,Se,  (instead  of  IniSe*),  and 
In;Se>.  with  only  Ift.Se>  melting  congruently  More  recent  work  by  (God  98]  reported 
IoiSeii  and  ImSe?  as  stable  phases  and  [i-IftSe.  as  stable  at  59.6  Se  at.%  in  contrary  to 


previous  study.  A self-consistent  set  of  phase  diagram  and  thermochemical  data  was 
selected  The  selected  database  was  combined  with  the  unary  data  for  indium  provided  by 
SGTE  and  the  assessed  Se  data  presented  in  Chapter  3 Nine  intermediate  solid  phases 
(IitjSc'..  InSe,  In.  Sc*.  I El  . S e I ■ . In.-Se:.  a-IniSe3,  P-  ln;Se',.  y-ln:Se,.  and  S-lnaSes)  were 
modeled  as  line  compounds.  The  association  model  [Som  82]  was  used  to  describe  the 
liquid  phase.  Missing  heat  capacity  C I ll'Se-7.  ln?Sen)  and  entropy  data  (InsSe*  ImSeu) 
were  first  estimated  by  UnSI's  [Kub  77]  rule  and  Latimer's  [Lat  75]  rule,  respectively,  and 
a re-interpreting  of  the  literature  heat  capacity  and  entropy  data  for  a-IivSe,.  p.  [n.Se.  y- 
ln_.Sc,,  and  6-In;Sc,  was  done  based  on  the  new  phase  diagram.  These  data  were  then 
assessed  through  optimization  using  ThermoCalc. 

It  was  found  that  the  reported  -ALCrtw  iik  value  for  liliSe,  [Mus74,  Cha  93]  is  too 
small  to  be  consistent  with  the  phase  diagram  data.  Using  the  reported  value  IiuSe? 
decomposes  at  a much  lower  temperature  (-  600K)  and  does  not  reach  the  peritectic 
temperature  (793K).  By  coupling  the  thermochemical  and  the  phase  diagram  data, 
optimized  G-Hsek  expressions  for  the  solid  compounds  and  the  liquid  phase  are  obtained 
together  with  the  major  vapor  phase  species  (Sen  (n~  I to  8),  In,  InSe,  In;Se,  IniSe?)  in  In- 
Se system.  The  T-x  and  P-T  projection  of  the  P-T-x  diagram  were  calculated.  The 
calculated  P-T  diagram  is  in  reasonable  agreement  with  the  reported  data  [Gri  72,  Sri  87, 
Bar  88,  Pas  89]  and  consistent  with  the  recent  T-x  phase  diagram  [God  98].  It  is  also 
suggested  that  the  experimental  data  of  [Sri  87,  Bar  88,  Pas  89]  represented  the  vapor 
pressure  in  equilibrium  with  the  fJ-ln2Ses,  fl-In;Se,  and  liquid  2,  and  liquid  2 phase 
domains,  rather  than  that  in  equilibrium  with  a-In-Se,  as  the  investigators  presumed. 


8 1.2  CuInSe;  thin  film  formation  by  RTP 

In  Chapter  6.  knowledge  of  thermochemistry  and  phase  diagram  was  applied  to 
engineering  of  alternative  reaction  pathway  for  CIS  thin  films.  A novel  bi-layer  precursor 
structure  of  a co-deposiled  In-Se  (50  to  55  Se  at.%)  top  layer  and  a Cu-Se  (50  to  55  Se. 
at.%)  bottom  layer  was  proposed  It  was  postulated  that  this  precursor  structure  should 
lead  to  CIS  formation  through  a liquid-solid  reaction  in  the  processing  temperature 
window  524°C  to  6I1”C.  Two  precursor  structures  were  deposited  in  a migration- 
enhanced  deposition  reactor.  The  first  structure  is  the  three  layered  structure 
IruSej/CuInSeyCujSe.  After  RTP,  it  was  found  the  precursor  film  underwent  a liquid- 
solid  reaction  and  formed  large  grained  CulnSci  films  at  a substrate  temperature  of  SOOT, 
as  suggested  by  SEM  and  XRD  analysis  It  was  believed  that  the  melting  of  top  layer 
liuSe*  phase  occur  through  a monotetic  reaction.  The  second  structure  consists  of  only 
two  alloy  layers,  InSe/CuSe.  This  precursor  structure  lead  to  the  formation  of  large 
grained  CIS(~5pm)  films  at  a substrate  temperature  of  SOOT,  consistent  with  the  pathway 
suggested  by  the  phase  diagram.  It  was  found  that  the  ramp  rate  has  a dramatic  effect  on 
the  film  quality  for  the  second  structure.  A slow  ramp  rate  (l°C/s)  resulted  in  bubble 
formation  and  de-lamination  of  the  film.  A higher  ramp  rate  (lOT/s)  produced  CuInSe2 
films  with  better  crystallinity  and  no  bubbke  formation  or  de-lamination  problems.  It  is 
believed  that  bubble  formed  by  evaporation  of  Se,  vapor  from  the  CuSe  layer.  This  is 
supported  by  RTP  experiments  on  CuSe  reference  films  under  the  same  conditions.  The 
CuSe  film  has  a lower  Se  content  and  transformed  to  Cu2.,Se.  With  InSc  layer  on  top, 
however,  the  overall  Se  content  did  not  change  significantly  (<0  5%)  according  to  the  ICP 


chemical  analysis.  The  presence  of  the  InSe  layer  provided  a barrier  lo  Se  evaporation  loss 
in  the  film. 

S.  1.3  XAFS  investigations  of  copper  indium  selenide  semiconductor  alloys 

XAFS  experiments  were  performed  to  study  the  local  structure  of  various  copper 
indium  selenide  semiconductor  alloys 

Fluorescence  Cu  K-edge  X-ray  Absorption  Fine  Structure  (XAFS)  measurements 
were  performed  on  Cu,InySe,  crystals  with  different  compositions  close  to  the 
stoichiometry  CulnSe2  and  CulnsSeg.  The  EXAFS  data  were  analyzed  using  the  same 
range  of  k-window  for  Fourier  transform  and  same  R-window  for  Fourier  filtering  of  the 
first  shell  EXAFS  oscillations  The  resulting  first  shell  Cue  k’yfk)  spectra  from  different 
stoichiometries  were  almost  identical.  This  result  indicates  the  local  structures  (i.e.  bond 
length,  deva.  and  number  of  nearest  neighbors.  N)  around  Cu  atoms  in  these 
semiconductor  alloys  arc  the  same.  Contrary  lo  other  X-ray  powder  diffraction  (XRD)  and 
Rietveld  refinement  work,  which  suggested  the  Cu-Se  and  In-Se  bond  length  values 
depend  strongly  on  Cu  occupation  number,  these  sludies  indicate  the  environment  around 
Cu  is  not  affected  by  small  changes  in  chemical  composition. 

Total  energy  calculations  were  also  performed,  and  show  that  the  Cu-Se  bond 
length  is  independent  of  the  Cu  occupation  number  while  In-Se  bond  lengths  decrease  only 
slightly  when  the  Cu  occupation  decreases.  The  calculation  also  suggested  the  formation 
of  other  polytype  structures  in  CuInSe.  in  addition  to  the  chalcopyritc  structure.  It  is  likely 
the  samples  used  in  the  measurement  by  [Mer  96]  showed  a degree  of  partial  disordering 
of  Cu  and  In  cations  or  the  formation  CulnSe;  polytypes  in  the  short  range.  The 
discrepancy  between  the  results  of  the  present  work  and  that  of  (Mer  96]  could  be 
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explained  by  the  different  nature  of  the  characterization  techniques:  XRD  is  not  suitable 
for  distinguishing  individual  bond  lengths  when  two  different  types  of  atoms  share  the 
same  crystallographic  site,  while  XAFS  is  capable  of  measuring  local  structure  around 
specific  types  of  atoms,  including  the  bond  lengths.  The  large  dependence  of  the  Cu-Se 
and  In-Se  bond  lengths  on  Cu  occupation  number  reported  in  the  literature  may  be  an 
artifact  of  the  fitting  procedure  used  A further  implication  of  this  result  suggests  the 
existence  of  short  range  order  CuInSe<  polytypes  in  CuInSe;  thin  film  photovoltaics. 

A detail  EXAFS  study  of  CuInSe;  and  CuIn.Se.  on  Cu-K,  In-K,  and  Se-K  edges 
was  performed.  It  was  found  that  the  Cu  and  In  first  nearest  neighbor  local  structures  in 
CuInjSej  arc  almost  identical  to  those  in  CuInSe;.  However,  the  Se  first  nearest  neighbor 
local  structures  are  quite  different  The  least-square  fitting  of  the  Sc-K  edge  EXAFS 
spectra  indicated  CuInSe;  consists  of  Se-centered  tetrahedron  with  an  average  2 Cu  and  2 
In  nearest  neighbors.  On  the  other  hand.  CulmSei  consists  of  Se-centered  tetrahedron 
with  an  average  of  0.8  Cu  and  2.4  In  as  nearest  neighbors  This  result  directly  proves 
CulnjSej  is  indeed  a defect  tetrahedral  structure  and  is  consistent  with  first-principle 
calculation  [Zha  98].  The  calculation  suggests  there  are  three  types  of  local  tetrahedral 
cationic  clusters  around  each  Se:  V+Cu+21n  (4=7),  2Cu+21n  (4=8).  and  V+3In  (4=9). 
where  4 denotes  the  sum  of  valence  electrons  for  cations  CuInSe;  consists  of  100%  lt=8 
clusters  and  CuIniSes  consists  of  20%  4=8  and  40%  4=7  and  40%  4=9  clusters 
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8.2  Future  work 

8.2.1  Thermodynamic  and  phase  diagram  assessment 

In  this  work,  a critical  thermodynamic  and  phase  diagram  assessment  of  the  Se 
unary.  Cu-Se.  and  In-Se  binary  systems  were  completed-  It  is  recommended  to  utilize  this 
database  for  analysis  of  binary  flux  calculation  during  the  binary  precursor  deposition. 

More  experimental  data  are  needed  for  the  Cu_'Se-lrt:Se<  pseudo-binary  section.  It 
is  recommended  to  perform  experiments  to  determine  the  important  phase  boundaries  and 
the  missing  energy  data  such  as  the  heat  capacity  data  for  most  of  the  ternary  phases  A 
complete  assessment  of  the  Cu-In-Se  ternary  system  will  be  the  next  step  after  these 
critical  data  are  available. 

8.2.2  CuInSe?  thin  film  formation  by  RTP 

In  this  work,  a novel  precursor  structure  for  RTP  synthesis  of  CuInSe,  film  was 
proposed  and  subjected  to  RTP  treatment.  The  experimental  results  showed  the  proposed 
precursor  structure  is  extremely  promising.  Large  grain  CIS  films  were  produced  without 
de-lamination  problems  in  a very  short  time  (-70  s),  It  is  recommended  to  optimize  the 
processing  conditions  based  on  the  solar  cell  performance,  by  a DOE  (Design  of 
Experiments)  methodology. 

It  is  also  recommended  to  study  the  effect  of  Na.  Preliminary  results  showed  the 
bubble  formation  is  much  more  prominent  when  Na  is  added  to  the  precursor.  The 
mechanism,  however,  is  not  clear. 

8.2.3  XAFS  investigations  of  copper  indium  selenide  semiconductor  alloys 

The  XAFS  results  of  this  work  suggest  the  coexistence  of  a polytype  structure  in 
CuInSe?  materials  in  addition  to  chalcopyrite  ordering.  Interestingly,  epitaxial  CuAu- 
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ordering  CuInSej  siruciure  was  produced  in  our  laboratory  [Sta  99]  for  the  first  time-  A 
theoretical  study  of  the  impact  of  the  polytypes  on  photovoltaic  devices  is  recommended. 

A detail  study  of  the  local  structure  of  CuImSej  was  reported  in  this  work.  A 
complete  structure  model,  however,  is  still  needed.  It  is  recommended  to  perform 
Convergent  Beam  Electron  Diflraction  (CBED)  in  combination  with  X-ray  difiraction 
experiments  to  resolve  the  OVC  issue. 
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